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Abstract 



We discuss exclusive radiative decays in QCD factorization within the Standard Model. 
In particular, we consider the decays B — > V7, with a vector meson K* or p in the 
final state, and the double radiative modes B° s — > 77 and B® — * 77. At quark level, 
all these decays are governed by the flavour- changing neutral-current b — > 57 or 6 — > 
transitions, which appear at the one-loop level in the Standard Model. Such processes 
allow us to study CP violation and the interplay of strong and electroweak interactions, to 
determine parameters of the CKM matrix, and to search for New Physics. The exclusive 
decays are experimentally better accessible, but pose more problems for the theoretical 
analysis. The heavy-quark limit m& Aqcd, however, allows to systematically separate 
perturbatively calculable hard scattering kernels from nonperturbative form factors and 
universal light-cone distribution amplitudes. 

The main results of this work are the following: 

• We apply QCD factorization methods based on the heavy-quark limit to the 
hadronic matrix elements of the exclusive radiative decays B — > V7 and B — > 77. 
A power counting in KqcD/ m b implies a hierarchy among the possible transition 
mechanisms and allows to identify leading and subleading contributions. In partic- 
ular, effects from quark loops are calculable in terms of perturbative hard-scattering 
functions and universal meson light-cone distribution amplitudes rather than being 
generic, uncalculable long-distance contributions. Our approach is model indepen- 
dent and similar in spirit to the treatment of hadronic matrix elements in two-body 
non-leptonic B decays formulated by Beneke, Buchalla, Neubert, and Sachrajda. 

• For the B — > V'-f decays we evaluate the leading AqcD/^b contributions complete 
to next-to-leading order in QCD. We adopt existing results for the hard-vertex 
corrections and calculate in addition hard-spectator corrections, including also QCD 
penguin operators. 

• Weak annihilation topologies in B — > V7 are shown to be power suppressed. We 
prove to one-loop order that they are nevertheless calculable within QCD factor- 
ization. Because they are numerically enhanced we include the 0(a®) annihilation 
contributions of current-current and QCD-penguin operators in our analysis. 

• The double radiative B — > 77 decays are analyzed with leading-logarithmic accu- 
racy. In the heavy quark limit the dominant contribution at leading power comes 
from a single diagram. The contributions from one-particle irreducible diagrams 
are power suppressed but still calculable within QCD factorization. We use these 
corrections, including QCD penguins, to estimate CP asymmetries in B — > 77 and 
so-called long-distance contributions in B and D — > 77. 

• We predict branching ratios, CP and isospin asymmetries, and estimate U-spin 
breaking effects for B — > K*^ and B — > p'-f. For the B — > 77 decays we give numer- 
ical results for branching ratios and CP asymmetries. Varying the individual input 
parameters we estimate the error of our predictions. The dominant uncertainty 
comes from the poorly known nonperturbative input parameters. 
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Chapter 1 
Introduction 



Elementary particle physics represents man's effort to answer the basic question: "What 
is the World Made of?" In search of the fundamental building blocks of matter physicists 
penetrated to smaller and smaller constituents, which later turned out to be divisible. The 
primary matter of Anaximenes of Miletus, the periodic table of Mendeleev, Rutherford's 
alpha scattering experiments, the detection of the neutron by Chadwick and that of the 
positron by Anderson, the discovery of nuclear fission by Hahn and Meitner, or Reines 
finding the neutrino were just some of the milestones on this way. What underlies our 
current theoretical understanding of nature is quantum field theory in combination with a 
gauge principle. Electromagnetism, weak and strong nuclear forces, and their interaction 
with quarks and leptons are described by the Standard Model of particle physics 
Combined with general relativity, this theory is so far consistent with virtually all physics 
down to the scales probed by particle accelerators, roughly 10" 16 cm, and also passes a 
variety of indirect tests that probe even shorter distances. 

In spite of its impressive successes, the Standard Model is believed to be not com- 
plete. For a really final theory it is too arbitrary, especially considering the large number 
of sometimes even "unnatural" parameters in the Lagrangian. Examples for such pa- 
rameters, that are largely different from what one naively expects them to be, are the 
weak scale compared with the Planck scale or the small value of the strong CP-violation 
parameter #qcd- Questions like: "Why are there three particle generations?", "Why is 
the gauge structure with the assignment of charges as it is?" , or "What is the origin of the 
mass spectrum?" demand an answer by a really fundamental theory, but the Standard 
Model gives no replies. Furthermore, the union of gravity with quantum theory yields a 
nonrenormalizable quantum field theory, indicating that New Physics should show up at 
very high energies. 

The ideas of grand unification, extra dimensions, or supersymmetry were put forward 
to find a more complete theory. But applying these ideas has not yet led to theories 
that are substantially simpler or less arbitrary than the Standard Model. To date, string 
theory ||, the relativistic quantum theory of one-dimensional objects, is a promising, and 
so far the only, candidate for such a "Theory of Everything" . 

Within this thesis we work exclusively in the Standard Model, where particularly 
in the flavour sector beside the aforementioned problems some unsolved issues remain. 
Among them are the mechanism of electroweak symmetry breaking and the hierarchy 
problem that comes along with it, the generation of fermion masses, quark mixing, or 
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the violation of the discrete symmetries C, P, CP, and T. Especially CP violation is of 
particular interest as it is one of the three vital ingredients to generate a cosmological 
matter-antimatter asymmetry [|J. In 1964, James Cronin, Val Fitch, and collaborators 
discovered that the decays of neutral kaons do not respect CP symmetry pL Only recently 
CP violation was established undoubtedly also in the neutral 5-meson system || . 

We will deal with the bottom quark system, which is an ideal laboratory for studying 
flavour physics. The history of B physics started 1977 with the "observation of a dimuon 
resonance at 9.5 GeV in 400 GeV proton-nucleon collisions" at Fermilab M. It was bap- 
tized the "T resonance" and its quark content is bb. With the start of BaBar [|J and Belle 



in 1999, dedicated B factories add a wealth of data to the results of CLEO 10 and the 



CERN |TT[ and Fermilab |T^] experiments. The upcoming 5-physics experiments at the 



Tevatron Run II fl3| and LHC |T4| will bring us ever closer to the main goal of B physics, 



which is a precision study of the flavour sector with its phenomenon of CP violation to 
pass the buck of being the experimentally least constrained part of the Standard Model. 
This is not only to pin down the parameters of the Standard Model, but in particular 
to reveal New Physics effects via deviations of measured observables from the Standard 
Model expectation. Such an indirect search for New Physics is complementary to the 
direct search at particle accelerators. It invites both experimenters and theoreticians to 
work with precision. We need accurate and reliable measurements and calculations. The 
calculational challenge we will meet for this thesis are exclusive radiative decays of B 
mesons. 

But why investigate B meson decays? Due to confinement quarks appear in nature not 
separately, but have to be bound into colourless hadrons. Considering constituent quarks 
only, the simplest possible of such objects consists of a quark and an antiquark only and 
is called a meson. The bound states with a b quark and a d or u antiquark are referred 
to as the B° and B~ meson, respectively. Those B mesons containing an s or c quark 
are denoted B s and B c , respectively. So "meson" decays because mesons are the simplest 
hadrons. But why of all mesons U B" mesons? Apart from the T resonances, the B mesons 
are the heaviest mesons, because the top quark decays before it can hadronize. The fact 
that B mesons are heavy has two weighty consequences: B decays show an extremely rich 
phenomenology and theoretical techniques using an expansion in the heavy mass allow 
for model-independent predictions. The rich phenomenology is based on the one hand on 
the large available phase space allowing for a plethora of possible final states and on the 
other hand on the possibility for large CP-violating asymmetries in B decays. The latter 
feature is in contrast to the Standard Model expectations for decays of K and D mesons. 
In D decays only the comparably light d, s, and b quarks can enter internal loops which 
leads to a strong GIM suppression of CP-violating phenomena. CP violation in K — > tttt 
is small due to flavour suppression and not because the CP violating phase itself is small. 
Actually, the sin 2(3 measurement in B — > J/ipKs showed that the CP- violating phase 
is large. Furthermore, especially the B — > J/ipKs decay mode is theoretically extremely 
clean as opposed to the large theoretical uncertainties in the kaon system. The pattern 
of CP violation in the K and B system just represents the hierarchy of the CKM matrix. 
The B meson system offers an excellent laboratory to quantitatively test the CP-violating 
sector of the Standard Model, determine fundamental parameters, study the interplay of 
strong and electroweak interactions, or search for New Physics. 

We will concentrate on a subgroup of B decays: exclusive radiative B decays, i.e. 
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we are interested in the exclusive decay products of a B meson containing at least one 
photon. The quark level decay is a flavour-changing-neutral-current process b — > 57 or 
b — > d'y. Such decays are rare, i.e. they come with small exclusive branching ratios of 
0(lCr 5 ), because they arise only at the loop level in the Standard Model. Thereby they 
test the detailed structure of the theory at the level of radiative corrections and provide 
information on the masses and couplings of the virtual SM or beyond-the-SM particles 
participating. Among the rare B decays the b —>■ sj modes are the most prominent ones 
because they are already experimentally measured. 

Primarily we are interested in the underlying decay of the heavy quark, which is 
governed by the weak interaction. But it is the strong force that is responsible for the 
formation of the hadrons that are observed in the detectors. If we want to do our exper- 
imental colleagues a favour, we let them measure the easier accessible exclusive decays, 
i.e. those where all decay products are detected. Then, however, we impose the burden 
of a more difficult theoretical treatment on ourselves. The "easier" option for theorists is 
to consider the inclusive decay, where e.g. for B — > all final states with a photon and 
strangeness content —1 have to be summed over. But quark-hadron duality allows us to 
consider instead of all the B decays only the parton-level b — > 57 decay, which amounts to 
a fully perturbative calculation. Much effort was put into the inclusive mode to achieve a 
full calculation at next-to-leading order in renormalizat ion-group- improved perturbation 
theory. Yet, for the exclusive decays we have to dress the b quark with the "brown muck", 
the light quark and gluon degrees of freedom inside the B meson, and have to keep strug- 
gling with hadronization effects. Despite the more complicated theoretical situation of the 
exclusive channels it is worthwhile to better understand them. Especially in the difficult 
environment of hadron machines, like the Fermilab Tevatron or the LHC at CERN, they 
are easier to investigate experimentally. In any case the systematic uncertainties, both 
experimental and theoretical, are very different for inclusive and exclusive modes. A care- 
ful study of the exclusive modes can therefore yield valuable complementary information 
in testing the Standard Model. 

The field-theoretical tool kit at our disposal for this analysis includes operator product 
expansion and renormalization group equations in the framework of an effective theory. 
Herewith the problem of calculating transition amplitudes can be separated into the per- 
turbatively calculable short distance Wilson coefficients and the long distance operator 
matrix elements. In principle, the latter ones have to be calculated by means of a non- 
perturbative method like lattice QCD or QCD sum rules. For the exclusive decays of B 
mesons, however, one can use additionally the fact that the 6-quark mass rnj, is large com- 
pared to Aqcd, the typical scale of QCD. This in turn allows one to establish factorization 
formulas for the evaluation of the relevant hadronic matrix elements of local operators 
in the weak Hamiltonian. Herewith a further separation of long-distance contributions 
to the process from a perturbatively calculable short- distance part, that depends only 
on the large scale m&, is achieved. The long-distance contributions have to be computed 
non-perturbatively or determined from experiment. However, they are much simpler in 
structure than the original matrix element. This QCD factorization technique was devel- 
oped by Beneke, Buchalla, Neubert, and Sachrajda for the non-leptonic two-body decays 
of B mesons [T5[] . We apply similar arguments based on the heavy-quark limit m& ^> Aqcd 
to the decays B — > K*j, B — > p7, and B — > 77. This allows us to separate perturbatively 
calculable contributions from the nonperturbative form factors and universal meson light- 
cone distribution amplitudes in a systematic way. With power counting in Aqcd/^6 we 
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can identify leading and subleading contributions. 

We have organized the subsequent 105 pages as follows: In the first part we fill our 
toolbox with the necessary ingredients. After a mini review of the Standard Model we 
present the basic equipment: operator product expansion, effective theories, and renor- 
malization group improved perturbation theory. A discussion of the effective b — > 57 
Hamiltonian and a short survey of its theoretical status quo concludes this chapter. Our 
best tool so far to treat the tough nut of exclusive B decays is QCD factorization. We de- 
vote Chapter [3] to the description of this useful technique. In order to be able to appreciate 
its merits we first present its predecessors "naive factorization" and generalizations. We 
then discuss rather in detail the QCD factorization approach, give a sample application 
to the calculation of the pion form factor, and mention limitations of QCD factorization. 
Finally, we comment on other attempts to treat hadronic matrix elements in exclusive 
non-leptonic B decays. 

Part II and III contain the main subject of this work: the treatment of B — > V'-f 
(V = K* or p) and B — ► 77 in QCD factorization. For both types of decays we first 
present the necessary formulas and basic expressions and then give numerical results and 
phenomenological applications. 

We derive for B — > V7 the decay amplitude complete at next-to-leading order in 
QCD and leading power in Aqcd/^6- Hard- vertex and hard-spectator corrections are 
discussed separately. Annihilation contributions turn out to be power suppressed, but 
nevertheless calculable. As they are factorizable and numerically important we include 
them for our phenomenological analysis. Doing so we become sensitive to the charge of the 
light spectator quark inside the B meson and can estimate isospin breaking effects. These 
turn out to be large. The most important phenomenological quantity we predict is the 
branching ratio. The NLL value is considerably larger than both the leading logarithmic 
prediction when the same form factor is used and the experimental value. Our calculation 
also allows us to estimate CP-asymmetries and U-spin breaking effects. We want to stress 
that this thesis is the first totally complete next-to-leading-logarithmic presentation of 
exclusive B — > V r y decays, and it is a model-indpendent one. 

The double radiative B — > 77 decays are analyzed with leading logarithmic accuracy 
in QCD factorization based on the heavy-quark limit m& ^> Aqcd- The dominant effect 
arises from the one-particle-irreducible magnetic-moment type transition b — > s(d)^ where 
an additional photon is emitted from the light quark. The contributions from one-particle 
irreducible diagrams are power suppressed but still calculable within QCD factorization. 
They are used to compute the CP-asymmetry in B — > 77 and to estimate so-called long- 
distance contributions in B and D — > 77. Numerical results are presented for branching 
ratios and CP asymmetries. 

We give our conclusions and an outlook in chapter [8[ 

Some more technical details are discussed in the Appendices. We give the explicit 
formulas for the Wilson coefficients, transformation relations between the two operator 
bases for b — > sj, and a one-loop proof that weak annihilation contributions to B — > V^j 
are calculable within QCD factorization. 



Part I 
Fundamentals 



Chapter 2 

The Basic Concepts 



In this chapter we briefly introduce the basic concepts needed for doing calculations in 
elementary particle physics. We present the Standard Model in a nutshell and introduce 
the concepts of renormalization, renormalization group, operator product expansion, and 
effective theories. We assume that the reader is familiar with quantum field and gauge 
theories and refer to the pertinent textbooks 
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2.1 The Standard Model 

As mentioned already in the introduction, the Standard Model (SM) is a comprehensive 
theory of particle interactions. Its success in giving a complete and correct description of 
all non-gravitational physics tested so far is unprecedented. In the following we give a short 
introduction into this beautiful theory. We want to introduce the "Old Standard Model," 
i.e. the one where neutrinos are massless. The recent evidence for neutrino masses, 
coming from the observation of neutrino oscillations |T7[ , has no direct consequences for 
our work. 

The Standard Model is made up of the Glashow-Salam- Weinberg Model |T| of elec- 
troweak interaction and Quantum Chromodynamics (QCD) ||. It is based on the princi- 
ple of gauge symmetry. The Lagrangian of a gauge theory is invariant under local "gauge" 
transformations of a symmetry group. Such a symmetry can be used to generate dynamics 
- the gauge interactions. The prototype gauge theory is quantum electrodynamics (QED) 
with its Abelian U(l) local symmetry. It is believed that all fundamental interactions are 
described by some form of gauge theory. 

The gauge group of strong interactions is the non-Abelian group SU(3)c which has 
eight generators. These correspond to the gluons that communicate the strong force be- 
tween objects carrying colour charge - therefore the "C" as subscript. Since the gluons 
themselves are coloured, they can directly interact with each other, which leads to the 
phenomena of "asymptotic freedom" and "confinement." At short distances, the coupling 
constant a s becomes small. This allows us to compute colour interactions using per- 
turbative techniques and turns QCD into a quantitative calculational scheme. For long 
distances on the other hand, the coupling gets large, which causes the quarks to be "con- 
fined" into colourless hadrons. In the words of Yuri Dokshitzer: "QCD, the marvellous 
theory of the strong interactions, has a split personality. It embodies 'hard' and 'soft' 
physics, both being hard subjects, the softer ones being the hardest." |T 
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Electroweak interaction is based on the gauge group SU(2)l <8> U(l)y - "L" stands 
for "left" and "Y" denotes the hypercharge - which is spontaneously broken to U(V)qed- 
This is achieved through the non vanishing vacuum expectation value of a scalar isospin 
doublet Higgs field [[HJ 

(2.1) 



Three of the four scalar degrees of freedom of the Higgs field give masses to the W and Z 
bosons. The remaining manifests itself in a massive neutral spin zero boson, the physical 
Higgs boson. It is the only particle of the Standard Model which lacks direct experimental 
detection. The current lower limit on its mass is 114.1 GeV at the 95% confidence level 
20fl . From electroweak precision data there is much evidence for a light Higgs. But as soon 



as such a light Higgs is found, this gives birth to the hierarchy problem. A scalar (Higgs) 
mass is not protected by gauge or chiral symmetries so we expect m R pa A pa 10 16 GeV if 
we do not want to fine-tune the bare Higgs mass against the mass aquired from quantum 
effects. Why should m# be much smaller than A? 

Fermions are the building blocks of matter. In the SM they appear in three generations 
which differ only in their masses. The two species of fundamental fermions are leptons and 
quarks. With regard to the gauge group SU{2) L the quarks and leptons can be classed 
in left-handed doublets and right-handed singlets. 



Quarks: 



u \ I c \ ft 



Leptons: 



d' I V s' I V V „ 

UR Cr t R 

d R s R b R 



V e \ ( Vn \ ( Vt 



L 



V Jr. \ T 



dR UR TR 



Within the Old Standard Model there are no right-handed neutrinos. The quarks carry 
colour charge and transform as SU(3)c triplets whereas the colourless leptons are SU(3)c 
singlets. 

Fermion masses are generated via a Yukawa interaction i/j(x)4>(x)il)(x) with the Higgs 
field fl2.ip . Using global unitary transformations in flavour space, the Yukawa interaction 
can be diagonalized to obtain the physical mass eigenstates 

d'\ ( V ud V us V ub \ ( d 

s' = V cd V cs V cb • Is | (2.2) 

v V v td v ts v tb \b 



= v< 



CKM 



The non-diagonal elements of the Cabibbo-Kobayashi-Maskawa matrix |H] allow for tran- 
sitions between the quark generations in the charged quark current 

d 

J t = («,c,t) L 7 M VcKM [ s I (2.3) 



2.1 The Standard Model 
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Figure 2.1: Example of a box and penguin diagram in the full theory. CP violating effects 
are possible because quarks of all three generations are involved in the processes. 

Unitarity of Vqkm guaranties the absence of flavour-changing-neutral-current (FCNC) 
processes at tree level. This Glashow-Iliopoulos-Maiani (GIM) mechanism [^2| would 
forbid FCNC transitions even beyond the tree level if we had exact horizontal flavour 
symmetry which assures the equality of quark masses of a given charge. Such a symmetry 
is in nature obviously broken by the different quark masses so that at the one-loop level 
effective b — > s, b — > d, or s — > d processes like B — ► X s 7 or K — > i\vv can appear. 

A unitary complex N x N matrix can be described by N 2 real parameters. If this 
matrix mixes N families each with two quarks one can remove 2N — 1 phases through 
a redefinition of the quark states, which leaves the Lagrangian invariant. Because an 
orthogonal N x N matrix has iV(iV — l)/2 real parameters (angles) we are left with 
N 2 - (2N - 1) - iV(iV - l)/2 = (N — 1)(N - 2)/2 independent physical phases in the 
quark mixing matrix. Therefore, it is real if it mixes two generations only. The three- 
generation CKM matrix, however, has to be described by three angles and one complex 
phase. The latter one is the only source of CP violation within the Standard Model if 
we desist from the possibility that #qcd 7^ 0. But these CP-violating effects can show 
up only if really all three generations of the Standard Model are involved in the process. 
Typically this is the case if one considers loop contributions of weak interaction, like box 



or penguin diagrams as in Fig. 2.1 



In principle there are many different ways of parametrizing the CKM matrix. For 



practical purposes most useful is the so called Wolfenstein parametrization [23 



1-f A A\ 3 (p-ir]) 

Vckm =| -A 1-f AX 2 | (2.4) 

AX*{\- p-in) -AX 2 1 

which is an expansion to O (A 3 ) in the small parameter A = \ V US \ ~ 0.22. It is possible to 
improve the Wolfenstein parametrization to include higher orders of A |[24|| . In essence p 
and rj are replaced with p = p(l — A 2 /2) and f] = r)(l — A 2 /2), respectively. 

For phenomenological studies of CP-violating effects, the so called standard unitarity 
triangle (UT) plays a special role. It is a graphical representation of one of the six unitarity 
relations, namely 

V ud V,: b + V cd V* b + V td V t t = (2.5) 

in the complex (p, rj) plane. This unitarity relation involves simultaneously the elements 
V u b, V c b, and Vtd which are under extensive discussion at present. The area of this and 
all other unitarity triangles equals half the absolute value of Jcp = Im(V us V cb V* b V* s ) , the 
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2. The Basic Concepts 




Figure 2.2: Unitarity Triangle 



Jarlskog measure of CP violation po[ . Usually, one chooses a phase convention where 
V cd V^ is real and rescales the above equation with |V^V^| = A\ 3 . This leads to the trian- 
gle in Figure |2.2| with a base of unit length and the apex (p,f]). A phase transformation 
in (|2.5|) only rotates the triangle, but leaves its form unchanged. Therefore, the angles 
and sides of the unitarity triangle are physical observables and can be measured. Much 
effort was and is put into the determination of the UT parameters. One tries to measure 
as many parameters as possible. The consistency of the various measurements tests the 
consequences of unitarity in the three generation Standard Model. Any discrepancy with 
the SM expectations would imply the presence of new channels or particles contributing 
to the decay under consideration. So far, all experimental results are consistent with the 
Standard Model picture p^j27|] . The state-of-the-art "frequentists 
ity triangle from the 2002 Winter conferences is displayed in Fig. 



result for the unitar- 



3fJ. Actually the 

good agreement of measurements with the Kobayashi-Maskawa mechanism gives rise to 
some theoretical puzzles: the KM mechanism for example does explain neither the cosmic 
baryon asymmetry nor the smallness of #qcd an d basically all extensions of the Standard 
Model introduce a large number of new CP-violating phases. 



2.2 Renormalization and Renormalization Group 

Given the Lagrangian of a theory one can deduce the Feynman rules by means of which 
amplitudes of the processes occuring in this theory can be calculated in perturbation the- 
ory. In Feynman diagrams with internal loops, however, one often encounters ultraviolet 
divergences. This is because the momentum variable of the virtual particle in the loop 
integration ranges from zero to infinity. The theory of renormalization is a prescription 
which allows us to consistently isolate and remove all these infinities from the physically 
measurable quantities. A two-step procedure is necessary. 

First, one regulates the theory. That is, one modifies it in a way that observable 
quantities are finite and well defined to all orders in perturbation theory. We are then free 
to manipulate formally these quantities, which are divergent only when the regularization 
is removed. The most straightforward way to make the integrals finite is to introduce a 
momentum cutoff. But this violates for example Lorentz invariance or the Ward identities. 



2.2 Renormalization and Renormalization Group 
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IP" 




Figure 2.3: The Range-fit result of the global CKM fit from CKMfitter (status May 2002) 
26|. 



A regularization method that preserves all symmetries of a gauge theory is dimensional 
regularization |2"E| , |2Tj| | . The basic idea is to compute the Feynman diagram as an analytic 
function of the dimensionality of space time D = 4 — 2e. For sufficiently small D, any 
loop-momentum integral will converge. The singularities are extracted as poles for e — > 0. 

Potential problems of dimensional regularization concern the treatment of 75 in D ^ 4 
dimensions. The definition 

(2.6) 



75 



4! 



with £ K \fiv the completely antisymmetric tensor in four dimensions, cannot straightfor- 
wardly be translated to D 7^ 4 dimensions. In the so called "naive dimensional regular- 
ization" (NDR) scheme |50| the metric tensor is generalized to D dimensions: g£ = D, and 
the 7 matrices obey the same anticommuting rules as in four dimensions. Even if these 
rules are algebraically inconsistent [31|, the NDR scheme gives correct results provided 
one can avoid the calculation of traces like tr (757^7^7^7^) [[321 . 

The scheme originally proposed by 't Hooft and Veltman (HV scheme) p8l allows a 



consistent formulation of dimensional regularization even when 75 couplings are present 
3j|. Besides the D- and 4-dimensional metric tensors g and g one introduces the — 2e- 



dimensional tensor g. One can then split the D-dimensional Dirac matrix 7 M into a 4- and 
a — 2e-dimensional part 7 M and 7 M which separately obey anticommutation relations with 
the appropriate metric tensors. A 75 can be introduced which anticommutes with 7 but 
commutes with 7. The price one has to pay for a consistent dimensional regularization 
scheme is a substantial increase in the complexity of calculations. 

The second step in our programme to eliminate the infinities from a theory is renormal- 
ization. This is the process of relating the unphysical (bare) and physical (renormalized) 
parameters like couplings g or masses m and rewrite observables as functions of the phys- 
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ical quantities. The renormalization procedure hides all divergences in a redefinition of 
the fields and parameters in the Lagrangian, i.e. 

g = Z g gfi e m = Z m m 

qo = zy\ Aq = zfA* (2 - 7) 

and thus guaranties that measurable quantities stay finite. The index "0" indicates bare 
quantities. Introducing the parameter /i with dimension of mass in (|2.7|) is necessary to 
keep the coupling dimensionless. The factors Z are the renormalization constants. The 
renormalization process is performed recursively in powers of the coupling constant g. If 
at every order of perturbation theory all divergences are reabsorbed in Z's, the theory is 
called "renormalizable" . Theories with gauge symmetries, like the Standard Model, are 
renormalizable. This is true even if the gauge symmetry is spontaneously broken via the 
Higgs mechanism because gauge invariance of the Lagrangian is conserved ||33|| . 

Renormalization can be straightforwardly implemented via the counter-term method. 
According to ( |2.7| ) the unrenormalized quantities are reexpressed through the renormalized 
ones in the original Lagrangian. Thus 

£() = £ + ^counter (2-8) 

The counter terms £ C ounter are proportional to (Z—l) and can be treated as new interaction 
terms. For these new interactions Feynman rules can be derived and the renormalization 
constants Z{ are determined such that the contributions from these new interactions 
cancel the divergences in the Green functions. This fixes the renormalization constants 
only up to an arbitrary subtraction of finite parts. Different finite parts define different 
renormalization schemes. In the Minimal Subtraction (MS) scheme only the divergences 
and no finite parts are subtracted 3J|. The modified MS scheme (MS) JEJ defines the 



finite parts such that terms ln47r — 7^, the artifacts of dimensional regularization, vanish. 
This can be achieved if one calculates with 

fi e lE/2 

V47T 



instead of \i and performs minimal subtraction afterwards. We will exclusively work with 
the MS scheme in the following. 

Every renormalization procedure necessitates to introduce a dimensionful parameter 
\x into the theory. Even after renormalization the theoretical predictions depend on this 
renormalization scale \i. At this momentum scale the renormalization prescriptions, which 
the parameters of a renormalized field theory depend on, are applied. One "defines 
the theory at the scale //." The bare parameters are //-independent. To determine the 
renormalized parameters from experiment, a specific choice of /1 is necessary: g = g(fi), 
m = m(/i), q = q(fi). Different values of \i define different parameter sets g(/i), m(fj,), 
q{n). The set of all tranformations that relates parameter sets with different /i is called 
renormalization group (RG). 

The scale dependence of the renormalized parameters can be obtained from the /i- 
independence of the bare ones. In QCD we get from ( [2.7|) the renormalization group 
equations (RGE) for the running coupling and the running mass 



dlnfi 
GHn/i 



(3(g(fi),s) (2.10) 
-lm{g{^))m{p) (2.11) 
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with the /3-function 



P{g{fj) i e) = -eg- ' (,Z " 



Z g d In /i 

and the anomalous dimension of the mass operator 



(2.12) 



im{g(fj>)) 

Calculating to two-loop accuracy we get 

0(g) 

7m(a s ) 



167T 



1 dZ m 
Z m d In fi 



;0o 



(2.13) 



16tt 



2^ 2 



01 



%(o) + f^V 7 (i) 



where 



A) 



11JV-2/ 



7<? = 6C F 



Pi 

7 



|w 2 - Hjv/ - 2C F f 



4tt 



N 2 



2N 



(2.14) 
(2.15) 

(2.16) 
(2.17) 
(2.18) 



with N the number of colours and / the number of active flavours. The solutions for 

ft lnln^/A 2 ^) 



a s (//) and m(/x) then are |35j 

4:71 



a 8 (ji) 



/3 ln(/i 2 /Ay 



1 - 



/?o 2 ln(/i 2 /A 



— ) 

MS 1 



m(n) = m(/i ) 



7^ - 

2/3 



1 + 



7m A 7m \ a s (/x) - a s (fi ) 



20 o 201 



4tt 



(2.19) 
(2.20) 



Here is a characteristic scale both for QCD and the used MS scheme and depends 
also on the number of effective flavours present in O and 0\. An cta(Mz) = 0.118±0.005 
corresponds to A^L = 225^™ MeV in NLO. It is interesting to note that such a mass 
scale A emerges without making reference to any dimensional quantity and would be 
present also in a theory with completely massless particles. In QCD with three colours, 
even for six active flavours, both O = 7 and /20 o = 4/7 are positive. This leads 
to asymptotic freedom as the coupling tends to zero with increasing fi. The pole at 
^■Ws signals the breakdown of perturbation theory but gives a plausible argument for 
confinement. Similarly, m(fi) decreases with \x getting larger. 

A particularly useful application of the renormalization group is the summation of 
large logarithms. To see this we reexpress ot s of Q2.19| ) as 



1 - 



0o 4vr v{h) . 



[2.21] 
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with 

^ = 1-/30^ In 4 (2-22) 
If we expand the leading order term of ( p.21| ) in a s (fi ) we get 

a s (fi) = a s (fi )^U^ln^) (2.23) 

m=0 ^ ^ " ' 

Thus the solution of the RGE automatically sums the logarithms ln(//o/ fi 2 ) which get 
large for fi <C fi$. Generally, solving the RGE to order n sums in a s (fi) all terms of the 
form 

(2 \ k 
a s (^ )ln^|J , 0<m<n, fcGff„ (2.24) 

This is particularly useful if, though a s (/io) is smaller than one, the combination 
a s (fj,o)\ia(fil/ n 2 ) is close to or even larger than one. Then the large logarithms would 
spoil the convergence of the perturbation series. 



2.3 Operator Product Expansion 

Up to now we treated processes of strong and electroweak interaction separately. But 
all weak processes involving hadrons receive QCD corrections, which can be substantial 
especially for non-leptonic and rare decays. The underlying quark level decay of a hadron 
is governed by the electroweak scale given by M W)Z = 0(100 GeV). On the other hand, 
the available energy inherent in a B meson decay is of O(mg). In dimensional regular- 
ization for example we encounter logarithms of the ratio of either of these scales with the 
renormalization scale fi. If the scales involved are widely separated it is not possible to 
make all the logarithms small by a suitable choice of the renormalization scale. As we 
have seen in the last section these large logarithms can be summed systematically using 
renormalization group techniques. 

But we have yet another energy scale in the problem. A priori we cannot consider 
the decay of free quarks. Due to confinement quarks appear in colourless bound systems 
only. The binding of the quarks inside the hadron via strong interaction is characterized 
by a typical hadronic scale of (9(1 GeV). Here, even without large logarithms the strong 
coupling a s is too large for perturbation theory to make sense. Unfortunately, in many 
cases the non-perturbative methods we have at hand nowadays are not yet developed 
enough to give accurate results. 

Coming back to the typical energy msinaB decay. Do we have to know at all what is 
really going on at energies of (9(100 GeV) or the corresponding extremely short distances? 
In fact we do not. We also do not bother general relativity to calculate the trajectory of 
an apple falling from a tree or QED and QCD to learn something about the properties 
of condensed matter. Instead, we employ Newtonian mechanics or the laws of chemistry 
and solid state physics, respectively. They are nonrelativistic approximations or effective 
theories appropriate for the low energy scale under consideration. This is exactly what 
we want to achieve for the weak interaction of quarks as well. The theoretical tool for 
this purpose is operator product expansion (OPE) |36| which we shall introduce in the 
following. 
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Figure 2.4: Replacing a W propagator with an effective four-fermion vertex. 

For small separations, the product of two field operators A(x) and B(y) can be ex- 
panded in local operators Qi with potentially singular coefficient functions q as 

A(x)B(y) = Ci(x - y)Qi(x) (2.25) 

i 

The Wilson coefficient c\(z) accompanying the operator with lowest energy dimension 
is the most singular one for z — > and the degree of divergence of the q decreases 
for increasing operator dimension. Furthermore, for dimensional reasons, contributions of 
operators with higher dimension are suppressed by inverse powers of the heavy mass (small 
distance) scale. In principle, we have to consider all operators compatible with the global 
symmetries of the operator product AB. The physical picture is that a product of local 
operators should appear as one local operator if their distance is small compared to the 
characteristic length of the system. One can systematically approximate the behaviour of 
an operator product at short distances with a finite set of local operators. This is exactly 
what is done for the theory of weak decays. 

Here, the mass of the W boson M\y ~ 80 GeV is very large compared to a typical 
hadronic scale. Therefore, the W propagator A Mly (x, y) is of very short range only. In 
the amplitude of weak decays it connects two charged currents J^(x) and J„(y), which 
hence interact almost locally so that we can perform an OPE. Let us consider the quark 
level decay b — > csu for definiteness. The tree-level ^-exchange amplitude for this decay 
is given by 



G M 2 
A{b^csu) = --J=V cb V* S k2 _ w M ^ {su) V - A (cb) v _ A 



l^ i^V^ +O^ ) (2-26) 
local operator 



Since the momentum transfer k through the W propagater is small as compared to Myy, 
we can safely neglect the terms 0{k 2 / M"^). The W propagator then quasi shrinks to a 
point (see fig. |2.4p and we obtain an effective four-fermion interaction. This is the modern 



formulation of the classical Fermi theory of weak interaction with Gp = 1.166- GeV -2 
the Fermi constant. The Wilson coefficient in this example is simply one. The notation 
(QiQ2)v-A in ( 2.26 ) is a practical shorthand for a left-handed charged quark current with 



the chiral vector minus axialvector structure 



(<?i<?2)v-a := 9i7m(1 ~ 75)92 



(2.27) 
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2.4 Effective Theories 

The result (|2.26|) can also be derived from an effective Hamiltonian 
Gp 

7i c s = —/=V c bV* s (su)v-A(cb)y_A + operators of higher dimension (2.28) 
V2 

where the operators of higher dimensions correspond to the terms Oik 2 jM^) in ( |2.26| ) 
and can likewise be neglected. In the effective theory the W boson is removed as an 
explicit, dynamical degree of freedom. It is "integrated out" or "contracted out" using 
the language of the path integral or canonical operator formalism, respectively. One can 
proceed in a completely analougous way with the heavy quarks. This leads to effective 
f quark theories where / denotes the "active" quarks, i.e. those that have not been 
integrated out. 

If we include also short distance QCD or electroweak corrections more operators have 
to be added to the effective Hamiltonian which we generalize to 



^E y cK M a(/i)Q^) (2.29) 



Here the factor Vq KM denotes the CKM structure of the particular operator. If we want 
to calculate the amplitude for the decay of a meson M = K, D, B, . . . into a final state 
F we just have to project the Hamilton operator onto the external states 



A(M — > F) = (F\H cS \M) 
G p 

= 71 



% E V&cMCifci) (F\Qi{»)\M) . (2.30) 



The Wilson coefficients Cj(/x) can be interpreted as the coupling constants for the effective 
interaction terms Qi(fi). They are calculable functions of a s , Myy, and the renormalization 
scale jji,. To any order in perturbation theory the Wilson coefficients can be obtained by 
matching the full theory onto the effective one. This simply is the requirement that the 
amplitude in the effective theory should reproduce the corresponding amplitude in the 
full theory. Hence, we first have to calculate the amplitude in the full theory and then the 
matrix elements (Qi). In this second step the resulting expressions may, even after quark 
field renormalization, be still divergent. Consequently we have to perform an operator 
renormalization 

Q? ] = Z i3 Q 3 (2.31) 

where Qf^ denotes the unrenormalized operator. This notation is somewhat sloppy and 

misleading. What actually is renormalized is not the operator but the operator matrix 

elements, or, even more exactly, the amputated Green functions (Qi). Then we have to 

i li 

include the renormalization constant Z q for each of the four external quark fields: 

(Qi)V> = Z- 2 Z l3 (Q 3 ) (2.32) 

In general, the renormalization constant is a matrix so that operators carrying the 
same quantum numbers can mix under renormalization. Operators of a given dimension 
mix only into operators of the same or of lower dimension. Again, the divergent parts of 
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the renormalization constant are determined from the requirement that the amplitude in 
the effective theory is finite. The finite part in Z^ on the other hand defines a specific 
renormalization scheme. In a third step we extract the Wilson coefficients by comparing 
the full and the effective theory amplitude. These are the Wilson coefficients at some fixed 
scale fiQ. A caveat here is that the external states in the full and the effective theory have 
to be treated in the same manner. Especially the same regularization and renormalization 
schemes have to be used on both sides. 

As the Wilson coefficients appear already at the level of the effective Hamiltonian, 
they are independent of the external states this Hamiltonian is projected onto to ob- 
tain the complete amplitude. When determining the Wilson coefficients, any external, 
even unphysical, state can be used. The coefficient functions represent the short- distance 
structure of the theory. Because they depend for example on the masses of the par- 
ticles that were integrated out, they contain all information about the physics at the 
high energy scale. The long-distance contribution, on the other hand, is parametrized 
by the process-dependent matrix elements of the local operators. This factorization of 
SD and LD dynamics is one of the salient features of OPE. We can calculate the Wilson 
coefficients in perturbation theory and the hadronic matrix elements by means of some 
non-perturbative technique like 1/N expansion, sum rules, or lattice gauge theory. Es- 
pecially to use the latter one, a separation of the SD part is essential for today's lattice 
sizes. The factorization can be visualized with large logarithms h^M^/m 2 ) being split 
into ln(M^//i 2 ) + ln(/z 2 /m 2 ). In doing so, the first logarithm will be retrieved in the 
Wilson coefficients and the second one in the matrix elements. From this point of view 
the renormalization scale 11 can be interpreted as the factorization scale at which the full 
contribution is separated into a low energy and a high energy part. 

A typical scale at which to calculate the hadronic matrix elements of local operators 
is low compared to M\y- For B decays we would choose /i = 0(i71b). Therefore, the 
logarithm h^M^/yU 2 ) contained in the Wilson coefficient is large. So why not use the 
powerful technique of summing large logarithms developped in section |2.2| ? In order to 
do so we have to find the renormalization group equations for the Wilson coefficients 
and solve them. But so far the Wilson coefficients were not renormalized at all. If 
we remember, however, that in the effective Hamiltonian the operators, which have to be 
renormalized, are accompanied always by the appropriate Wilson coefficent we can shuffle 
the renormalization as well to the Wilson coefficients. Let us start with the Hamiltonian 
of the effective theory with fields and coupling constants as bare quantities, which are 
renormalized according to 



i 



(0) 




(2.33) 
(2.34) 



Then the Hamiltonian Q2.29| ) is in essence 



H eS oc Cl 0) Qi(q®) 



CiQi + (Z q Z?j — Sij)CjQi 



(2.35) 



i.e. it can be written in terms of the renormalized couplings Cj and fields Qi plus countert- 
erms. The g^ -* indicates that the interaction term Qi is composed of unrenormalized fields. 
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If we calculate the amplitude with the Hamiltonian ( 2.35 ) including the counterterms, we 
get the finite renormalized result 



[2.36) 



Comparing with ( |2.32| ) we read off the renormalization constant for the Wilson coefficients 

Zti = Zji (2-37) 

So we can think of the operator renormalization in terms of the completely equivalent 
renormalization of the coupling constants Cj, as in any field theory. If we again demand 
the unrenormalized Wilson coefficients not to depend on [i we obtain the renormalization 
group equation 

,/r ' (/ " - 7i<(M)C,-(A*) , (2.38) 



c?ln/i 

with the anomalous dimension matrix for the operators 



i dZ, 



ik 



(iln/x 



(2.39) 



Let us simply state here that the numerical values for the 7^ can be determined directly 
from the divergent parts of the renormalization constants Zij. In ( |2.38 ) the transposed 
of this anomalous dimension matrix appears. It is only the sign and the fact that the 
anomalous dimension is a matrix instead of a single number that distinguishes the RGE 
for the Wilson coefficients from that of the running mass in ( 2.1 1| ) . Therefore, we could 
use the solution ( 2.20Q with the appropriate changes. To leading order this is in fact 
possible. But if we want to go to next-to-leading-order accuracy we run into problems, 
because the matrices 7-^ and 7^- in the perturbative expansion 



_ (0) «£_ (1) / fts 



+ 0{ct s 



(2.40) 



do not commute with each other. Let us instead formally write the solution for the Wilson 
coefficients with an evolution matrix U ((/,, Hq) 



Ci(fi) = Uijifi, /io)Cj(Ak)) 
The leading order evolution matrix can be read off from (|2.20|) 



(2.41) 



^ (0 Wo) 




(2.42) 



where V is the matrix that diagonalizes ^^> T 



7 g> = v-'j^v 



(2.43) 
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and 7^ is the vector containing the eigenvalues of 
solution we make the clever ansatz 



For the next-to-leading order 



i + j 

47T 



t/ ( °Wo) 



47T 



J 



which proves to solve (2.38[) if P?| 



where the elements of H are 



Hi 



a 



J = VHV- 1 



(0) A 



m 2/3o + 7] 



(o) 



7 



(0) 



(2.44) 



(2.45) 



;2.46) 



with 



G = V-^WV (2.47) 
As we have mentioned in section |2.2|, the procedure of renormalization allows to subtract 



arbitrary finite parts along with the ultraviolet singularities. Whereas physical quantities 
must clearly be independent of the renormalization scheme chosen, at NLO unphysical 
quantities, like the Wilson coefficients and the anomalous dimensions, depend on the 
choice of the renormalization scheme. To ensure a proper cancellation of this scheme 
dependence in the product of Wilson coefficients and matrix elements the same scheme 
has to be used for both. In order to uniquely define a renormalization scheme it is not 
sufficient to quote only the regularization and renormalization procedure but one also has 
to choose a specific form for the so-called evanescent operators. These are operators which 
exist in D 7^ 4 dimensions but vanish in D = 4 P21RS,R5|. 



So what do we have achieved so far? We have determined the Wilson coefficients at 
a scale /io via a matching procedure. These are the initial conditions for the evolution 
from /i down to an appropriate low energy scale /x via U(fj,,fjbo) which sums large loga- 
rithms. Herefore, we had to determine the anomalous dimensions of the operators and 
solve the renormalization group equation for the Wilson coefficients. We thus arrive at a 
RG improved perturbation theory and officially don't speak any more of "leading" (LO) 
and "next-to-leading order" (NLO) but rather of "leading" (LL) and "next-to-leading- 
logarithmic order" (NLL). Yet, we might carelessly use the terms synonymously. In our 
task to evaluate weak decay amplitudes involving hadrons in the framwork of a low en- 
ergy effective theory we then only lack the calculation of the hadronic matrix elements 
(Qi (/■*)}■ This, however, is a highly non-trivial problem which this and many other works 
are devoted to. 

Looking at the complicated NLO formulas one might ask why at all going to next-to- 
leading order accuracy. After all these calculations imply the evaluation of two or even 
more loop diagrams which are technically very challenging. But they are very important. 
First of all we can test the validity of the renormalization group improved perturbation 
theory. Then, of course, we hope that the theoretical uncertainties get reduced. One 
particular issue is the residual renormalization scale dependence of the result. The scale 
H enters for example in ck s (/j) or the running quark masses, in particular m t (/j), m&(/x), 
and m c (/i). In principle, a physical quantity cannot depend on the renormalization scale. 
But as we have to truncate the perturbative series at some fixed order, this property is 
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broken. The renormalization scale dependence of Wilson coefficients and operator matrix 
elements cancels only to the order of perturbation theory included in the calculation. 
Therefore, one can use the remaining scale ambiguity as an estimate for the neglected 
higher order corrections. Usually one varies \i between half and twice the typical scale 
of the problem, i.e. m b /2 < p, < 2m b for B decays. Going to NLO significantly reduces 
these scale ambiguities. Furthermore, the renormalization scheme dependence of the 
Wilson coefficients appears at NLO for the first time. Only if we properly match the long 
distance matrix elements, obtained for example from lattice calculations, to the short 
distance contributions, these unphysical scheme dependences will cancel. Another issue 
is that the QCD scale A^g, which can be extracted from various high energy processes, 
cannot be used meaningfully in weak decays without going to NLO. 



2.5 The Effective b —> S7 Hamiltonian 



In this section we want to discuss the effective Hamiltonian necessary for the calculations 
to follow. For b — > S7 transitions it reads 



Ti-cs = —7= ^ 



s/2 



p=u,c 



C1Q1 + C2Q2 + CiQi 

i=3,...,8 



where 



a? } = v; s v ph 

The operators originate from the diagrams in Fig. |2.5| and are given by 



2.48) 

2.49) 

2.50) 
2.51) 
2.52) 

2.53) 

2.54) 

2.55) 

2.56) 
2.57) 



with e and g s the coupling constants of electromagnetic and strong interaction and 
and G^y the photonic and gluonic field strength tensors, respectively. In Q 7 and Q 8 
we neglected the very small m s (l — 75) contribution. The i,j are colour indices. If no 
colour index is given the two operators are assumed to be in a colour singlet state. The 
operator basis fl2.5C - 2.57 ) consists of all possible gauge invariant operators with energy 
dimension six with the following properties: they have the correct quantum numbers to 
contribute to b — > 57, they are compatible with the symmetries of electroweak interaction, 



Qi = 
Ql = 
Qs = 


(sp) v _ A (pb) v _ A 

(siPj) V -A(pjbi)V-A 

(sb)v- A ^(qq)v- A 
g 


Qa = 


(s i b j ) v _ A ^2(q j q i ) v _ A 


Q 5 = 


(sb) V -A^2(qq)v+A 
q 


Qe = 


(sibj) V -A /iQjQi)v+A 
q 


Qi = 


g 

SiCr MI/ (l + 75)6, F^v 


Q 8 = 
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(a) (b) (c) (d) 




(e) (f) (g) (h) 



Figure 2.5: The diagrams where the operator basis for b — > S7 originates from: Current- 
current diagram (a) with QCD corrections (b), (c), (d); gluon penguin diagram (e), mag- 
netic photon penguin diagrams (f), (g), and magnetic gluon penguin diagram (h). The 
cross in (f), (g), and (h) denotes that the mass of the external b quark has to be kept. 
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and they cannot be transformed into each other by applying equations of motion. As a 
consequence, our operator basis is only the correct one if all external states are taken 
on-shell [fiO| |. Additional operators have to be considered if an off-shell calculation is 



performed. There is yet another operator basis used in the literature. It was introduced 



by Chetyrkin, Misiak, and Miinz (CMM) |f41|j42| because then no Dirac traces containing 



75 arise in effective theory calculations, wich allows to use fully anticommuting 75 in 
dimensional regularization. It reads 

(2.58) 
(2.59) 
(2.60) 

(2.61) 

(2.62) 

(2.63) 

(2.64) 
(2.65) 

where T a stand for SU(3) co \ OT generators and L = (1 — 7s)/2 and R = (1 + 7s)/2 for 
the left and right-handed projection operators. We denote the corresponding Wilson 
coefficients with Zj. The operator basis (|2.58H2T65|) is in principle the more natural one 



p p = 

A 1 


(sr'y„T a vr)(vT'y fl T a bT) 


pp _ 


(sL7 M PL)fe7 M ^) 


P 3 = 


(sl7A) X^^) 

q 


Pa = 


(s L l,T a b L ) ^(g 7 M T a g) 
q 


P 5 = 


q 


P 6 = 


(sLl, 1 l, 2 l,J a b L )J2( ( ll >11 l tl2 l fl3 T a q) 
q 


P 7 = 


167T Z 


Ps = 


^m b {s L a^T% R )Gl v 



as the operators appear exactly in this form when calculating the diagrams of Fig. |2.5f In 
our normal operator basis ( |2.50| - |2~57|) , the following four-dimensional identity was used 
to "simplify" P 5 and P 6 

liflvlp = g^ulp + gvpln - g^plu + i^^ P l a l5 in D = A (2.66) 

This step, however, requires to introduce several more evanescent operators and leads to 
problematic traces with 75 in two-loop calculations. The choice of the operator's colour 
structure is more natural in the CMM basis, too, as it is the one emerging from the 
diagrams in Fig. J!5". 

In the common nomenclature Q\ and Q2 are called current- current operators, Q3...6 
QCD penguin operators and Qj and Q$ electromagnetic and chromomagnetic penguin 
operator, respectively. For historical reasons the numbering of Q\ 2 is sometimes reversed 
in the literature ||3|| . The sign conventions for the electromagnetic and strong couplings 



correspond to the covariant derivative = + ieQfA^ + igT a A^. The coefficients C 7 
and C§ then are negative in the Standard Model, which is the choice generally adopted 
in the literature. The effective Hamiltonian for b — > d'-f is obtained from equations ( |2.48| - 
2.65|) by the replacement s — > d. 



Let us summarize the status quo in determining the effective Hamiltonian for b 57. 
At leading order only Q-j contributes to b — > 57. The corresponding Wilson coefficient was 
first calculated by Inami and Lim in 1980 ||44fl. For the leading logarithmic RG improved 
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calculation the 8x8 anomalous dimension matrix was obtained bit by bit. Mixing of the 
four-quark operators Qi..s among each other was considered already earlier for the analysis 
of nonleptonic decays [[45] . The 2x2 submatrix for the magnetic penguins was obtained by 
Grinstein, Springer, and Wise JRJ. Because the mixing of Qi..s into Qj and Q 8 vanishes 
at the one- loop level, one has to perform two- loop calculations if the LL result for CV(jU&) 
is wanted. This technical complication delayed the first completely correct result for the 
8x8 anomalous dimension matrix until 1993 E7j followed by independent confirmations 



Going to next-to-leading order accuracy was highly desirable, because the leading 
logarithmic expression for the branching ratio B(B — > X s 7) suffers from sizable renor- 
malization scale uncertainties at the ±25% level. The step from leading order to leading 
logarithmic order already increased the branching ratio B(B — > X s ^y) by almost a fac- 
tor of three p9|-|5T|. The high-flying NLL enterprise was a joint effort of many groups. 
QCD corrections affect both Wilson coefficients and operator matrix elements. The 0{ot s ) 
matching for the current- current operators was calculated already a long time ago [|52|| , the 
one for the QCD penguins more than ten years later [53fl. Two-loop matching is necessary 



for Qi and Qs, which was achieved first by Adel and Yao |3J and subsequently checked 
thoroughly |55|j56|| . The 0(a 2 s ) contributions to the 6x6 anomalous dimension matrix 
were calculated in p2| , |52| , p3| , p7| , the submatrix for Q 7 and Q% in [j5"8] . Chetyrkin, Misiak, 



and Munz finally succeeded in determining the three-loop mixing of four-quark operators 
into magnetic penguin operators [IT]. The explicit formulas for the Wilson coefficients 



needed subsequently can be found in Appendix and in Appendix |B| we comment on 
the transformation properties of the two operator bases used in this work. 

For the inclusive decay B — > X^, i.e. the radiative decay of a B meson into the sum 
of final states with strangeness S = — 1, the operator matrix elements can be computed 
perturbatively employing the heavy-quark expansion (HQE) fS"9]| . This method consists 
of an OPE in inverse powers of the large dynamical scale of energy release ~ m?, followed 
by a nonrelativistic expansion for the b field. Using the optical theorem, the inclusive 
.B-decay rate can be written in terms of the absorptive part of the forward scattering 
amplitude (B\T\B). The transition operator T is the absorptive part of the time-ordered 
product of H c s(x)H e G(0) 



T = Im 



d 4 xTH efi (x)H e{i (0) 



(2.67) 



If we insert a complete set of states inside the time-ordered product we see that this was 
just a fancy way of writing the standard expression for the decay rate 

F(B - X) = J_ ^(27r) 4 5 4 (p B -p x ) \(X\H cS \B)\ 2 (2.68) 
a x 

However, the formulation in terms of the T product allows for a direct evaluation using 
Feynman diagrams. Because of the large mass of the b quark, we can construct an OPE in 
which T is represented as a series of local operators containing the heavy-quark fields. The 
operator of lowest dimension is bb. Its matrix element is simplified by a nonrelativistic 
expansion in powers of 1/rrit, starting with unity. In the limit of an infinitely heavy b 
quark the B meson decay rate is therefore given by the b quark decay rate. This quark- 
hadron duality is of great use for many inclusive calculations and justifies for example 
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the spectator model. Corrections to this relation appear only at 0(l/m^) because the 
potential dimension 4 operator bJpb can be reduced to mjjb. However, the OPE only 
converges for sufficiently inclusive observables. But for B — > X s ^ decays experimental 
cuts are necessary to reduce the background from charm production. This restricts the 
available phase space considerably leading to complications of the OPE based analysis. 
The virtual corrections to the matrix elements (s~f\Q\j$\b) were first calculated by 



Greub, Hurth, and Wyler |60| and checked by Buras, Czarnecki, Misiak, and Urban using 
another method pl|. The latter authors presented recently also the last missing item 
in the NLO analysis of B — > X s 7, namely the two-loop matrix elements of the QCD- 
penguin operators Yet, their contribution is numerically very small because the 

corresponding Wilson coefficients almost vanish. The Bremsstrahlung corrections, i.e. 
the process b — > s^g, influence especially the photon energy spectrum. At leading order 
this spectrum is a 6 function, smeared out by the Fermi motion of the b quark inside the 
B meson, whereas it is broadened substantially at NLO |)3[ . 



Even higher electroweak corrections have been calculated for B — > X s 7 |63J]. Finally, 
there are non-perturbative contributions which can be singled out in the framework of 
HQE. The corrections mainly account for the fact that in reality a B meson and 

not a b quark is decaying. Additionally, one has to consider long distance contributions 
originating in the photon coupling to a virtual cc loop. This Voloshin effect is proportional 
to 1/ml and enhances the decay rate by 3% p3HB^j. 



c 

All this effort was taken to reduce the theoretical error and account for the ever in- 
creasing experimental precision. The current experimental world average for the inclusive 
branching fraction is 

B(B -» A s7 ) exp = (3.23 ± 0.42) • 10" 4 (2.69) 



combining the results of pB|-[7D|. For the theoretical prediction 

(2.70) 



of pq-[7D|. For the theoretical prediction |Jj 
B(B - X.y)% >1 . 6GeV = (3.57 ± 0.30) • lO" 4 



there is an ongoing discussion about quark mass effects. The matrix elements of Q\ 
depend at two- loop level on the mass ratio m 2 c /m^. Gambino and Misiak [71| argue that 



the running charm quark mass should be used instead of the pole mass, because the charm 
quarks in the loop are dominantly off-shell. Strictly speaking, this is a NNLO issue and 
could be used to estimate the sensitivity to NNLO corrections. Numerically, however, it 
increased the branching ratio by 11%. A more conservative error estimate would rather 
add this shift to the theoretical error. 

As b — > sj decays appear at one-loop level for the first time, they play an impor- 
tant role in indirect searches for physics beyond the Standard Model. Effects from new 
particles in the loops could easily be of the same order of magnitude than the Standard 
Model contributions. The excellent agreement of experimental measurement ( |2.69| ) and 
theoretical prediction (|2.70p therefore places severe bounds on the parameter space of 



New Physics scenarios, like multi-Higgs models [56,72], Technicolor [73], or the MSSM 
7J. 



The treatment of the matrix elements for exclusive b — > 37 decays as for example B — > 
1T*7 is in general more complicated. In this case, bound-state effects are essential and 
need to be described by non-perturbative hadronic quantities like form factors. Exactly 
these exclusive radiative decays of B mesons are the subject of this work and will be dealt 
with in detail in part II and III. 



Chapter 3 
Factorization 



The idea of factorization in hadronic decays of heavy mesons is already quite old. Let us in 
the following sections describe "naive factorization" and its extensions, state the problems 
and shortcomings of these models, and then introduce QCD factorization, which allows 
for a systematic and model-independent treatment of two-body B decays. In the last 
section of this chapter we will shortly present some other approaches used to tackle the 
difficulties with these decays. 



3.1 Naive Factorization and its Offspring 

For leptonic and semi-leptonic two-body decays, the amplitude can be factorized into the 
product of a leptonic current and the matrix element of a quark current, because gluons 
cannot connect quark and lepton currents. Pictorially spoken, the Feynman diagram falls 
apart into two simpler separate diagrams if we cut the W propagator. For non-leptonic 
decays, however, we also have non-factorizable contributions, because gluons can connect 
the two quark currents and additional diagrams can contribute. 

In highly energetic two-body decays, hadronization of the decay products takes place 
not until they have separated already. If the quarks have arranged themselves into 
colour-singlet pairs, low-energetic (soft) gluons cannot affect this arrangement (colour 
transparency) |75| , |76| . In the naive factorization approach, the matrix element of a four- 



fermion operator in a heavy-quark decay is assumed to separate ( "factorize" ) into two 
factors of matrix elements of bilinear currents with colour-singlet structure |77|]78|| , e.g. 

(D + n- 1 (c6) iV -A) (du) iV -A)\B d ) -> (tT | (du) {V - A) \0) (D + \(cb) {V -A)\B d ) 

oc UF B ^ D (3.1) 

In general, the complicated non-leptonic matrix elements (M]_M2\Qi\B) are decomposed 
into a, form factor F B ^ Ml and a meson decay constant fM 2 - The factorized matrix elements 
of Qi and Q2 are dressed with the parameters 

= CM + (3.2) 

respectively, to give the amplitude. In the literature one distinguishes three classes of non- 
leptonic two-body decays. The first class contains only a± such that the meson generated 



26 



3. Factorization 



from the colour-singlet current is charged as in ( |3.1| ). The second class involves only a-i 
and therefore consists of those decays in which the meson generated directly from the 
current is neutral. The third class finally covers decays in which the a\ and a 2 amplitudes 
interfere. 

Already two of the prominent proponents of factorization in heavy meson decays, 
Dugan and Grinstein, admitted that "at first sight, factorization is a ridiculous idea." 
|76| . In (|3.1| ) the exchange of "non-factorizable"Q gluons between the 7r~ and the (BdD + ) 
system was completely neglected, which consequently does not allow for rescattering in 
the final state and for the generation of a strong phase shift between different amplitudes. 
Yet, the main problem of naive factorization is that it reduces the renormalization scale 
dependent matrix element to the form factor and decay constant, which have a rather 
different scale dependence. This destroys the cancellation of the scale dependence in 
the amplitude and is therefore unphysical. Naive factorization cannot be correct exactly 
and gives at most an approximation for one single suitable factorization scale fif. The 
value of this particular scale is not provided by the model itself, but usually expected 
to be 0(rrib) and 0(m c ) for B and D decays, respectively. Another problem with naive 
factorization arises beyond the leading logarithmic level. Here, the Wilson coefficients 
become renormalization scheme dependent whereas the factorized matrix elements are 
renormalization scheme independent such that no cancellation of the scheme dependence 
in the amplitude can take place. This is unphysical again. 

The concept of "generalized factorization" tries to solve these problems by introducing 
non-perturbative hadronic parameters, which shall quantify the "non-factorizable" con- 
tributions and herewith cancel the scale and scheme dependence of the Wilson coefficients 
|79H3T]. Neubert and Stech |8(J replace a x and a 2 by 



,eff _ (n(„\ 1 C M\n 1 JBn,D),, A] 1 n _(„\AB-K,D), 



af = \CM + ^J[l + er^\^)} + CMer^^) (3-3) 

Due to the aforementioned renormalization scheme dependence, however, it is possible 
to find for any chosen scale fif a renormalization scheme in which the non-perturbative 
parameters S\(fif) and £s (/■*/) vanish simultaneously p2| . This leads back to naive factor- 
ization instead of describing the "non-factorizable" contributions to non-leptonic decays 
properly. Another variant of generalized factorization achieves the scale and scheme in- 
dependence via effective Wilson coefficient functions and an effective number of colours 
0: 



5 = Cf 2 + —,Cll (3.4) 



1 

Again, these effective parameters are nicely renormalization scale and scheme indepen- 
dent. Yet, Buras and Silvestrini |82| showed that they depend on the gauge and the 



infrared regulator and are thus unphysical. 

All these shortcomings are resolved in the QCD factorization approach. 



■""We put quotes on "non-factorizable" if we mean the corrections to naive factorization to avoid 
confusion with the meaning of factorization in the context of hard processes in QCD. 



3.2 QCD Factorization 
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The typical scale for a B meson decay is of order m# and therefore much larger than 
Aqcd, the long-distance scale where non-perturbative QCD takes over. QCD factoriza- 
tion, introduced by Beneke, Buchalla, Neubert, and Sachrajda, uses exactly this fact |T5| . 
This time, something is factorized in the general sense of QCD applications: Namely 
the long-distance dynamics in the matrix elements and the short-distance interactions 
that depend only on the large scale m#. And again, the short-distance contributions 
can be computed in a perturbative expansion in the strong coupling a s iO{m b )). The 
long-distance part has still to be computed non-perturbatively or determined experimen- 
tally. However, these non-perturbative parameters are mostly simpler in structure than 
the original matrix element and they are process independent. 



3.2.1 The factorization formula 

Let us here explain QCD factorization for exclusive non-leptonic two-meson B decays 
B — > M\M 2 following ||15|| . Applications of these methods to exclusive radiative decays 
will be the subject of the following two parts. We consider B — > M\M 2 in the heavy- 
quark limit and differentiate between decays into final states containing a heavy and a 
light meson H\ and L 2 or two light mesons L\ and L 2 . A meson is called "light" if 
its mass m remains finite and "heavy" if its mass scales with m& such that m/rrib stays 
fixed in the heavy-quark limit rrib — * oo. Up to power corrections of order Aqcd /n^b, for 
the transition matrix element of an operator Qi in the weak effective Hamiltonian the 
following factorization formula holds: 

(LiLalftlB) = J2 F r Ll ^ fduT^u)^ L2 {u) 

+ E F ^ 2 ( m i) fdvT{ k (v)$ Ll (v) 
k Jo 

+ [ d^dudvT t n (^ U ,v)^ B (O^M^L 2 (u) (3.5) 

JO 

(H x L 2 \Qi\B) = TFf^iml) [ duT^u) <M«) (3-6) 

3 J ° 

F?-* M denotes a B — > M form factor and <3>m is the light-cone distribution amplitude 
(LCDA) for the quark-antiquark Fock state of meson M. These non-perturbative quan- 
tities are much simpler than the original non-leptonic matrix element. The LCDA $m 
reflects universal properties of a single meson state and the form factors refer only to a 
relatively simple B — > meson transition matrix element of a local current. Both can be 
calculated using some non-perturbative technique, like lattice QCD or QCD sum rules, 
or they can be obtained from experimental results. The hard-scattering functions T^(u), 
Tl k (v) and 7^ 7/ (^,M,f) are perturbatively calculable functions of the light-cone momen- 
tum fractions u, v and £ of the quarks inside the final state mesons and the B meson, 
respectively. We distinguish "type I" or "hard vertex" and "type II" or "hard spectator" 
contributions. The idea that the decay amplitude can be expressed as a convolution of a 
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Figure 3.1: Graphical representation of the factorization formula for a B meson decaying 
into two light mesons, e.g. B~ — > tt°K~. At leading order in AQcr>/ m b there are no 
long-distance interactions between the system of B meson and the meson that picks up 
the spectator quark, and the other final state meson. 



hard-scattering factor with light-cone wave functions of the participating mesons is anal- 
ogous to more familiar applications of this method to hard exclusive reactions involving 
only light hadrons [53]. A graphical representation of ( |3.5| ) is given in Fig. |3.1| . 



When the spectator quark in the B meson goes to a heavy meson as in Bd — > D + n~, 
the spectator interaction is power suppressed in the heavy-quark limit and we arrive at 
the simpler equation (|3.6|). For the opposite situation where the spectator quark goes to a 
light meson, but the other meson is heavy, factorization does not hold, because the heavy 
meson is neither fast nor small and cannot be factorized from the B — > Mi transition. 
Annihilation topologies do not contribute at leading order in the heavy-quark expansion. 

The hard spectator interactions appear at 0(a s ) for the first time. Since at O{a s ) 
the functions Tt, k are independent of u and v, the convolution integral results in a meson 
decay constant and we see that ( |3.5| ) and ( |3.6| ) reproduce naive factorization. We now, 
however, can systematically compute radiative corrections to naive factorization. This 
immediately solves the problem of scale and scheme dependence in the conventional fac- 
torization approaches. Because the form factors are real quantities, all strong rescattering 
phases are either generated perturbatively or are power suppressed. 

In principle, QCD factorization is nothing else but a consistent formalization and 



generalization of Bjorken's colour transparency argument [75]. This is most obvious for 



the decay Bd — > D + n~ . The spectator quark and the other light degrees of freedom 
inside the B meson can easily form a D meson after the weak b —>■ c decay. This B — ► D 
transition can be parametrized by a set of form factors. The other two light quarks are 
very energetic. To form a pion they must be highly collinear and in a colour-singlet 
configuration, the probability of which is described by the leading-twist pion light-cone 
distribution amplitude. Such an energetic "colour-transparent" compact object can leave 
the decay region without interfering with the D meson formation because soft interactions 
decouple. 



3.2.2 The non-perturbative input 

We will now discuss the necessary non-perturbative input, namely form-factors and light- 
cone distribution amplitudes, in more detail. 



3.2 QCD Factorization 
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Form factors 

A form factor is a function of scalar variables accompanying the independent terms in the 
most general decomposition of the matrix element of a current consistent with Lorentz 
and gauge invariance. In the literal meaning a form factor goes with the matrix element 
of some current where initial and final state are one and the same particle. In the non- 
relativistic limit a particle has two electromagnetic form factors, which are simply the 
Fourier transforms of its charge and magnetic moment distributions. They therefore 
indeed give information on the form of the particle under consideration. 

In the context of QCD factorization we often need the matrix element of the vector 
current which is conventionally parametrized by two scalar form factors 



(P(k)\q^b\B(p)) = F^(q 2 )(p» + kn+ F^tf) ~ F™tf) ™ fl 2 ™V (3.7) 

L J q 1 

where q = p—k. Strictly speaking F£~£ p are transition form factors, which do not describe 
the form of neither the B meson nor the pseudoscalar meson P, but rather their overlap 
during the weak decay. For q 2 = the two form factors coincide, F p ^ p (0) = F p ^ p (0). 
The scaling behaviour of the form factors is 

Fjfc D (0) ~ 1 (3.8) 



vWqod) (^f ) 



3/2 



^o/+ ; w(0) ~ ^(vWqod) l-^J (3.9) 

/A \ 3 / 2 

Ffc?**® ~ l^f) ( 3 - 10 ) 



for the B — > D and B — > n transition, respectively |15| j8^| . Both the heavy-to-heavy 
and heavy-to-light B — > M form factors receive a leading contribution from soft gluon 
exchange. This is why the form factor has to enter the factorization formula as a non- 
perturbative input. 



Light-cone distribution amplitudes for light mesons 

Let us now consider the other non-perturbative ingredients which are the meson light-cone 
distribution amplitudes. We define light-cone components 

k°±k 3 , 
k± = — (3.11) 

for any four vector k^ = (k + ,k_,k±). These variables naturally distinguish between 
a particle's longitudinal and transverse degrees of freedom. Let us construct a light 
pseudoscalar meson out of the on-shell constituent quarks in a spin singlet state and with 
no net transverse momentum: 

im) = 1^^7=2 (H^ " a! A) |0) ^' l±) (3 ' 12) 

where aj* creates a (anti)quark with momentum I and spin up. Their momentum 
shall be a fraction v and v of the meson momentum k: if = vk + and if = vk + plus a 
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momentum l± perpendicular to the meson momentum direction, which adds up to zero 
for both quarks. Here and in the following we use the short-hand notation 

v = l-v (3.13) 

As long as vv = 0(1), i.e. not at the endpoints of the spectrum, l\jvv can be neglected 
compared to E 2 for a meson with large energy E. If we are interested in the leading twist 
contributions only, we can perform the l± integration and define, with fp being the meson 
decay constant, via 

d 2 U 1 . , if f 



--ip(v,l±) = 



:<$>p(v) 



(3.14) 



16?r 3 y/2N T AN 

the light-cone wave function $p. The latter is normalized as f Q dv<&p(y) = 1 and has the 
asymptotic form Qp(v) = 6vv. Using the explicit expressions for Dirac spinors we can 
then write a light pseudoscalar meson via its light-cone distribution amplitude in position 
space for light-like separations of the constituent quarks as 



(P(*)|^) a ^(0) /J |0> = ^(7fe 



la/3 



p{v) 



(3.15) 



55,36 



with i and j colour and a and (3 spinor indices. In the following we will need also the 
leading twist LCDA for a vector meson with polarization vector 77 

-1 

dv e mk z & 



<v(MM*)«fc(o)/>|o) 



f, 



v 



8N 
fyrny 
AN 



i ftapTl-Z 



(V) 



(3.16) 



Sij 



dv e 



ivk-z 



+ 



dve™ k - z gy {v \v) 



a/3 



~dve™ k - z gy {a \v) 



The light-cone wave function $y has an expansion in terms of Gegenbauer polynomials 
Cl /2 (2v-l) 



&v(l 



l+J2*nMC 3 n /2 (2v 



n=l 



(3.17) 



3x , C< 



3/2 



x 



(5x 2 — 1), etc. The Gegenbauer moments oc^(fi) are 



where C^ 2 (x) 

multiplicatively renormalized. They vanish logarithmically as the scale \x — > 00. In this 
limit, $y reduces to its asymptotic form $>y(v) = 6vv, which often is a reasonable first 
approximation. The remaining leading-twist light-cone wave functions for light vector 
mesons, $y, g^ v \ and gy^ do not contribute at leading power if the mesons are trans- 
versely polarized. An example of how to use the light-cone distribution amplitudes in an 
actual calculation will be given in the following subsection |3.2.3| . 

We conclude the discussion of light meson LCDAs with the counting rules for the wave 
functions. Using the asymptotic form &x{ v ) — 6iw for X = P and X =_L we count the 
endpoint region, where v or v is of order AqcD/ m b, as order AqcT>/ m b- Away from the 
endpoint the wave function is 0(1): 

for v, v ~ Aqcd/^6 
1 for v away from endpoint 



$x(i>) 



(3.18) 
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B meson light-cone distribution amplitudes 

The I?-meson light-cone distribution amplitude appears in the hard-spectator interaction 
term in (|3.5|) , because a hard gluon can probe the momentum distribution of b and 
spectator quark inside the B meson. Compared to the light meson LCDAs we now have 
one very heavy quark as constituent: The h quark carries the largest part of the B mesons 
momentum p: p^ = £p + w p + , whereas for the spectator quark we have l + = £p + with 
£ = 0(Aqcd/"t.&)- From the explicit expression for a heavy quark and light antiquark 
spinor we can pull out a factor (ft + 777.5)75. In the 5-meson rest frame the remainder 
is essentially proportional to f since the spectator quark is neither energetic nor heavy, 
and thus no restrictions on the components of its spinor exist. However, we still want to 
perform the l± integration and absorb it into the 5-meson light-cone wave functions. To 
this end we decompose finto a part proportional to fi + and ^_ where n is an arbitrary 
light-like vector which can be chosen for example in the direction of one of the final state 
particle momenta. Our standard choice is n = n_ = (1, 0, 0, —1). For this choice of n we 
define two scalar wave functions $m/2 in the B meson light-cone distribution amplitude 
as 0,0,^ 



(0|fc(0)«fc(*)/»|5(p)> = ^W+mJ^J ^e-^-[$ m (£)+ 7«M0] 7/J (3.19) 

This is the most general decomposition of the leading-power LCDA only if the transverse 
momentum of the spectator quark lj_ can be neglected. The wave functions describe 
the distribution of light-cone momentum fraction £ = l+/p+ of the spectator quark with 
momentum I inside the B meson. They depend on the choice of n and are normalized as 



/ de$Bi(0 = l / d^ B2 (0 = (3.20) 
Jo Jo 

Unfortunately, the 5-meson wave functions are poorly known, even theoretically. At 
scales much larger than mj they should approach a symmetric form as for light mesons. 
At scales of order m& and smaller, however, one expects the distribution to be highly 
asymmetric with £ = (9(AQCD/^fe)- 

What we will need later on besides the normalization conditions is the first negative 
moment of $m(C)> which we parametrize by a quantity \b = O(Aqcd), i-e. 

dt^f 1 = T (3-21) 

The counting for B mesons is different from the counting for light mesons. We use 
the normalization condition (|3.20|) to get 



a (t\ S «V A QCD for £ ~ A QCD /m 6 , , 

$ B!(e)~j for£~l (3 ' 22) 

This represents the fact that it is practically impossible to find the light spectator quark 
with momentum of order m&, because there is only a small probability for hard fluctuations 
that transfer a large momentum to the spectator. The same counting applies to other 
heavy mesons. 
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q i vk 
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Figure 3.2: The Feynman diagrams contributing to the 7* — > tt° + 7 form factor. 



3.2.3 A simple application: F r 



777T U 



We now want to apply the pion light-cone distribution amplitude in the calculation of 
the 7* —>■ 7r° + 7 form factor. At lowest order in the hard scattering this is a pure QED 



process where the Feynman diagrams in Fig. |T2] contribute. We assign momenta vk and 
vk to the quark and antiquark in the outgoing pion, respectively. The Feynman rules in 
momentum space give for on-shell final-state particles 



V iq) 



qj{vk)^ 



-ieQ q <f) 
+(-ieQ q j u 



i5ji 



i5 



j 1 



vk 1 - 



-ieQ„ 



qi{vk) c 



J 0a 



iQ q e ^j%q%{vk) a qj{vk)p 



vq z 



vq 2 



(3.23) 



0a 



where we neglected the quark and pion mass and simply rearranged the quark fields to 
enable the replacement 

-1 



qi(vk) a qj{vk)p 



(3.24) 



according to ( [3.151 ) for q = u or q = d. We perform the trace in fl3.23j ) and use the phase 
convention 



7T 



mm — dd 

~7T 



for the neutral pion. This leads to 

_ £>( rf ) 



V 



V2 



■ie 2 e KX ^hq x e^ u 



U 1 



1 , *«(v) 
_ , , dv — — 

3V2 -q 2 Jo vv 



which defines the neutral pion form factor F 777r o. For the asymptotic form <& n (v 
of the pion wave function we obtain its asymptotic value 



V2U — 



(3.25) 



(3.26) 



Qvv 



(3.27) 
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B„ D + 



Figure 3.3: The leading-order contribution to the hard-scattering kernel in Bd — > 
D + 7i~ . The pion's momentum is k. 

If one measures the form factor experimentally one can get information on the pions 
valence quark distribution. 

3.2.4 One-loop proof of factorization in B —* Dir 



In order to prove the validity of the factorization formulas ( |3.5|) and ( |3.6|) , one has to 
show that the hard scattering amplitudes T~, k and T* are free of infrared singularities 
to all orders in perturbation theory. For heavy-light final states such a proof exists at two 
loops fll5f . It has subsequently been extended to all orders |j89| . We here only want to 



sketch the proof of factorization for B — > Dn at the one-loop order as given in the second 
reference in |L5| . 



The leading order 

The relevant quark level decay for B^ — > D + n~ is b —>■ cud so that no penguin operators 
or diagrams can contribute. At lowest order there is only a single diagram with no hard 



gluon interactions, shown in Fig. |3.3| . The spectator quark is soft and absorbed as a soft 
quark by the recoiling D meson which is described by the B — > D form factor. The 
hard subprocess is just given by the insertion of the colour-singlet four-fermion operator 
S i n gi e t = C7 M (1 — r y 5 )bdj ll (l — 75)w. Therefore, it does not depend on the longitudinal 
momentum fraction v of the two quarks that form the emitted n~ . Consequently, with 
T^Av) being v -independent, the f-integral reduces to the normalization condition for the 
pion wave function. So (|3.6|) reproduces the result of naive factorization if we neglect 



gluon exchange. In the heavy-quark limit the decay amplitude scales as 

A(B d -> D+Tr-) ~ G F m 2 b F B ^ D ^)f^ ~ GVm'AqcD (3.28) 

Factorizable contributions 



We now have to show that radiative corrections to Fig. are either suppressed by a s 
or or are already contained in the definition of the form factor or the decay 

constant of the pion. For example the first three diagrams in Fig. [T3| are part of the form 
factor and do not contribute to the hard-scattering kernels. In particular, the leading 
contributions from the region in which the gluon is soft (first and second diagram in 
Fig. p.4|) are absorbed into the physical form factor. This is the one that appears in (|3.5|) 
and ( |3.6|) and is also the one directly measured in experiments. The fourth diagram in 
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Figure 3.4: Factorizable diagrams at d?(a s ) that are part of the form factor (first three 
diagrams) or renormalize the conserved light-quark V — A current (last diagram) . 




2 





(vl) (v2) (v3) (v4) 



Figure 3.5: "Non-factorizable" vertex corrections. 

Fig. |3.4| simply renormalizes the conserved (ud) light-quark V— A current and cancels the 
wave-function renormalization of the quarks in the emitted pion. 

"Non-factorizable" vertex corrections 

The diagrams containing gluon exchanges that do not belong to the form factor for the 
B — > D transition or the pion decay constant are called "non-factorizable" . 

The vertex corrections of Fig. [T5] violate the naive factorization ansatz ( |3.1| ). Nev- 
ertheless they are calculable and form an 0(a s ) correction to the hard-scattering kernel 
T?Ju). Although the diagrams in Fig. |3.5| exhibit seperately a collinear and soft infrared 
divergence, we can show that these divergences cancel in the sum of all four diagrams. The 
cancellation of the soft divergences is a manifestation of Bjorken's colour-transparency ar- 
gument and the collinear divergences cancel due to collinear Ward identities. So the sum 
of the four diagrams involves only hard-gluon exchange at leading power. 

At 0(a s ) the only vertex correction in — > D + n~ comes from the colour-octet 
operator oc tet = cj^(l — ^)T a b d^ ' M (l — r y^)T a u. We choose the quark and antiquark 
momentum in the pion as 

l 1 = vk + l± + t4^- h = vk-l ± + -±=n- (3.29) 

with k = E(l, 0,0,1) the pion momentum, E = pi, ■ k/rriB the pion energy, and n_ = 
(1,0,0,-1). We want to show that the transverse momentum l± can be neglected at 
leading power and that the contributions from the soft-gluon region and from gluons 
collinear to k are power suppressed. As we have to keep l± for the moment we use (|3.15|) 
before the l± integration in ( |3.14| ): 

MAOK^O^IO) = / dv f l± 3 I pfalJijsflape** (3.30) 
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By this means the diagrams in Fig. give for the O ctet insertion 

(D + n-\0 octct \B d ) 0(as) = (3.31) 
9 CV f d^q . , S7l =.l f 1 , d 2 l± ip*(vJ±) r „ , , 



where 



2p c - q- q 2 2p b - q + q 2 



Mh , h , q) = I^±^_^ii±|^ (3.33) 
2/ 2 • g + <r 2Zi - g + <r 

Here, = 7 M (1 — 75), and p& and p c are the momenta of the b and c quark, respectively. 
We now show that the integral over q does not contain infrared divergences at leading 
power in A QC D/^b- 

In the soft region all components of q become small simultaneously, which we describe 
by the scaling q ~ A. Power counting then shows that each diagram in Fig. [33] is 
logarithmically divergent. But exactly the smallness of q allows the following simplification 
of A 2 : 

A 2 (h, h, q) = ~ 4^ FM + °( X ) = °W ( 3 - 34 ) 

v k ■ q vk ■ q 

where we used that a ft, to the extreme left or right of an expression gives zero due to 
the on-shell condition for the external quark lines. So the soft contribution is of relative 
order Aq CD /m b or smaller and hence suppressed relative to the hard contribution. 
For q collinear with the light-cone momentum k of the pion we have the scaling 

q + ~ A , q ± ~ A, q~ ~ A 2 (3.35) 

so that 

A~A 4 , fc-g~A 2 , g 2 ~A 2 (3.36) 
The divergence is again logarithmic. Substituting q = ak + . . . we simplify A 2 as 

2(v + a)k x T» 2(v + a)k x T» 

A 2 (h, k, q) = — -3 T] — — oc k\ 3.37 

2/ 2 ■ q + q 2 2li ■ q + q 2 

using again the on-shell relations. To leading power, A 2 is proportional to k\ in the 
collinear region. Therefore, we obtain with k 2 = 



2ap c ■ k 2apb ■ k 

This completes the proof of the absence of infrared divergences at leading power in the 
hard-scattering kernel for B^ — > D + tt~ to one-loop order. For the proof at two-loop 
order one can similarly analyze the Feynman integrands corresponding to the diagrams 
in momentum configurations that can give rise to singularities. This proof can be found 
in [|l5j for B — > Dtt. A simpler two- loop order proof will be performed in part III for the 
one-particle irreducible contributions in B — > 77. 
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Figure 3.6: Penguin diagram and chromomagnetic dipole operator contribution to B 
MiM 2 . 



\ 



s: 



Figure 3.7: "Non-factorizable" spectator interactions. 



As | Z_i_ | <C E we can expand A 2 in \l±\/E, which amounts to simply neglecting l± 
altogether. Then li = vk and l 2 = vk, and we again can use the l± integrated pion wave 
function ( 3.14 ) to simplify the matrix element of oc t e t in (|3.31| ) to 



(D- 



7T 



\o 



octet 



IN 
ocF B ^ D 



Bd)o(a a ) 

d 4 q 



9l% I ■i^i{D + \cA x {q)b\B i )^U I dv$ v (v)ti[<y B }6A 2 (vk,vk,q)] (3.39) 



(0) / dvT ociet (v,z)$ n (v) 



(3.40) 



where z = m c jm^^ F B ^ D (0) is the form factor that parametrizes the (D + \c[. . ]b\Bd) 
matrix element, and T octet (v,z) is the hard-scattering kernel whose epxlicit expression 



can be found in 15 



Penguin diagrams 

In general, a non-leptonic B decay can also get contributions from penguin diagrams 
and the chromomagnetic dipole operator Q 8 as depicted in Fig. |3.6j They would both 
give calculable 0{a s ) corrections to the hard-scattering kernel T^(u). For Bd — > D + n~ 
however, neither of the two contributions exists because of the flavour structure of the 
final state. 



Hard spectator interaction 



In ( |3.6| ) we claimed that there is no contribution from hard-spectator interaction to heavy- 
light final states at leading power in the heavy-quark expansion. Soft contributions from 
diagrams in Fig. EO cannot be absorbed into the B — » D form factor or the pion decay 
constant. For factorization to hold true, no such soft contributions may appear at leading 
power. Although they are present in each of the two diagrams separately, to leading power 
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Figure 3.8: Annihilation topologies. 



they cancel in the sum over the two gluon attachments to the ud pair by the same colour- 
transparency argument that was applied to the "non-factorizable" vertex corrections. 
Power counting for the summed amplitude of the two diagrams shows explicitely the 
power suppression relative to (|3.28j) . 



Annihilation topologies 

Finally, we will recall the power suppression of both hard and soft part of annihilation 
topologies as in Fig. |3]^. According to (|3.22|) the light quark that goes into the D meson 
must always be soft. Therefore, the virtuality of the gluon can be at most m&AQCD- The 
potentially largest contribution can be achieved if both quarks the gluon splits into are 
soft. This configuration, however, is power suppressed via the pion wave function integral. 



3.2.5 Limitations of QCD factorization 



The factorization formulas ( |3.5[ ) and fl3.6| ) are exact in the heavy-quark limit — > oo. 
This is a theoretically very well defined limit. In the real world, on the other hand, mj 
is fixed to about 4 — 5 GeV and one may question the accuracy of the heavy-quark limit. 
The generic size of power corrections is 0(AQCD/"if>) ~ 0(10%). Yet, sometimes power 
corrections turn out to be numerically larger than suggested by this naive estimate. 

Due to a conspiracy of small parameters, like Wilson coefficients, CKM elements, or 
colour suppression factors, the leading-power, factorizable term can be somehow sup- 
pressed, making power corrections disproportionately large. A related issue is that the 
strong-interaction phases appearing at 0(a s ) in QCD factorization might receive large 
corrections from soft phases of C^qcd/" 7 ^), because a s and Aq C d/ ' m b are of comparable 
size. 

Especially for B decays into two two light mesons there sometimes is another enhance- 
ment of power-suppressed effects connected with the curious numerical fact that 

2m 2 

2/^ = 5— « 3 GeV (3.41) 

m u + m d 

is much larger than its naive scaling estimate Aqcd- The Qq contribution to the diagram in 
Fig. |3l| for Bd — > 7r + 7r~ is an example for such a chirally enhanced power correction. Such 
numerically large corrections introduce a substantial uncertainty in some decay modes. 

The last point we want to mention is the hybrid role of the charm quark. In our power 
counting we treated the charm quark as heavy, taking the heavy-quark limit for fixed 
m c /mb- However, in reality m c is not markedly large compared to Aqcd- In particular, 
the power suppression of the hard-spectator term in Q3.6| ) for B decays into a heavy and 
a light meson is Aq CD /i7i c only, if one employs the limit m h ^> m c ^> Aqcd- 
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Although the issue of power corrections still deserves thorough investigation we believe 
that the QCD factorization approach is a big step forward in the treatment of exclusive 
B decays. For the first time we really have a method based on solid theoretical grounds, 
namely the heavy-quark limit. Time will show if its phenomenological value is comparable. 
Proposals of how to estimate or even measure the up to now unknown power corrections 
are already on the market |]90|,|9Tf and more work on this subject is expected. 



3.3 Other Approaches to Non-Leptonic Decays 

We conclude this chapter with a mini review of attempts to treat hadronic matrix ele- 
ments in non-leptonic B decays. Some different formulations and generalizations of the 



naive factorization idea were presented already in section |3.1| . In these approaches a 
small set of phenomenological parameters is introduced in order to parametrize impor- 
tant non-factorizable effects. Yet, no attempt is made to calculate these parameters from 
first principles. Another class of approaches aims at a dynamical understanding of non- 
leptonic weak decays starting from QCD and making a controlled set of approximations. 
The third class, finally, is based on classifications in terms of flavour topologies or Wick 
contractions. We now will briefly discuss some representatives of the second and third 
class of approaches. 



3.3.1 Dynamical approaches 

QCD is a beautiful and most probably even the correct theory to describe strong interac- 
tion effects in weak decays of hadrons. Still, the main drawback of QCD is that it is not 
yet solved in the nonperturbative domain. So, especially for the study of non-perturbative 
properties of QCD, many more or less well-defined approximations were invented to get 
insight into QCD dynamics at least in some limiting situations. 



Large-N expansion 



QCD simplifies substantially in the limit of a large number N of colour degrees of free- 
One can treat SU(N) QCD as a theory of weakly interacting mesons, and 



dom 



an expansion in 1/N then amounts to a loop expansion in this meson theory [p3fl . Most 



important for us is that the factorization of non-leptonic decay amplitudes becomes exact 
in the large- N limit. Hence, an expansion in powers of 1/N provides a natural framework 
in which to discuss the structure of non-factorizable corrections. For kaon decays a com- 
bination of the 1/N expansion with chiral perturbation theory led to important insight 
into non-factorizable effects 
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For B decays, however, chiral perturbation theory does 
not apply. Still, one can use the large-iV limit to show that non-factorizable corrections 
to class-I B decays are generally suppressed by two powers of 1/N, whereas for class-II 
decays they have the same 1/N scaling as the leading factorizable contribution. 



Lattice gauge theory 

Although lattice field theory is the carrier of hope for many non-perturbative problems, 
no real results have been obtained so far in the evaluation of exclusive non-leptonic matrix 
elements. The Maiani- Testa No-go theorem p5[ even forbade lattice QCD to calculate 



3.3 Other Approaches to Non-Leptonic Decays 



39 



the matrix elements for non-leptonic two-body decays from first principles until Lellouch 
and Liischer found a way to possibly still do so [|9(J. Yet, there are still many technical 
difficulties with lattice QCD and non-leptonic B decays, in particular concerning final- 
state interactions and the inclusion of penguin diagrams. 

QCD sum rules 

Although intrinsically limited in their numerical accuracy, in many areas QCD sum rules 
have been established as a serious competitor to lattice gauge theory computations. In 
particular for the form factors at large recoil, QCD sum rules might be equal or even 
superior to lattice results for still some time. Some attempts were made to apply QCD 
sum rules to non-leptonic weak decays |97fl . Non-factorizable effects were studied for 
B — > J/tpKs O and B<i — > D°n° [p8|| . Recently, Khodjamirian suggested to study power 
corrections for two-body B decays in the framework of light-cone QCD sum rules |90| . An 



important conceptual limitation for QCD sum rules in non-leptonic decays is that it will 
not be possible to obtain a realistic description of final-state interactions if the projection 
on the external hadron state is performed using an ad hoc continuum subtraction. 

Large-energy effective theories 

Dugan and Grinstein tried to formalize the concept of colour transparency by introducing 
a large- energy effective theory (LEET) to describe the soft interactions of gluons with 



a pair of fast-moving quarks inside a pion [76]. Yet, due to the presence of collinear 



singularities, the LEET is not the correct low-energy theory for non-leptonic B decays. 
In the context of the LEET effective Lagrangian no light-cone distribution amplitudes can 
be introduced to absorb the factorizable collinear singularities. We refer to fl88| , |100| , pIJT[ 
for further developments of "LEET." 

PQCD 

The goal of perturbative QCD (pQCD) or other hard-scattering approaches is likewise 
the separation of soft and hard physics in the -B-decay matrix element ||102|| . Similar 
to QCD factorization the decay amplitudes are expressed as a convolution of a hard- 
scattering amplitude and meson wave functions. Yet, there are fundamental differences 
in the details. In the pQCD approach the B —>■ M form factors are assumed to be 
perturbatively calculable. It is postulated that soft contributions to the form factors are 
strongly suppressed by Sudakov effects. The form factor is then counted as being of 0(a s ) 
in perturbation theory, and thus the hierarchy of the various contributions to the decay 
amplitudes in pQCD is very different from that in QCD factorization. One consequence 
hereof is that naive factorization is not recovered in any limit. 

Especially for the claim that Sudakov effects suppress long-distance contributions suf- 
ficiently for the form factors to be calculable reliably and precisely in perturbation theory, 



serious doubts were raised ||103|| . Furthermore, the estimate of the strong-interaction phase 
within pQCD is both model dependent and numerically sensitive to effects that are poorly 
under control [0 . 
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A P PA 



Figure 3.9: The six building blocks of the diagrammatic way: Tree emission of an external 
W boson (T); internal colour-suppressed W emission (C)\ Exchange of the W boson 
between the quark lines of the initial state meson (E)\ annihilation of the B meson and 
subsequent decay of the W boson into the final state particles (A); penguin diagram with 
a "horizontal" W exchange (P); penguin annihilation (PA). 




3.3.2 The diagrammatic way and flavour symmetry approaches 

There exist several strategies to classify the various contributions to non-leptonic weak 
decay amplitudes in a convenient way, and then to use symmetry arguments to derive 
relations between different decay processes. We shall briefly discuss the diagrammatic 
approach, a parametrization introduced by Buras and Silvestrini, and a flavour SU(3) 
decomposition of decay amplitudes. 



The diagrammatic way 

Especially for the description of charmless non-leptonic two-body decays of B mesons the 



six diagrams in Fig. p79] were used [104]. Because the amplitudes A, PA, and E are sup- 



pressed with at least one factor of 1/ (mass of the decaying heavy meson) compared with 
T, C, and P, they are often neglected altogether. The fact that no explicit conntection to 
the effective theory of weak interaction exists is also a disadvantage. Instead, Feynman 
diagrams with propagators in the full theory are used. But this describes the situation 
at large scales 0(M W ) rather than at the much better suited mass scale of the decaying 
heavy meson. 



The parametrization of Buras and Silvestrini 

The basis for the parametrization of Buras and Silvestrini (BS) are Wick contractions of 
the effective Hamiltonian 'H^ =1 ||105 |. The different topologies arise from the different 
possibilities to connect the quark lines of the initial and final state mesons, taking into 
account the operator insertion. One obtains seven different topologies, which appear 
in a connected and disconnected version each. This leads to fourteen renormalization 
scale and scheme independent parameters, which allow to describe all amplitudes for non- 
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Figure 3.10: Disconnected and connected emission topology leading to the parameters 
DE and CE, respectively. The dashed lines represent insertions of the operators Qi. 

leptonic two-body decays of B mesons. An example for the simplest topology leading to 
the parameters DE and CE is given in Fig. |3.10| . The matrix elements of the operators 
Qi depend on the renormalization scale and scheme for the operators which has to be 
cancelled by a corresponding dependence of the Wilson coefficients. Usually, several 
operators are involved in this cancellation as the operators mix under renormalization. 
The fourteen effective parameters of Buras and Silvestrini are linear combinations of 
operator matrix elements and corresponding Wilson coefficients. Insertion of Q± and Q 2 
into DE and CE topologies of Fig. |3 . 10| yields the first two scale and scheme independent 
combinations of Wilson coefficients and operator matrix elements: 

E 1 = C 1 DE 1 + C 2 CE 2 

E 2 = C 1 CE 1 + C 2 DE 2 (3.42) 



which are generalizations of af c (Qi) and (i2 S (Q 2 ) in and (|3.3j ). Since all BS param- 



eters are combinations of Wilson coefficients and special contractions of local operators, 
they are in principle calculable in QCD. Unfortunately, so far the BS parametrization is in 
particular because of the large number of parameters only of limited use for phenomeno- 
logical analyses. In combination with a dynamical approach such as QCD factorization, 
however, one could obtain model-independent predictions for special decay channels or 
make statements about the relative importance of the various BS parameters. 

SU(3) flavour symmetry 

Compared with the typical scale of QCD A^j ~ 325 MeV, the masses of the light quarks 
m u ~ 5 MeV, m<2 ~ 10 MeV, and, less exactly, m s (m c ) w 130 MeV are small and approx- 
imately degenerate. Because the strong interaction is flavour blind, it therefore satisfies 
a flavour-S'C/(3) symmetry, which is explicitely broken via the inexact mass degeneracy. 
This symmetry property was used both for D- and 5-meson decays to determine for 
instance CKM elements from experimental data ||106| , |1071 . 



Let us consider the decay B — > PP as an example. We follow the conventions of ||108| . 
The Hamilton operator 7~C^ =1 connects the B mesons in the initial state i with the two 
pseudoscalar mesons P in the final state /. In view of the S'?7(3)-flavour symmetry the 
most general interaction Hamiltonian is composed of the representations R appearing in 
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f ® i. The i and / denote here both states and SU (3) representations. A completely 
reduced SU(3) amplitude can be denoted by (/||Rj||z) where I denotes the eigenvalue 
of the isospin Casimir operator I 2 . Due to Bose statistics we only have to consider the 
completely symmetric representations for a B decay into two light pseudoscalar octet 
mesons: 

(8 ® 8) 5 = 1 © 85 © 27 (3.43) 

The physical amplitude for the decay of a particle of the £77(3) representation R c with 
quantum numbers v c into two particles of the representation R a and R?, with quantum 
numbers v a and v>b, respectively, can be decomposed into the reduced matrix elements as 
follows: 

M>2: (3-44) 



:_!,<*♦*♦*>* E [ «• * « ) ( 5 * ? ) cw« 

R',u',R,u V 7 X 



where the triality T guarantees the reality of the phase factor. The quantities in round 
brackets are Clebsch-Gordan coefficients. They connect, via the Wigner-Eckart theorem, 
the irreducible tensor operators in-between two states, which are base vectors of irreducible 
representations. Because the Clebsch-Gordan matrices are orthogonal, the reduced £77(3) 
amplitudes build an orthonormal basis. In a decomposition of amplitudes, ordered by 
particle content of the final state, and then by the number of units of strangeness changed, 
as e.g. 

U B^PP — ^B^PP* B^PP yd.td) 

the matrices Os^f thus are orthogonal. This implies that the reduced SU(3) amplitudes 
Vb->pp are related with the decay amplitudes UgZ~p P via the transposed of O^^pp. If 
there is a physical reason for a reduced SU(3) amplitude to vanish, the corresponding 
combination of physical amplitudes also vanishes. For A£ = B —>■ PP decays for 
example, no matrix elements with representations 24 or 42 can contribute in the limit of 
exact £?7(3)-flavour symmetry. There are five such reduced ££7(3) amplitudes, leading to 
five amplitude relations, the simplest of them being the isospin relation 

- V2A(B+ -> 7T + 7T ) + A(B° d -> 7T + 7T-) + V2A(B° d -> 7T 7T ) = (3.46) 



Similar relations can be derived for D decays 
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Part II 

The Radiative Decays B — > V7 at 

NLO in QCD 



Chapter 4 

Basic Formulas for B — > V7 



The second part is devoted to the analysis of the exclusive radiative 5-meson decays with 
a vector meson V in the final state, i.e. in particular the decays B — > K*^ and B — > p7. 
We present a factorization formula for the B — ► V^7 matrix elements, which is based on 
the heavy-quark limit m& ^> Aqcd- This allows us to factorize perturbatively calcula- 
ble contributions from nonperturbative form factors and universal light-cone distribution 
amplitudes. Our results are complete to next-to-leading order in QCD and to leading 
power in the heavy-quark limit and include the branching ratios and estimates of direct 
CP violation, isospin and U-spin breaking effects. In chapter |] we present the structure 
of the QCD factorization formulas. These formulas are applied to B — > K*^ and B — > p7 
in chapter |5[ The results given here have to a great extent already been presented in 
|110|| . We go more into detail in the work on hand, add a proof for the factorizability of 



annihilation contributions, and, above all, include the effects of QCD penguin operators. 



4.1 What is the Problem? 



The radiative transitions b — > 37 and b — > d'y are among the most valuable probes of 
flavour physics. Proceeding at rates of order G 2 F a, they are systematically enhanced over 
other loop-induced, non-radiative, rare decays, which are proportional to G 2 F a 2 . In fact, 
the Cabibbo-favoured b — > 57 modes belong to the small number of rare decays that 
are experimentally accessible already at present. The experimental measurement of the 
inclusive branching fraction was quoted already in ( 2.69|) . The branching ratios for the 



exclusive channels have been determined by CLEO |58| and more recently also by BELLE 
7T| and BaBar fTTIJ : 



B(B° -> ir 7 ) ex P- 



(4.55 ± 0.70 ± 0.34) • 10^ 5 [68] 
(4.96 ± 0.67 ± 0.45) • 10~ 5 @ 
(4.23 ± 0.40 ± 0.22) • 10~ 5 flTTT 



B(B + fT+7) 



exp. 



(3.76 ±0.86 ±0.28) • 10~ 5 
(3.89 ±0.93 ±0.41) • 10~ 5 
(3.83 ±0.62 ±0.22) • 10~ 5 



[68] 

E3 



(4.1) 



(4.2) 
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Although for the B — > fq decays so far only upper limits exist |112j |, they are already 
very close to the Standard Model expectations: 

B(B°^p° 1 ) < 1.5 -KT 6 at90%C.L. 

B(B+ -> p+ 7 ) < 2.8 • KT 6 at 90% C.L. ^ } 

On the theoretical side, the flavour-changing-neutral-current (FCNC) reactions b — > s(d) r y 
are characterized by their high sensitivity to New Physics and by the particularly large 
impact of short- distance QCD corrections [f47] -|5"T| . 



As we already mentioned in section |2J| the inclusive b — > sj mode can be calculated 
perturbatively, whereas the treatment of the exclusive channels is in general more com- 
plicated. Seen from the experimental point of view, however, the situation is the opposite 
way round. Measurements of inclusive rates are rather difficult to perform, as one has 
to sum over all final states containing e.g. a photon and a strange quark. Furthermore, 
a lower cutoff in the photon energy has to be imposed in order to get along with the 
background coming mainly from the nonleptonic charged current processes b — > c/uqq'j. 
Theoretical input on the photon energy spectrum is then needed to translate the measured 
"kinematic" branching ratio into the total one. The exclusive decays, on the other hand, 
show a very clean experimental signature, which is particularly important for the difficult 
environment of hadron machines as the Fermilab Tevatron or the LHC at CERN. Once 
again the conservation law, that "the product of theoretical difficulty and experimental 
effort needed to solve a problem is constant," seems to hold. 

The main difficulty in calculations of exclusive hadronic decays is how to evaluate the 
matrix elements of the operators sandwiched between final V7 and initial B state. These 
matrix elements lump together all the nonperturbative long-distance dynamics below a 
scale of 0(1713)- From Lorentz covariance and gauge invariance, the B — ► V'y transition 
matrix elements have the following general form 

(V(M)7(g,e)lQi|S) = --^m B F y (0) [e^ p t^ v k x q p + i(e ■ rj k ■ q - e ■ k V ■ q)] (4.4) 

This is just a convenient parametrization in terms of the form factor Fy(q 2 = 0), which 
constitutes now the challenge to be determined. For q 2 7^ 0, i.e. an off-shell photon, 
we have two form factors: one in front of each of the two terms in ( |4.4|) , which are 
gauge invariant separately. The form factors are calculable in principle from lattice QCD 
113||. Since the lattice calculations provide results over a limited region of high q 2 , i.e. 



small recoil, extra assumptions are needed to scale down the form factors to the case of 
an onshell external photon with q 2 = 0. Via the light-cone sum-rule (LCSR) approach 
84] , 114 1 a more direct access to the q 2 = region is opened [|S6|]115|| . But the theoretical 



accuracy of sum rule calculations is, due to the duality approximation made, irreducibly 
restricted to 0(10%). There are further prescriptions on the market how to include bound 
state effects in B — > [|116| , |1171| . However, hadronic models are used which do not allow 
a clear separation of short- and long-distance dynamics and a clean distinction of model- 



dependent and model-independent features. In section |5.2| we comment on these and 
related calculations in the context of QCD factorization. 



4.2 The Factorization Formula 
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Figure 4.1: Contribution of the magnetic penguin operator Q? described by B — > V form 
factors. All possible gluon exchanges between the quark lines are included in the form 
factors and have not been drawn explicitly. 



4.2 The Factorization Formula 

In the heavy-quark limit, a systematic treatment of the hadronic matrix element is possible 



110| , |118| , |119|| . In this case, the following factorization formula is valid 



(V 7 (e)\Q t \B) 



F B ~> V (0) T 1 + I di dv 7f (£, v) $ fl (0 *v(v) 
Jo 



(4.5) 



where e is the photon polarization 4- vector. With the B — > V transition form factor F B ^ V 
and and $y, the leading twist light-cone distribution amplitudes of the B meson and 
the vector meson, respectively, we encounter universal, nonperturbative objects. They 
describe the long-distance dynamics of the matrix elements, which is factorized from 
the perturbative, short- distance interactions expressed in the hard-scattering kernels T/ 
and T? 1 . The QCD factorization formula ( f4.5|) holds for all operators in (|2.50| - pl35"|) up 



to corrections of relative order AQ CD /m b . We include the calculable contributions of 
all operators unless they are both of subleading power and suppressed by a s (fi b ) with 
fi b = 0(m b ). 

In the leading-logarithmic approximation and to leading power in the heavy-quark 
limit, Qm gives the only contribution to the amplitude of B — > V7, and the factorization 
formula ( (4.5| ) is trivial. The matrix element is simply expressed in terms of the standard 
form factor, T/ is a purely kinematical function, and the spectator term T/ 7 is absent. 
An illustration is given in Fig. [4.1|. The matrix element reads 



(V{k, r))j(q, e)\Q 7 \B) = -^rn b c v F v [e^ x % Vu k x q p + i(e • rj k • g — e • krj • q)] (4.6) 

where Cy = 1 for V = K*, p~ and cy = l/y/2 for V = p°. Our phase conventions coincide 
with those of |86| , 115|| . In particular, we have Fy > 0, e 0123 = — 1 ; and the phases of V 



(with flavour content qq') and B are such that 

(V(M)I^V?'I°> = ~^K-V,K)fv (4-7) 
(0\u W b\B(p)) = +if BP , (4.8) 

with positive fg, fy. 

Only at 0(a s ) the matrix elements of the other seven operators start contributing. 
In this case the factorization formula becomes nontrivial. The diagrams for the "type 
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Figure 4.2: 0(a s ) contribution to the hard-scattering kernels Tf from four-quark oper- 
ators Qi..s- The crosses indicate the places where the emitted photon can be attached. 
Diagrams where the gluon is exchanged in the loop vanish for the same reason as the 
mixing of Qi..a into Q-j vanishes at one-loop level: chirality conservation in the quark 
loop and QED gauge invariance. 



XlX I X X ■ X g X 

Figure 4.3: 0(a s ) contribution to the hard-scattering kernels from chromomagnetic 
penguin operator Q 8 . 



I" or "hard vertex" contributions are shown in Fig. |4.2j for Qi...q and in Fig. [O] for 
Qs- The calculation of these diagrams leads to infrared-finite functions of m^/m?,, where 
m q is the mass of the internal quarks in the loop diagram. The "type II" or "hard 
spectator" contributions are those where the I?-meson-spectator quark participates in 
the hard scattering. The non-vanishing contributions to T/ 7 are shown in Fig. |4.4|. The 



hard gluon therein probes the momentum distribution of the quarks in the B meson and 
the light vector meson. To evaluate these diagrams we need the light-cone distribution 
amplitudes (LCDA) of the participating mesons. In the case of the B meson we have to 
leading power the formula given in fl3.19| ) . If we choose n in ( |3.19|) appropriately, we need 
besides the normalization conditions only the first negative moment of which 
we parametrized by the quantity A# in ( |3.21| ). The vector meson in B — > Vy decays 
is transversly polarized to leading power. This can be described by the leading- 
twist and leading-power distribution amplitude for a light vector mesons with transverse 
polarization, defined by the first line in equation ( |3.16|) : 



(ViKvUizUMMO) = -ffiq [rf, fl^J dve^ ± (v) (4.9) 
The diagrams in Figs. [O-fO are not yet the complete set of possible contributions 



to B — > decays. One of the further possibilities is weak annihilation, depicted in 
Fig. |4~5| . If we calculate this contribution for Q1...4 using the leading-power projection 



onto the meson V in ( 4.9|) we get a vanishing result because the trace over an odd number 
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Figure 4.5: Annihilation contribution to the B — > Vj decay. The dominant mechanism 
for operators Q1...4 is the radiation of the photon from the light quark in the B meson, as 
shown. This amplitude is suppressed by one power of Aq C D/ m b, but it is still calculable 
in QCD factorization. Radiation of the photon from the remaining three quark lines is 
suppressed by (Agcr>/ m &) 2 for operators Qi..a- For operators Qn t 6, however, radiation 
of the final state vector meson quarks is power suppressed diagrammatically but gives a 
contribution from the leading power projector. 



of Dirac matrices is zero. Yet, a non-vanishing result arises from the projector 

{V(k,rj)\q lv q'\0) = fvmvVv (4-10) 

which, however, is suppressed by one power of Ag C D/ m b compared to (fOj), fl4.9|) for 
transverse polarization r\ v (f v , fy, m v ~ Aq C d, k ~ m h ). For operators Q 5j6 one gets a 
contribution from the leading power projector when the photon is emitted from the final 
state vector meson, which is power suppressed diagrammatically. Hence, weak annihila- 



tion contributions are power suppressed. When the photon is not emitted as in Fig. |4.5| , 
the annihilation contributions are either vanishing or even stronger suppressed because 
the quark propagators are scaling as l/m b instead of 1/Aq C d and only the sub leading 
projector appears.. 

Despite their power suppression, the dominant annihilation amplitudes can be com- 
puted within QCD factorization. For the explicit proof to 0(a s ) we refer to Appendix y. 
Physically, the annihilation contributions factorize because the colour-transparency argu- 



ment applies to the emitted, highly energetic vector meson in the heavy-quark limit [15 



A similar observation was already made in [ 120|| . Grinstein and Pirjol, however, claim 



that the 0(a s ) corrections to Fig. [O)] would vanish identically in the chiral limit and to 
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Figure 4.6: Other contributions that are power-suppressed in the heavy-quark limit. The 
penguin-annihilation diagrams (left) can be calculated in QCD factorization, but for the 
subleading-power Q 8 matrix elements (right) where the photon is emitted from the final 
state vector meson, factorization breaks down. 



leading- twist, because the projector in Q4.9| ) gives zero when applied to a current with an 



odd number of Dirac matrices. But, by the same argument, the diagram in Fig. ^4.5| would 
vanish even at leading order in a s , which is not the case. The proper treatment of the 
0(a s ) correction for Q 12 should employ the subleading-power projections related to the 
wave functions $y, gy^ v \ gy |85|j86[1 , corresponding to the use of (|4.10| ) at O(a s ). 



Since weak annihilation is a power correction, we will content ourselves with the lowest 
order result {O{a s )) for our estimates below. In particular, we shall use the annihilation 
effects to estimate isospin-breaking corrections in B — > Vj decays. The reason for includ- 
ing this class of power corrections is that they come with a numerical enhancement from 
the large Wilson coefficients [C\ ~ 3 1 C 7 1 ) or are numerically enhanced in case of the 
Q$fi penguin operator contributions. 

Finally, Fig. WM displays the remaining diagrams with a power suppression. The 



penguin-annihilation diagrams (left) are power suppressed in a way similar to the ordinary 
annihilation topologies. They are additionally suppressed by a s . The spectator diagrams 
from Q 8 (right) lead to amplitudes that are either manifestly power-suppressed (photon 
emitted from the quark line in the upper right) or that are superficially of leading power, 
but vanish when the leading-order projections are performed (photon emitted from either 
of the quarks forming the B meson). Recall, however, that photon emission from the 
light-quark line from Q$ is a leading-power effect (see Fig. pO]). 



4.3 The Hard Vertex Corrections 

To 0(a s ) the amplitude for B — > reads, with the short-hand notation (Qi) = 
(V^f\Qi\B) , for our standard or the CMM operator basis 

A(B^V 7 ) = ^EA( s )^ 1 (QS p) ) + C 2 (gf) + ^Q(g i )') (4.11) 

* p=u,c ^ i=3 ' 

= % E x l s) { z ^ p) ) + z ^ ] ) + E zw) 

* p=u,c ^ i=3 
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The leading-order matrix element (P7) = (Q7) is given in ( |4.6| ). At subleading order in ct ( 
the matrix elements (Qi..s,s) need to be computed from the diagrams in Figs. |4]^ - ^L4 



We have to combine results of the matrix elements obtained in different operator bases: 
ours and the one of Chetyrkin, Misiak, and Miinz (CMM). The final amplitude has to 
be independent of the choice of an operator basis, of course. Even partial results which 
have a physical meaning have to be independent of the chosen operator basis. This will 



be used subsequently. We will treat the type I contributions in the CMM (|2 . 58^pT65|) and 
the type II contributions in our operator basis ( |2.5Q - |2~57| ). 



Let us now have a closer look at the type I or hard vertex contributions of Figs. (L2 
and 13. The results for these diagrams can be inferred from [|6O-0^]. In these papers the 
diagrams were computed to obtain the matrix elements for the inclusive mode b — > 57 at 

L3| represented the virtual corrections 
§2 



In our case they determine 



next-to-leading order. In this context Figs. [4.2| and 
to the inclusive matrix elements of Qi g |6(| or Pi 

the kernels T( 6 and T-. As required for the consistency of the factorization formula, 
these corrections must be dominated by hard scales ~ and hence must be infrared 
finite. This is indeed the case. Only recently, Buras, Czarnecki, Misiak, and Urban 
completed the calculation of the NLO matrix elements including the penguin operators 
||62|| . However, they calculated in the CMM basis. Nevertheless, we can re-interprete 
their results as the perturbative hard-scattering kernels for the exclusive process. This is 
because the NLO corrections to the matrix elements have a physical effect on the inclusive 
rate. For the exclusive decay the type I corrections are completely independent of the 
type II contributions, which appear at 0(a s ) for the first time. Therefore, we can use 
different operator bases for type I and type II corrections. The results from [62] imply 



(Pi) 1 = {Qi)^r ± G i i + l 



4tt 



(4.12) 
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where Cf = (N 2 — l)/(2iV), with iV = 3 the number of colours, and 

Gi(z p ) = ^-\n-^-+ gi (z p ) 

^ < x 104, /j . , 

G 2 (^) = -— In— - + g 2 {z p ) 

r 44 i ^ a. 

27 m& 

G 4(^ c ) = ^ln — + 5- 4 (2; c ) 
81 m b 

1568 n 
G 5 = — — In — +g 5 

27 mi, 

, x 1156 , (j, , , 

G 6 (^ c ) = — ^-m hfi-el^c) 

81 m fc 

8 , A* 
G 8 = -\n^+g 8 
3 m b 

i \ 833 1 r / x 7 / \ i 10i7r 

= 972 - 4 1^)+^)] + -^- 

/ \ 833 3 r . , 20iir 

92(z) = -_ + -[ fl (*) +&(*)]- — 

598 2tt 8 v 3 3 I/1X 14ivr 

93 = ^r + 7f + 3 At -i a(l) + 2 ,,(l) + ^ 

, . 761 7r 4 T ^ 1 5,, . 37i% 

* w = -972 -wf-9- Y ' + 8 a(1) + 4 6(2) -^r 

14170 8tt 32 v _ „, 224i7r 
95 = — + + + 246(1) + — 

9 ° w = W-^-f Xs -5°W + m W + 9aW + 156W -^ 

11 2tt 2 2ivr 

where 

X 5 = f dx f dy f dvxy \n[v + x(l - x)(l - v)(l -v + vy)} « -0.1684 (4.13) 
Jo •/ o •/ o 

a(z) = - j dx j dy j dv {[2 — v + xy(2v — 3)]ln[vz + x(l — x)(l — v)(l — v + vy)} 
9 Jo Jo Jo 

43 4i7r 

+ [1 - v + xy(2t; - 1)] ln[z - ie - x(l - a;)H} + — + — (4.14) 

4 , 16 o 224 92 4ivr 

Hz) = — \nz-\ z H ^ 1 

v ' 81 27 81 243 81 

-48z 2 -64z + 4 r—— 8 2 /2 \ 2 

+ ^ y/l-Azf{z) --z 2 (-z-l) f(z) 2 (4.15) 
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1 |l/(l/ - l)x(l - x) + (2 - y)ui lnwi + (2y 2 - 2y - l)tt 2 In w 2 

* (^W 
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with Uk = y k x(l — x) + (1 — y)z and 



1 - ^l^Ai 
The results for a(l) and 6(1) read 



27T 



2i6(Az — l) arctan 



a(l) 
6(1) 



4.0859 + — 
9 



4tt 632tt 2 
+ 



320 



81 3v/3 1215 



8 

45 



d 2 lnT(x) 

dx 2 



Am 



(4.16) 



(4.17) 



4ot 



+ _ ^0.0316 + — (4.18) 



For z < 0.3 the functions a(z) and b(z) are very accurately given by their expansion in z: 
a(z) = f {[|-^-3C(3)+(|-^)ln^ + iln 2 z + 1 Lln 3 ^ 



+ 
+ 
+ 
+ 
+ 



f^-^lnz-iln 2 z+^ln 3 z 



■l-^ + 2\nz 



fln 2 z 



216 18 72 111 ^ 6 ^ 



35101 _ 35tt 2 
8640 

67801 



12 



8000 



72 

21tt 2 
20 



185 J 35 j 2 

144 111 Z 24 111 Z 



3303 i n _ _ 63 ] n 2 
800 111 Z 20 111 Z 



+17T 



(2 - £ + \ In 2 + \ In 2 * + (± - £ - ln 2 + ± In 2 

, r 3 , 5 r 4 , 49 r 5 , 231 61 \ 

-i-z -r 9 ^ ~r 72 ^ -1- 200 ^ J / 
+C(^ 7 ln 2 2) 



b(z) 



-3 + ^-ln^-^z 3 / 
+ (| + 7T 2 -21n^-iln 2 2) 



+ 


1 25 TT 

^ 12 9 


+ 


1376 , 
225 ~ l ~ 


+ 


131317 


11760 


+ 


2807617 


97200 



137 
30 



In z + 2 hfz + 

2 



2zi! 
3 



fi-ln^ + 51n 2 S + ^f 



+ 



16597 
540 



lnz + 14 In 2 ,2 + 



14?H 



+in [-2 + (1 - 2 In z)z 2 + (-¥ + $ ln z) + z 4 + jz 5 + |^ 6 ] } 
+C(^ 7 ln 2 z) 



We denote with 



m„ 



z i = 



mi 



(4.19) 



the squared ratio of the mass of the internal quark in the loop with respect to the b quark 
mass. 
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Figure 4.7: The building block for the Q1...6 contributions. The indices denote the colour 
of the quarks. The other possible topology is a sq x qb instead of the sb x qq operator 
insertion leading to diagrams with a continuous fermion line. It is equivalent to inserting 
Fierz-transformed four-quark operators into the displayed topology. 

4.4 The Hard Spectator Corrections 

We now turn to the mechanism where the spectator participates in the hard scattering. 



To find the correction for (Qi...e) we compute the left diagrams in Fig. |4.4j using the 



light-cone projectors in (|3.19f) and (fO|). We first calculate the "building blocks" B a p of 



Fig. P~T| . We obtain to leading power, with the photon on-shell and the gluon off-shell, 

B^- AMVTA) = T^^^tS (« (r' - <f ) + ^ - flgr a )) (1 - 75 )Ai 

(4.20) 

for the insertion of & (V — A) x {V — A) operator, like in Q1...4, and a (V— A) x (V + A) 
structure, needed for Qs^, respectively. The relative minus sign of the building block for 
(V— A) x (V— A) and (V—A) x (V+A) insertions can be easily understood using Furry's 
theorem. Because a loop diagram with an odd number of vector currents vanishes, it is 
the sign of the axialvector current inserted in the loop that determines the overall sign. 



The integral occuring in (|4.20D is 



Ai = 8q-r [ dx [ dy— (4.21) 

Jo Jo m q ~ r z xx — 2r ■ q xy — 10 

The 1/e pole of the building block vanishes in the naive dimensional regularization scheme. 
The use of this scheme is a priori not consistent for the topologies in figure [O] because 
traces with 75 appear. Using another regularization scheme or calculating the correspond- 
ing topologies with a continuous fermion line will give a non-vanishing pole contribution. 
This, however, is mass independent and possible finite anomaly terms emerging from it 
would be cancelled via the GIM mechanism. Therefore the result is indeed finite and we 
can safely calculate in D = 4. We checked that the building block B \ A ^ xt - VTA "> j s the 
same when calculated with the Fierz-transformed operators inserted into the other possi- 
ble topology of Fig. [17|. Apart from the fact that we set D = 4 already, kept the leading 
power part only, and an overall factor of 4 due to the definition of the operator basis, our 
building block £^ _j4 ) x ( y ~ A ) j s the same as J a p in p0| . The index a is contracted with 
the photon polarization vector e a and (3 is the gluon index. 

For the QCD penguin operators we also have to consider the insertion of the Fierz- 
transformed operators if q = s or q = b. One has the choice of inserting either both 
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possible Fierz versions of a four-quark operator into one topology of a given diagram or 
one Fierz version into the two possible topologies. The result in the end is the same. For 
the topology in Fig. |4.7| , however, we can relate the (V—A) x (V— A) and (V—A) x (V+A) 
operators via Furry's theorem, which saves oneself the calculation of the diagrams with a 
continuous fermion line. We still performed both calculations to cross-check. 

A Fierz transformation of a (V—A) x (V—A) operator again gives a (V—A) x (V—A) 
operator. Therefore, we can use B^L A > x ( v A ) f or insertions of Q3 and Q4. Yet, Fierz- 
transforming Q 5 = (sibi) V -A^l q (qjqk)v+A5ii$kj leads to 



Q 



Fierz 
5 



Siqk)s+p(qjk)s-p SuS kj 



(4.22) 



with the shorthand (qiq2)s±p '■= <7i(l i 75)^2- An insertion of such an operator (without 
the sum over q) into the building block diagrams of Fig. [O leads to 



B 



(S+P)x(S-P) 
a/3 



-m g (l + 75) 



2tt 2 



qpr a 
q ■ r 



- 9ap Ai 



(4.23) 



+ (? • rgap - r a q p - r, 7„, 7^)) Aj 



with 



Aj 



dx I dy- 



r 2 xx 



2r ■ q xy — iS 



(4.24) 



'0 JO ---q 

As the result in (4.23) contains an overall m q to leading power only the case q = b has to 
be considered wheras q = s is of subleading power. 

If we insert all B a p into the left diagrams of Fig. 44 we see that to leading power 
r = vk, the momentum of the spectator quark in the vector meson, i.e. the former light 
spectator quark of the B meson which was hit by the gluon and absorbed into the vector 
meson. We can express the result for the (Qi) again in terms of the Qj matrix element 
and have 



with 



(Qi) 



a s (Vh)C F 
4tt 



HY(z q )(Q 7 ) 



(4.25) 



v 



H. 



Hi 



Hl(z c ) 



H 



Hl(z c 



v 



27T 2 hft 

3N F v m 2 B 







H v x (z c )--H v x (l) 
2HY(1) 



dvh(v, z p )$j_(v) 



-HX(z c 



+ - 2 H v Ai) 



-Hl(z c 



(4.26) 

(4.27) 
(4.28) 

(4.29) 
(4.30) 
(4.31) 



where the factor —1/2 in H% is the ratio Qd/Qu of down-type and up-type quark charges 
and we set the masses of the light quarks u, d, and s to zero. To 0(a s ) the matrix element 
of Q2 is zero because of its colour structure. Q5 contributes to leading power to the hard 
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b. In H^Q the contributions of an up, down, and strange 



spectator scattering only for q 
quark in the loop cancel out if their masses are neglected. The first integral in (|4.26 ) is 
just parametrized by of eq. ( p.21|) and the hard-scattering function h(u,z) = Ai/u in 
the second integral is given by 



h{u,z) = — 



u—Az+ie 



+ L 2 



u—Az+is 



2 

u 



(4.32) 



L 2 is the dilogarithmic function 



L 2 (x) 



-fdt^ 
Jo t 



t) 



(4.33) 



The function h(u, z) is real for u < Az and develops an imaginary part for u > 4z. At 
small values of u it has the expansion 

,2\ 

(4.34) 



1 



h(u,z) = — 
bz 



45z 



u + 0[ — 



It is also regular for z — ► 



h(u, 0) = — 
u 



(4.35) 



and real for z — > 1. Consequently only the charm quark running in the loop leads to an 
imaginary part. The function h(v, z c ) for this case is displayed in Fig. |4.8| . 

Inserting the building block B^ +P ^ 5 P ^ of ( |4.23| ) into the diagrams in Fig. [4.4| and 
performing the projections onto the B meson and vector meson wave function gives to 
leading power a non- vanishing result only for the Ai term in (|4.23 ). This insertion con- 
tributes due to the colour structure only to H§ (and not to Hq ), which in turn can be 
expressed via HY as done in ( |4.30| ) . 

The correction to (Qs) from the hard spectator interaction comes from the second 



diagram in Fig. |4.4j. One finds 



(Qs) 



11 



(Qr f MCF HY 



where 



H, 



v 



4tt 2 f B h 



47T 



dv 



(4.36) 



(4.37) 



3A^ F v m 2 B J £ Jo v 

Here, the hard-scattering function is simply 1/v such that the second integral can be 
performed analytically, leading to the combination 3(1 — a\ + + ■ ■ ■) of Gegenbauer 
moments in (|3.17|) . 

Because the hard spectator corrections are not proportional to the form factor, we had 
to remove in H\ % {z) and H% the Fy contained in (Qr). The gluons in Fig. |4.4| transfer a 



momentum of order fih ~ v/ '^QCD^nb- Therefore, we set a s = a s (/i/ l ) in ( f4.25| ) and ( f4.36p 

0(m b ). 



For our numerical analysis we shall use ^ = v^Wi with = 0.5 GeV and \x 

Finally, we can combine the results obtained so far and write, summing the up- and 
charm-quark contribution 

Gp \ ^ 



A(B -> \/ 7 ) 



X^a^iVjlQrlB) 



(4.38) 
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Figure 4.8: The hard-scattering kernel h(v, z c ), for a charm quark running in the loop, as 
a function of v. 



where at NLO 



a p 7 (Vj) = C 7 + 



47T 
47T 



,i=l,2 j=3...6,8 / 

J] C s {ji h )Hj\ 

K i=3...6,8 / 



(4.39) 



Here, the NLO expression for C 7 has to be used while the leading order values are sufficient 
for the other Wilson coefficients. The explicit formulas for the Wilson coefficients can be 
found in f|l],fG| and in Appendix |A|. 

The scale dependence of the matrix element of Qr is reflected in the running of the 
product of 6-quark mass and form factor, which is explicitely given as 



(////,• /V)M • (»,,,■ /-V )[//>,,] (l-^8C F lnA 



(4.40) 



This dependence on the scale u. has to be taken into account when the residual scale 
dependence of physical quantities is investigated. 



4.5 The Annihilation Contribution 



So far, the flavour of the light quark in the B meson and therefore the B meson charge were 
irrelevant. Yet, if we consider annihilation contributions we become sensitive to whether 
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we are dealing with the decay of a B° or B meson and thus to isospin breaking effects. 
In [ITDf we considered the annihilation contributions from operators Qi,2 only. They are 



CKM suppressed for B — > K*^, but give a sizeable contribution to B — > pj decays. 
Kagan and Neubert found that in particular the operator Q 6 gives a large contribution 



to B — > K*j because for it the CKM suppression is lifted [121]. So in the following we 



will calculate annihilation contributions from Qi...6- O n the other hand, we decided not 



to include the subleading power contributions in Fig. [4.6| because they are additionally 
suppressed by a s or not calculable. We will further comment on this in the following 
section. 

If the photon emission is from the light quark in the B meson the annihilation ampli- 
tude contains 

h y = 2^ f B m V f V 

F v m B m b X B 

whereas the Q^q insertion with the photon emitted from one of the vector meson con- 
stituent quarks leads to 

€ ^Mt fdv (442) 



v 



Recalling that F p ~ , f B ~ in the heavy-quark limit, we note that b v ,d v ~ 

This shows explicitly the power suppression of weak annihilation. The integral 
in Q4.42D again can be performed and leads to 3(1 =F + a 2 + • • •) for f and v in the 



denominator, respectively. The generic $5,6 annihilation matrix elements therefore have 
the opposite sign and an approximately twice as large absolute value than the ones from 

Ql..A- 

The annihilation amplitude for B~ — > p~7 for example is 



G F 

A ann (B -> p 7) = -j=Qu 



X^aibf +(Ai d) + A^) ( ad> p +%-a s d p v + a 6 d p . 



(p-i\Q 7 \B-) 



(4.43) 

where the dj denote the following combinations of Wilson coefficients, evaluated in leading 
logarithmic approximation 

ai, 2 = C li2 + — C 2 ,i (4.44) 

a 4 = C 4 + ^C 3 (4.45) 

a 6 = Ce + ^Cs (4.46) 

The annihilation components are included in the decay amplitudes by substituting 

a? - a u 7 +a u ann (4.47) 
a c 7 -> a c 7 + a c ann (4.48) 



where 



re (iT 7 ) = Q d [a 4 b K * + a 6 (df + df)] (4.49) 
AK*--y) = Qu[aib K * +a 4 b K * +Q s /Q u a 6 d«* + a e df] (4.50) 
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Cn(A) 


= [~a 2 & p + a 4 & p + a 6 « + 4)] 


(4.51) 


<Cn(P"7) 


= Q u [ai6 p + a 4 b p + Q s / Q u a 6 d p + a 6 <] 


(4.52) 




= Q4a 4 6^ + a 6 (df + df )] 


(4.53) 


<nn(^*"7) 


= Q u [a±b K * + Q s /Q u a 6 df* + a 6 df*] 


(4.54) 


<nn(A) 


= Q d [a 4 b p + a 6 (d p + d P ] )} 


(4.55) 


<nn(P"7) 


= Q u [a<ib p + QJ Q u a 6 d p + a 6 d£] 


(4.56) 



In B° — > p°7 the minus sign of a<i comes from the relative sign between the up-quark and 
down-quark components of the p° wave function. It is the up quarks only that produce 
the p° in the annihilation process while only the down quarks are relevant in (Q7). 

We include annihilation contributions in our numerical analysis unless otherwise 
stated. 



4.6 A Comment on Power Corrections 

Here, we want to comment on the role of power corrections in B — > V'y decays. Let us 
first repeat that the annihilation effect from operator Qi gives a numerically important 
power correction because it comes with a relative enhancement factor of \C\/C-j\ ~ 3. 
This leads to a 40% correction in the up- and a 4% correction in the total amplitude of 
the charged mode B~ — ► p~7- Note that the parameter Q u b v , describing the generic effect 
of the Q1...4 annihilation terms, is 18%, which is consistent with a Agco/mb correction 
of canonical size. The contribution Q u d v from Q$fi, however, can be as large as 40% 
which is larger than naively expected. But it is accompanied by the very small Wilson 
coefficient combination m —0.03. 



The power-suppressed penguin annihilation diagrams of Fig. [4.6| can be calculated 
within QCD factorization, but they contain an a s {0{m^j) and hence their effect is safely 
negligible. For the same reason we do not include QCD corrections to the annihilation 



diagrams in Fig. (O) 



The subleading power contributions for the matrix element of the chromomagnetic 
dipole operator Q$ deserve some special inspection. For one of them factorization breaks 
down. Kagan and Neubert already calculated the diagram in Fig. |41] where the photon is 
emitted from the vector meson spectator [ |121|| . This diagram is explicitely suppressed by 
power counting because both the gluon and quark propagator are hard. So it is sufficient 
to consider the projection onto the leading-power part of the vector meson wave function. 
We arrive at a result containing the integral 

a 8 (ji) / dv^p- (4.57) 
Jo v 

which suffers from a logarithmic divergence as v — *■ 1 if *&y(v) vanishes only linearly at 
the endpoints. Exactly such infrared divergences signal the breakdown of factorization at 
subleading power. As this contribution is dominated by soft physics we cannot bother the 
not necessarily small a s in (|4.57|) to simply neglect it. Although the endpoint singularity 
can be regulated by introducing a cutoff in the convolution integral [|121|| we refrain from 



including this contributions in our phenomenological analysis. The reason is that no 
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explicit numerical enhancement, as for the annihilation contributions, is present and that 
there are other sub-leading power corrections which are not (yet?) calculable. These we 
discuss next. 

At the same level of power counting we also have contributions from the photon emis- 
sion from the other three possible quark lines in the Q$ diagram. As already mentioned 
in section fO|, diagrammatically these contributions are of leading power. Yet, the lead- 



ing order projection gives a non-vanishing result only for the contribution in Fig. 4.4 



But if we start to consider the contribution from the subleading power diagram with the 
projection onto the leading power wave function, we also have to treat the leading power 
diagrams with projections onto subleading power LCDAs. These are calculable within 
QCD factorization, but the Qg diagram of Fig. f4.4j involves not only the first negative 
moment of $#1 but also that of $#2- This is unknown and would have to be parametrized 
by a further non-perturbative quantity Xb2- But there are contributions which are even 
less known. Among these are the contributions from keeping the transverse momentum 
l± of the spectator quark inside the B meson. This would require to modify the 5-meson 
wave functions in (|3.19|) . Yet, compared with the annihilation contributions, the power 
corrections from Qg are numerically absolutely insignificant. Therefore we will not include 
them in our numerical analysis in the following chapter. 

In contrast to nonleptonic modes with pseudoscalar mesons in the final state, no 
chirally enhanced power corrections from the light-cone expansion of meson wave functions 
arise in our case. Finally, power corrections can also come from the loop effects with up- 
and charm quarks whose leading-power contributions were computed in ( |4.25|) . These 
power corrections correspond to the region of integration where the gluon becomes soft, 
that is v = 0(AQCD/ m b)- Their contribution is nonperturbative and cannot be calculated 
in the hard-scattering formalism. Nevertheless, the expression in ( f4.25|) can be used to 
read off the scaling behaviour of these power corrections in the heavy-quark limit. For the 
charm loop the kernel approaches a constant ~ vn^jvn? c in the endpoint region as shown 
in ( |4.34j ). Taking into account the linear endpoint suppression of the wave function $y, 
the integral in (f4.25|) over the region v ~ hQCD/ m b therefore contributes a term of order 
(Agc^/rrtf,) 2 x (mb/m c ) 2 = (Aq C D/ m c) 2 - It is interesting that we thus recover the power 
behaviour of soft contributions in the charm sector first pointed out in p5 . This was 



discussed for the inclusive decay b — > in |65H67| and for the exclusive mode B — > K* r y 
in | 122 |. Numerically, this correction is with ~ 3% in the decay rate very small. A 
similar consideration applies to the up-quark sector. In this case the endpoint behaviour 
of the kernel is singular Q4.35| ), which now leads to a linear power suppression of the form 
AQ C D/ m b- This coincides with the scaling behaviour derived in |67j in the context of the 
inclusive process. 

In view of these remarks, we see no indication for large uncalculable power corrections. 
Further studies on this issue are definitely desireable, but would go beyond the scope of 
this work. 



4.7 The Observables 

In this section we present the formulas for the observables of interest in the phenomenology 
of B — > V'y decays. These are the branching ratio, CP and isospin-breaking asymmetries, 
and U-spin breaking effects. 
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From the amplitude in ( J4.38 ) the branching ratio is obtained as 



B(B -> V-y) = r B 



G 2 F a m^ml 
32^ 



o?(V 7 ) 



(4.58 



where g = s for V = K* and g = d for V = p. 

The rate for the CP-conjugated mode B — > is obtained by replacing A 
We may then consider the CP asymmetry 



p 



(q)* 



A C p(V-f) 



T(B -> p 7 ) - T(B -> p 7 ) 

r( J B^p 7 ) + r(s^p 7 ) 



A 



(4.59) 



with T denoting the rate. A non-vanishing CP asymmetry appears at 0(a s ) only. Ex- 
panding Acp in ot s we obtain 



A C p(V-f) 



47T 



(9) 1 2 
1 

lev! 



|C 7 



(4.60) 



Zi Im (Gi^c) - (7,(0)) + Ci Im (ffft^) - fff(0)) 



8=1,2 

2rp\ 2 
2?7 



(1 -p) 2 + ?f 



for V 
for V 



K* 



P 



(4.61) 



where we neglected annihilation contributions and used only one a s . When expanding 
in a s we in particular considered the denominator in ( |4.59| ) to 0{a Q s ) only because the 
numerator starts with 0(a s ). In this approximation the effects from penguin operators 
cancel out completely in ( |4.61|) . We see that the CP asymmetry in B — > K*j decays is 



,Wm(s)|2 



suppressed to the one in B — > fry via the factor -\\r/K 
—0.04 and is thus very small. 

Further interesting observables are the isospin-breaking quantities 

4r(fi + 



A 2 [(1 - p 2 ) 2 + 77 



r,-l 



V 



A 



A y 



A(V 7 ) 



r(B° 

clT{B~ 



+ y- 7 ) 



- 1 



T(B° -> ^° 7 ) 



A v 4- A y 



(4.62) 
(4.63) 
(4.64) 



Within our approximations, isospin breaking is generated by weak annihilation. Isospin 
breaking was already discussed in [ [I23|| , partially including NLO corrections. It was 
considered in more detail for B — > K* / y in [|121|| . 

Another application of our results concerns an estimate of U-spin breaking effects in 
B — > V'j decays. U-spin symmetry, the symmetry of strong interactions under exchange 
of d and s quarks, has been advocated as a tool to control hadronic uncertainties in tests 
of the Standard Model |124l|T25| . Defining 



AB{B -> iT 7 ) 
AB(B -> p 7 ) 



B(B + 
B(B + 



K* + 7 ) - B(B- 
p + 7 ) - B(B~ 



- K*" 7 ) 
P~ 7 ) 



(4.65) 
(4.66) 
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the quantity 

AB(B -> ir* 7 ) + AB(B -> p 7 ) = (4.67) 

in the limit of U-spin symmetry and within the Standard Model. This has been discussed 
125,126] and was considered in more ail in |127| . 



in 



Because all the contributions to B — > V'y are proportional to the matrix element of Q-j 
in ( |4.6|) , the polarization of the final state photon is left-handed for a B and right-handed 
for a B decay. Small corrections of 0(m s jm^) come from the m s (l — 75) contribution 
in the operator Q7 which we neglected. In many extensions of the Standard Model the 
(1 — 75) part in Q 7 can be considerably enhanced due to a chirality flip along a heavy 
fermion line in the electroweak loop process. If the New Physics contribution appears not 
already at tree level, the only effect is a new local Q 7 operator with so^il — 75)6 F^ v 
structure and a corresponding Wilson coefficient C 7 . Such an operator would not affect 
our result for a 7 in (|4.39|) but simply give an additional contribution C' 7 (Q' 7 ) to ( f4.38|) . 



Here the Dirac structure of (Q 7 ) differs from the one given in ( |4.6| ) by a relative minus 
sign in front of the i(e ■ rj k ■ q — e ■ k rj ■ q) term. 



A real photon can convert inside the detector via the Bethe-Heitler process ||1 28|| into 
an electron-positron pair. Any significant measurement of an angular dependence of this 
conversion rate will be an indication for new physics. Yet, the experimental efficiency for 
this measurement is only at the few percent level ||129|| . 



Chapter 5 

Phenomenology of B — > V7 



Here we want to evaluate the results of the previous chapter numerically for the B 
and B — > p7 decays. Our choice of input parameters is summarized in table [TT 
default choice for the CKM angle 7 is 58°. 



The 



5.1 Numerical results 

We begin with the numerical result for the NLO QCD coefficients a 7 (K*'~f) as a typical 
example. For central values of all parameters, at \i — nib, and displaying separately the 
size of the various correction terms without annihilation effects, we find 



4(iT7) 



a?(iT 7 ) 



+0.0113 



-0.3221 

C 7 LO 

-0.4070 - 0.0257* 



AC? LO 



+0.0113 



-0.3221 

C 7 LO 

-0.4167 - 0.0684i 



-0.0796 - 0.0134i 
T 7 -contribution 



-0.1397 - 0.0682i 
T 7 -contribution 



-0.0166 - 0.0123i 
T 77 -contribution 



+0.0338 - 0.0002i 
T^-contribution 



(5.1) 



(5.2) 



We note a sizable enhancement of the leading order value, dominated by the T^-type 
correction. This feature was already observed in the context of the inclusive case in 
pD| . A complex phase is generated at NLO where the hard- vertex corrections T 1 and 
the hard-spectator interactions T n yield comparable effects for the charm amplitude. 
Although the hard-spectator interaction HY{0) for a massless up quark running in the 
loop is real, a tiny imaginary part is generated from the QCD penguin operators Q^q 
also in the up amplitude. The numerical results for the a^(pj) are very similar. Then, 
with the two equations above and ( |4.60J). w e can immediately estimate the CP asymme- 



tries A^p^— * (K*^) w —0.5% and A^j}— -* (pj) ~ +13% with annihilation contributions 
neglected. The strict expansion in a s with the neglect of annihilation contributions and 
NLO corrections to the denominator of the CP asymmetry overestimates its numerical 
value. 

For B —>■ K*^j the effect from the annihilation contributions to the charm-quark part 
of the amplitude is numerically small because only the penguin operators with tiny Wilson 
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CKM parameters and coupling constants 


V us 


v cb 




V u b/V C b\ 


MS 


a 




0.22 


0.041 


0.085 ±0.025 


(225 ± 25) MeV 


1/137 


1.166 x 10~ 5 GeV~ 2 


Parameters related to the B mesons 


m B 






T B + 




5.28 GeV 


(200±30) MeV 


(350 ± 150) MeV 


1.65 ps 


1.55 ps 


Parameters related to the K* meson 


86 




F K * 




m K * 


af 


af 


fx* 


0.38 ±0.06 


185 MeV 


894 MeV 


0.2 


0.04 


230 MeV 


Parameters related to the p meson ' 


36 




F P 


J P 


m p 


a p x 




f P 


0.29 ±0.04 


160 MeV 


770 MeV 





0.2 


200 MeV 


Quark and W-boson masses 


m b (m h ) 


m c (m b ) 


"\pole 


M w 


(4.2 ± 0.2) GeV 


(1.3 ±0.2) GeV 


174 GeV 


80.4 GeV 



Table 5.1: Summary of input parameters. 



coefficients can contribute: 

'"Vi = -- 
3 



,.(^*°7) 



<nn(^*~7) = ±3 



-0.0063 ±0.0138 ±0.0203 

a^) K * a 6 d K * a 6 d£* 

-0.0063 -0.0069 ±0.0203 

a^b K Q s /QuCted K a&d^ 



-0.0093 



±0.0047 



(5.3) 
(5.4) 



But as they come with opposite sign, they still lead to a 3.4% difference on the am- 
plitude level which results in a —7.5% isospin asymmetry. They also suffice to make 
the expectation for the branching ratio B(B° — > K* Qr y) = 7 AO ■ 10 -5 larger than that 
of B(B~ — > K*~j) = 7.25 • 10" 5 although t b o < t b +. The operators Qi,2 with large 
Wilson coefficients contribute to annihilation only in the up-quark sector. We have 
ai b K * = 0.2781 which gives a% nn (K*-j) = ±0.1901. This represents 45% of the abso- 
lute value of a"(i^*7). Yet, the strong CKM suppression |A!^ / Mc S ' ) | ~ 0.02 puts this large 
correction for B — > K*j into perspective. If we consider the isospin breaking due to Qi,2 
only, we get A(K*j) = —0.9%. For B — > fry the numerical values for the a^ nn are similar, 
but as both sectors of the effective Hamiltonian have the same order of magnitude the 
impact of annihilation contributions is larger. 

We now turn to analyze the implications of the numerics at amplitude level for our 
observables of interest. The net enhancement of aj at NLO leads to a corresponding 
enhancement of the branching ratios, for fixed value of the form factor. This is illustrated 
in Fig. |57i] where we show the residual scale dependence for B(B — > K*°^) and B(B~ — > 
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Figure 5.1: Dependence of the branching fractions B(B° — > K* 0r y) and B(B~ — > p~ r y) on 
the renormalization scale p. The dotted line shows the LO result, the dash-dotted line 
the NLO result including type-I corrections only, and the solid line shows the complete 
NLO result for central values of all input parameters. For these three results the effect of 
annihilation contributions was included. The dashed line, finally, shows the NLO result 
without annihilation contributions. 



p~7) at leading and next-to-leading order. As expected, the inclusion of the hard-vertex 
corrections T 1 reduces the scale dependence coming from the Wilson coefficients. The 
scale dependence of the complete NLO result, however, is "deteriorated" because the type- 
II contributions appear at 0(a s ) for the first time and therefore introduce a completely 
new scale dependence which would be cancelled at NNLO only. The effect of annihilation 
contributions on the branching ratios is small for our default choice of the CKM angle 
7 = 58° 

The sensitivity of the B —>■ K* 0/ y and B~ — > p~ r y branching ratios, and of the CP 
asymmetry Acp{p ±r )) to variations in the relevant input parameters is summarized in 



table |57!| Numerical differences compared to the results presented in table 2 of ||110|| are 
due to a change in the input value for the B meson decay constant from fg = 180 MeV 
to f B = 200 MeV (+0.8% for the K*° and +2.5% for the p~ mode branching ratio) and 
the inclusion of QCD penguin operator contributions (+3.6% and —2.4%, respectively). 
The uncertainty of the branching fractions is currently dominated by the form factors 
Fk*, F p . They are taken from QCD sum rules [^6| with an error of 15% and enter the 



branching ratio quadratically. This situation, however, can be systematically improved. 
In particular, our approach allows for a consistent perturbative matching of the form 
factor to the short- distance part of the amplitude. For B(B~ — > p~7) a considerable 
effect also comes from the variation of the hadronic parameter A# and the CKM angle 7. 
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Jd{Jd — > A 7J [1U J 


ByB -^/) 7) [10 J 


a i „±„,\ ro/i 


central 


7.40 


1.61 


10.11 


Fk* 


+2.35/-2.03 


— 


— 


F P 


— 


+0.42/-0.37 


+0.08/-0.06 


nib/2 < fi < 2m b 


+0.85/-1.13 


+0.05/-0.18 


+5.65/-2.51 


m b 


+0.57/-0.55 


+0.12/-0.12 


+0.16/-0.14 


Xb 


+0.40/-0.15 


+0.39/-0.13 


+0.31/-0.14 


m c 


+0.38/-0.44 


+0.11/-0.11 


+1.21/-1.29 


A (5) 

A JTS 


+0.22/-0.23 


+0.04/-0.04 


+0.17/-0.19 


\v ub /v cb \ 


+0.05/-0.05 


+0.15/-0.13 


+3.54/-3.33 


7 = (58 ± 24)° 


+0.08/-0.11 


+0.34/-0.22 


+0.12/-2.20 


fs 


+0.14/-0.14 


+0.06/-0.06 


+0.10/-0.05 


m b <-> TUb 


+0.24/-0.66 


+0.11/-0.14 


+0.00/-0.81 



Table 5.2: Predictions for branching ratios and CP asymmetries with the errors from the 
individual input uncertainties. 

The vector meson decay constants fy and the Gegenbauer moments af are not known 
very precisely either, but they appear always in combination with Fy and/or Xb, which 
are both generously varied. The distinction between m& and ms in Gi, H^, b v , and d v 
is, strictly speaking, a subleading power issue. An estimate of this effect is given in the 
last row of our output tables. 



Adding the uncertainties for the branching ratios in table |5.2| in quadrature leads to an 
error estimate of (9(35%) and 0(40%) for the B — > K*^ and B fry mode, respectively. 
Taking this sizable uncertainties into account, the results for B — > K*^ are compatible 
with the experimental measurements quoted in (|4.1|) and ( |4.2| ), even though the central 
theoretical values appear to be somewhat high. According to our predictions we expect 
the B —* fry decays to be seen soon at the B factories as their experimental upper limits 
on the branching ratios quoted in (|4.3j ) are only a factor of two above our results (we get 
B(B° p° 7 ) = 0.76 • 10~ 6 ). 

If these measurements of B(B — > py) are available one can also consider the ratio of 
B — > K*^( and B — >• fry branching ratios where for example of (|4.58| ) cancels. Fur- 
thermore, the ratio of form factors is better known because it are mostly SU(3) breaking 
effects that make it differ from 1. The ratio of the rates can also be employed to get 
constraints on the p — f] plane of the CKM unitarity triangle. 

The direct CP asymmetry, which is substantial for the fry modes, is much less de- 
pendent on the form factors because the \Fy\ 2 cancels. Here, the largest theoretical 
uncertainty comes from the scale dependence. This is to be expected because the direct 
CP asymmetry is proportional to the perturbative strong phase difference, which arises at 
0(a s ). If one wants to reduce the renormalization scale dependence for the CP asymme- 
try one could determine its NLL result. This would amount to calculating the absorptive 
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Figure 5.2: The CP asymmetry Acp{p ±f y) as a function of the CKM angle 7 for three 
values of the renormalization scale p, = rrib/2, rrib and 2m;,. 



parts of the diagrams in figures [4.2| - [4.4| with one additional gluon. But as the influence 
of power corrections on the CP asymmetry might be large as well, we see at present no 
urgent motivation to embark upon such a three-loop calculation. 

The decay B —>■ fry also depends sensitively on fundamental CKM parameters, such 
as I V^i,/V^b| and 7, and can thus in principle serve to constrain the latter quantities once 
measurements become available. This is further illustrated in Figs. |5.2| and |5.3| where the 



dependence on 7 is shown for Acp(p 7) and B(B~ — > p~7), respectively. Including 
or omitting penguin operators has practically no effect on the CP asymmetries as was 
mentioned already after equation ( [4.6 1| ). The Acp(p ±r y) gets smaller by about 20% if 
annihilation contributions are omitted whereas Acp{p 0, j) = 8.25% (for the CKM angle 
7 = 58°) is practically not affected over the whole range of 7. The maximal value of 
the CP asymmetry in Fig. is obtained for 7 rs 65°, near the value of 7 prefered by 
standard unitarity triangle fits. For B — > K*^ the CP asymmetry is negative and smaller 
by a factor |A t W /*t°l « 25. We get typically Acp(K*°i) = -0.3%. 

Let us now comment on Fig. |5.3| . The 7 dependence of the leading and next-to-leading 
order B~ — > p~7 branching ratios can be understood when looking at the amplitudes 
a"(p7) = —0.4154 — 0.0685i and cC nn (p~7) = 0.1884, which are accompanied by A« = 
I Kid I iKifele -27 . The leading order value for without annihilation effects is simply Cj° = 
—0.3221. Its absolute value gets reduced when the positive Oann(P~7) i s added and hereby 
also the dependence on 7 gets reduced. The same story is repeated at next-to-leading 
order. Among the a" amplitudes the NLO has the largest absolute value and thus 
the largest 7 dependence, which is again reduced by the inclusion of the annihilation 
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Figure 5.3: The branching fraction £?(.B — > p 7) as a function of the CKM angle 7 at 
leading and next-to-leading order with and without annihilation contributions. 
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Figure 5.4: The isospin-breaking asymmetries A(K*j) and A(fry) as a function of the 
CKM angle 7 at leading and next-to-leading order with and without the inclusion of QCD 
penguin operator effects. 



contribution. For our default choice of 58° for 7 we see that the effect of annihilation 
contributions on the branching ratio is small. Although for this 7 the real part of the 
total amplitude is increased when annihilation is included, this is compensated by the 
imaginary part getting smaller. 

Figure |5.4| finally displays the isospin breaking quantities A(K #r y) and A(/ry) as a 
function of the CKM angle 7. Within our approximations isospin breaking is generated 
by weak annihilation. For comparison with our earlier work [110] the dashed curves give 
the result without QCD penguin operator effects. As Kagan and Neubert ||121|| already 
realised, this has large impact on the prediction for the isospin asymmetry in B — > K*^. 
We comment briefly on the numerical results of Kagan and Neubert in the next section. 
The QCD penguin operators give a, albeit small, contribution to the not CKM suppressed 
charm amplitude in contrast to Qi 2 , which can contribute to the up amplitude only. 
Looking at equations fl5T3] ) and ( p.4j ) we see that the dominant effect is indeed coming 
from Qq. 

In table [O] we summarize the sensitivity of A(K*j) and A(fry) with and without 
the effects from QCD penguin operators to variations of the relevant input parameters. 
Especially the strong 7 dependence of A(jry) around 7 = 58° is striking. It is also rather 
insensitive to the inclusion of penguin-operator and even next-to-leading order corrections 
(see Fig. |5.4j) . Therefore the isospin asymmetry in the p7 channel, which is far easier to 
measure than the CP asymmetry, could give valuable direct information on the CKM 
angle 7. Yet, one has to keep in mind that the errors from uncertainties in other input 
parameters, in particular Ab, are sizeable as well. 

Averaging the branching ratio measurements in (|4.1|) and (|4.2|) gives B(B° — > K* 0/ y) = 
(AAA ± 0.35) ■ 10~ 5 and B(B + -> K*+i) = (3.82 ± 0.47) ■ 10" 5 . This implies, when 
corrected for the B- meson lifetime differences, an isospin violation of A(i^*7) exp ' = 
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a ^ i^^a ro/i 
A(A 7J [/oj 


A no pcn {K 7 J [/oj 


A(p7j [/oj 


A no pcn^P 7 J [/oj 


central 


-7.46 


-0.92 


2.02 


10.04 


Fk* 


+0.94/-1.26 


+0.12/-0.16 


— 


— 


F P 


— 


— 


+1.23/-0.79 


+2.54/-1.77 




+2.69/-5.29 


+0.05/-0.13 


+2.41/-4.37 


+1.89/-0.81 


m b 


+0.46/-0.52 


+0.04/-0.04 


+0.01/-0.02 


+0.60/-0.54 


Xb 


+1.71/-0.70 


+0.27/-0.63 


+15.56/-4.92 


+13.65/-4.13 


m c 


+0.20/-0.24 


+0.03/-0.03 


+0.19/-0.22 


+0.03/-0.04 


A (5) 

A-— 

MS 


+0.14/-0.13 


+0.01/-0.01 


+0.31/-0.28 


+0.14/-0.12 


\V u b/V c b\ 


+0.26/-0.26 


+0.27/-0.27 


+0.03/-1.62 


+0.03/-1.76 


7 = (58 ± 24)° 


+0.62/-0.46 


+0.66/-0.49 


+21.79/-14.10 


+23.25/-15.04 


f B 


+1.04/-1.01 


+0.13/-0.13 


+1.15/-0.96 


+2.37/-2.17 


rrib <-> me 


+2.08/-1.72 


+0.17/-0.22 


+2.34/-0.29 


+4.18/-2.45 



Table 5.3: Predictions for the isospin-breaking asymmetries A(K*j) and A(p7), including 
and not including Q CD-penguin-operator contributions, with the errors from the individ- 
ual input uncertainties. 

(—19.2 ± 11.8)%. It is not significantly different from zero and consistent with the value 
A(K*^y) = (— 7.46^5;gg) % from table |5.3j when adding all uncertainties in quadrature. Our 
corresponding prediction for the p7 mode is A(p7) = (2.02^^;gg)%. 

The last application of our results we want to discuss is an estimate of U-spin breaking 
effects in B — > V'y decays. Using our expressions and central values for all parameters 
we find for the differences of the charge conjugated branching ratios defined in ( |4.65|) and 

AB(B — > K*q) = -7.3 -10- 7 (5.5) 
AB(B -> p 7 ) = +4.3 • 10~ 7 (5.6) 

where we have chosen the CKM angle 7 = 7r/2 to maximise the effects. The two quantities 
indeed have opposite signs, but their sum only partly cancels, leaving a U-spin breaking 
remainder of —3.0 • 10~ 7 (for 7 = 58° we get —2.6 • 10~ 7 ). This effect is almost entirely 
due to the difference (Fx* — F p ). For form-factor values different from those in table |5.1| 
the U-spin breaking effect would approximately be rescaled proportional to (Fk* — F p ). 
For our choice the U-spin breaking remainder is of the same order of magnitude as the 
individual asymmetries. This example quantifies the limitations of the relation ( [4. 67 ) as 
a Standard Model test. 



5.2 Related Calculations 
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5.2 Related Calculations 

Finally, we want to comment on related work done for B — > Vj. Beneke, Feldmann, 
and Seidel treat exclusive B — > Vl + l~ and B — > decays also in QCD factorization 
||118|| . They use the CMM operator basis ( |2.58| - pl65D supplemented with the semileptonic 
operators Pg and Pio 

P 9 = ^{U) v (sb)v- A (5.7) 

OL - 

Pro = a {8b) V - A (5.8) 

87T 



Another technical difference is that they use the formalism of [ggj with our physical form 
factor Fk* replaced by its HQET analog ^ . In this form factor the hard corrections 
from diagrams as in Fig. [3.4| are taken explicitely into account. The analytic result of 



118|| for the B — > K*^y branching ratio perfectly agrees with ours. The numerical central 
value is 8% larger because a somewhat larger value for the b quark mass was used. The 
interpretation is slightly different as well. Although within the error bars still consistent 
with the experimental averages, Beneke, Feldmann, and Seidel get a little worried by the 
fact that the prediction for the central value of the branching ratio is nevertheless nearly 
twice as large as the experimental central value. So they consider the possibility that the 
form factors at q 2 = are substantially different from what is obtained with QCD sum 
rules or in quark models. Fitting their result for the branching ratio to the experimental 
value results in Fx* = 0.27 ± 0.04, which is almost 30% smaller than the value we gave 



in table FTT 



The paper of Kagan and Neubert ||121|| was mentioned already several times. Herein 
the leading isospin-breaking contributions to the B — > K*^ decay amplitudes are cal- 
culated. In particular a large effect from Qq is obtained. This brings their theoretical 
prediction in good agreement with the current experimental value. Yet, we have to keep 
in mind that there are still large error bars. The definition of the isospin-violating ratio 

a kn = r(P° -> jr°7) - t(b- -> 

0_ r(P° -> K*°~f) + T{B- -> K*~i) 1 ' ' 



Kagan and Neubert use differs by approximately a factor —1/2 from our A5J in ( 4.63 ). 
In addition to the isospin-sensitive annihilation diagrams in Fig. [4.5| they also calcu- 
lated those penguin-annihilation and Q$ hard-spectator contributions of Fig. WM where 



the photon is emitted from one of the upper light quark lines. As explained in section 
|4.6| we included hereof only the annihilation contributions because they are numerically 
enhanced, calculable, and not additionally suppressed by a s (0(mb))- With this choice 
we get A™ = (3.8~^2l)% for our input parameters. The main difference to the Kagan- 
Neubert value of (8.0l|; 2 )% comes from other values of the form factor and B meson 
decay constant used and from the different contributions included. 

In ||116yil7]| the basic mechanisms for the B — > Vj amplitudes at next-to-leading 



order were already discussed. However, a simple model was used to take into account the 
bound-state effects. Additionally, some heuristic methods had to be used to regulate the 
infrared divergences. These emerge from the quark self-energy diagrams and the b — ► 57 
vertex correction and are part of the form factor in our QCD factorization approach. 
The use of hadronic models does not allow a clear separation of short- and long-distance 
dynamics or a clean distinction of model-dependent and model-independent features. 
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Another related work is that of Ali and Parkhomenko where also B — > V'y is considered 



119|| . The first version of this work differed in some points from our results, but there is 
full agreement after revision of the paper. 

Last but not least we comment on a perturbative QCD study of the B — > K*^ decay 



by Li and Lin [|131|| . They calculate most of the relevant diagrams, in particular the 



spectator contribution from Q ± . Yet, in our language, they include type-II corrections 
only, but thereof also power-suppressed contributions. The numerically much smaller 
result for the branching ratio is used to determine the B meson wave function from a fit 
to the experimental data. 



Part III 
The Radiative Decays 



Chapter 6 

The Basic Formulas for B — > 77 



In the third part we discuss the double radiative 5-meson decays B s — > 77 and 1?^ — > 
77 in QCD factorization. We systematically discuss the various contributions to these 
exclusive processes in the heavy-quark limit m& ^> Aqcd- The factorization formula for 
the hadronic matrix elements of local operators in the weak Hamiltonian again allows 
us to separate perturbatively calculable hard scattering kernels from the nonperturbative 
S-meson light-cone distribution amplitude. With power counting in h.QCD/ m b we identify 
leading and subleading contributions to B — > 77. Only one diagram contributes at leading 
power, but an important class of subleading contributions can also be calculated. The 
inclusion of these corrections is used to estimate CP asymmetries in B — > 77 and one- 
particle irreducible two-photon emission from light-quark loops in B and D — > 77. These 
corrections represent the quark-level analogue of so-called long-distance contributions to 
these decays. In this chapter we give the basic formulas to leading logarithmic accuracy 
and discuss the leading-power and subleading-power contributions separately. Numerical 
results are presented and discussed in chapter [7[ Compared to our work in ||132|| we again 
include the QCD penguin operator contributions. 



6.1 Motivation 

For the double radiative decays B — > 77 experimentally so far only upper limits on the 
branching fractions exist: 

B(B° 8 -> 77) < 1.48 • 10" 4 at 90% C.L. [133] 

5(50^77) < 1.7- 10- 6 at 90% C.L. p| [b } 

The Standard Model expectations are roughly two orders of magnitude below these upper 
limits. B — > 77 decays have a rather clean experimental signature but are much better 
suited to be searched for at e + e~ B factories because at hadron colliders the combinatorial 
background is probably too large ||135|| . Furthermore these decays are of interest because 
they could provide useful tests of QCD dynamics in B decays. The B — > 77 modes 
realize the exceptional situation of nontrivial QCD dynamics related to the decaying B in 
conjunction with a completely nonhadronic final state and simple two-body kinematics. 
In principle they also probe the CKM parameters V ts and V t d and could allow us to study 
CP violating effects as the two-photon system can be in a CP-even or CP-odd state. 
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The theoretical treatments of _B S — > 77 performed so far ||136| - |146| all had to employ 
hadronic models to describe the B s meson bound state dynamics. A clear separation 
of short- and long-distance dynamics and a distinction of model-dependent and model- 
independent features were therefore not possible. This concerns especially the dynamics 
of the light-quark constituent inside the B, but also contributions from intermediate J/ip, 
r] c , 0, or D s meson states, which have been discussed as sources of potentially important 
long-distance effects. 



6.2 Basic Formulas 



The effective Hamiltonian for 
addition to the six four-quark 
are ten more gauge-invariant 
numbers to contribute to a b 
combinations of the covariant 
gluonic field strength tensors, 
three photon fields, Grinstein 

H 



b — > 577 is identical to the one for b — > 57 transitions. In 
operators Q1...6 in ( ^.5(J| - |233| ) or Pi. ..6 in ( |2.58 - |2"T6"3 ) there 
dimension-six operators with the right flavour quantum 
— > s transition. These are two-quark operators containing 
derivative = + ieQfA^ + igT a A a ^ the photonic and 
and their duals. Although these operators contain up to 
Springer and Wise showed that the equations of motion 



m q q 



and 



eQ q qYQ 



(6.2) 



can be used to reduce them to the four-quark operator basis and the magnetic penguin 
operators Q7 and Qs |46]. Up to corrections of order l/M^ from higher-dimensional 
operators the effective Hamiltonian for b — > 577 thus is just the one of ( |2.48D for b — > 
sj. We use the standard operator basis ( p.50| - |237|) . The most important operators 



are the magnetic penguin operator Q 7 and the four-quark operators Q\ 2 - The effective 



Hamiltonian for b — > d'yy is obtained from ( |2.48| - p37| ) by the replacement s 
The amplitude for the B — > 77 decay has the general structure 



d. 



A{B -7(fci,ei)7(fc2,e 2 )) = ^^f B \{ll\A + F, v F^ - iA_F^|0> 



Gf Ci 



f B [A + (2ki ■e 2 k 2 -e 1 - m%ei ■ e 2 ) - 2iA-e(kx, h,ei, e 2 



(6.3) 



Here F^ and F^ v are the photon field strength tensor and its dual, where 



2 " 

with e 0123 = — 1 again. The decay rate is then given by 

G F m%f 2 B a 2 



T(B -> 77) 



288vr 3 



(l^ + | 2 + l^-| 2 ) 



(6.4) 



(6.5) 



In the heavy-quark limit we propose a factorization formula for the hadronic matrix 
elements of the operators in the effective Hamiltonian (|2.48|): 



( 7 (ci)7(e 2 )|0i|S> 



(6.6) 



6.3 B 



77 at leading power 



77 



where the q are the polarization 4- vectors of the photons and <3>#i is the leading twist 
light-cone distribution amplitude of the B meson defined in ( |3.19|) . If the light-like vector 
n is chosen appropriately parallel to one of the 4-momenta of the photons only the first 
negative moment of $bi(£) appears, which we parametrize again via Xb = O(Aqcd) a s 



m 



Because there are no hadrons in the final state only one type of hard-scattering kernel 
T (type II or hard-spectator contribution) enters the factorization formula. The QCD 
factorization formula ( |6.6| ) holds up to corrections of relative order AQC£>/m&. The form 
of (|6.6|) with a simple convolution over the light-cone variable £ is appropriate for the 
lowest (leading logarithmic) order in a s , which we will use in the present analysis. A 
generalization to include transverse-momentum variables is likely to be necessary at higher 
orders in QCD Pj88| Jl0l|, [l47 . 



We conclude this section with a brief discussion of CP violation in B — > 77. The 
subscripts ± on A±, defined in (|6.3| ) for B — > 77, denote the CP properties of the corre- 
sponding two-photon final states, which are eigenstates of CP: A C p=+i is proportional to 
the parallel spin polarization t\ ■ e 2 and Acp=-\ is proportional to the perpendicular spin 
polarization t\ x e 2 of the photons. In addition to A± we introduce the CP conjugated 
amplitudes A± for the decay B — > 77 (decaying b anti-quark). Then the deviation of the 
ratios 

I A , I 2 - I A , I 2 

from zero is a measure of direct CP violation. We shall focus here on this direct effect, 
which is specific for the B — > 77 decay, and will not consider CP violation originating in 
B-B mixing. 

Due to the unitarity of the CKM matrix we can parametrize 

A± = \ u ate i0 £ + \ c afe ia ° (6.8) 
A± = Ka±e^ + Kafe^ (6.9) 



where a^ c are real hadronic matrix elements of weak transition operators and the a^ c are 
their CP-conserving phases. We then have 

r ± = 2 Im(\ u \* c )a±at sin(a^ - at) , g 1Q , 

|A u | 2 a^ 2 + |A c | 2 a^ 2 + 2 Re(A u A*)aJa^ cos(a^ — a J) 

The weak phase differences read 

Im(A^A^) = T |A^||AW|sin 7 (6.11) 

with the negative sign for q = s and the positive sign for q = d. The strong phase 
differences sin(o^ — a^) arise from final state interaction (FSI) effects, generated for 
instance via the Bander-Silverman-Soni (BSS) mechanism ||148|| . 



6.3 B —¥ 77 at leading power 

To leading power in AQCD/^b and in the leading logarithmic approximation of QCD only 
one diagram contributes to the amplitude for B — > 77. It is the one-particle reducible 



78 



6. The Basic Formulas for B — > 77 




b > 



s/d 




Figure 6.1: The leading power contribution to the B — > 77 amplitude given by the 
magnetic penguin operator Qj. The diagram with interchanged photons is not shown. 



Radiation from the 6-quark line (Fig. |6.2|) is power suppressed 



(1PR) diagram with the electromagnetic penguin operator Qj, where the second photon 



is emitted from the s-quark line. An illustration is given in Fig. |6.1| . Evaluating the 
factorization formula ( |6.6| ) for this diagram leads to the matrix element 

(6.12) 



( 1 (h,e l )^(k 2 ,e 2 )\Q r \B) 



i^-Q s I d£ _ gljj (2ki ■ e 2 k 2 ■ e 1 - m 2 B e x -e 2 -2i e(k 1} k 2 , e h e 2 )) 



7T 







Together with the effective Hamiltonian fl2.48|) and the definition of A± in (|6.3| ) this 
implies 

At = Aj 9) C 7 ^ (6.13) 



B 



—\[ q \ has been used. 



where CKM unitarity, \ { u q) + \i q) 

After summing over the photon polarizations the branching fraction to leading power 
becomes 



B{B q - 77) 



- Z_!gj ix(g)i2 



r 2 

° 7 3A 2 



(6.14) 



Bo 



where q = d,s for the decay of a B® or a B® meson, respectively. 

In the present approximation the strong-interaction matrix elements multiplying Xu^ 
and Ac 9 '' are identical and have no relative phase. The CP-violating quantities r^ p defined 
in (|6.7|) therefore vanish for the strictly leading-power result. 



6.4 Power-suppressed contributions to B — > 77 

Subleading contributions come from the 1PR diagram, where the second photon is emitted 
from the b quark line (Fig. |6.2j ), and from the one-particle irreducible diagram (1PI) (see 
Fig. U. 

To estimate the subleading 1PR contribution, we simply evaluate the graph in Fig. |6l2] 
using the B- meson projector in ( p.!9| ). The result equals Xb/tub times the leading-power 



expression in Q6.12 ), clearly showing the power suppression of this mechanism. Other 
corrections of the same order can arise from higher-twist terms in the i?-meson wave 
function. We shall neglect all those subleading terms for the matrix element of Qi, 



6.4 Power-suppressed contributions to B — > 77 
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Figure 6.2: The subleading power 1PR diagram of the magnetic penguin operator Qj. 
The diagram with interchanged photons is not shown. 




Figure 6.3: The subleading power 1PI diagram. Here the four-quark operators Qi can 
contribute. The diagram with interchanged photons is not shown. 

keeping in mind that they could naturally contribute relative corrections of order ~ 10%. 
We stress, however, that subleading effects in (Q 7 ) contribute equally to the up- and 
charm-quark components of the amplitude (see ( |2.48| )). Therefore they do not give rise, 
in particular, to relative FSI phases between these sectors, which would affect direct CP 
violation. The same is true for perturbative QCD corrections to (Q7). 



We next consider more closely the diagrams in Fig. 5.3. These 1PI contributions 



which come from the matrix elements of four-quark operators in (|2.48 ), are of special 



interest for two reasons. First, they are the basic effects responsible for a difference be- 
tween the up- and charm-quark sectors of the amplitude, including rescattering phases. 
Second, they represent the parton-level processes that are dual to B — > 77 amplitudes 
from Qi with hadronic intermediate states (J/*&, D^*\ etc.), which are commonly consid- 
ered as generic long-distance contributions. Investigating the amplitude in Fig. ^T3] could 
therefore shed light on this class of effects in B — > 77. 

It shall now be argued that the 1PI contributions are calculable using QCD factoriza- 
tion. To this end we show at 0(a s ) that there are no collinear or soft infrared divergences 
at leading power in l/m^. By leading power we mean here the lowest nonvanishing order 
in the power expansion, where the entire 1PI contribution starts only at subleading power 
with respect to the dominant mechanism in Fig. |6.1| . This suppression will be apparent 
from the explicit expressions given in ( |6.18| - |5T2"B| ) below. The relevant diagrams at 0(a s ) 
are shown in Fig. |6.4j . Following the approach explained in the second reference of [|15| 
for the case of B — > Dn decays, we examine the potentially IR-singular regions of the 
momentum-space integrals in Fig. p.4| . We demonstrate the absence of divergences or 
their mutual cancellation among different diagrams for the hard-scattering part of the 
amplitude at 0(a s ). As a consequence the matrix elements of four-quark operators can 
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s kj 




(a) (b) (c) 




Figure 6.4: The gluon corrections to the 1PI diagram. Quark labels are as in Fig. |6.3| 
and the momentum assignment is everywhere as in diagram (a). The diagrams with 
interchanged photons are not shown. 



be written in factorized form as in (|6.6|) , where in the present case the kernel Tf v is 
independent of £: 

(7(ei)7(e 2 )|gi... 6 |S> = / fl r^ 6 e lM e 2l/ (6.15) 



Power counting shows that all diagrams in Fig. |6.4j are IR finite in the soft-soft and 
hard-soft region as well as when the gluon is collinear with one of the photon momenta 
and the quark in the loop is soft or supersoft. But for a soft or supersoft gluon emit- 
ted from the B meson spectator quark in addition to a collinear quark in the loop we 
encounter logarithmic IR-divergences. To leading power, however, setting the loop mo- 
mentum proportional to the photon momentum, the diagrams vanish individually. The 
last region of phase space we have to consider is where both the gluon momentum s and 
the momentum of the quark in the loop t are collinear with one of the photon momenta, 
for example 

s = otk 2 + s± t = /3k 2 + t± (6.16) 

where a, (3 ~ 1 and sj_, t± ~ Aqcd <C m^. This again leads to a logarithmic IR divergence. 
The leading t = (3k 2 contribution vanishes again as in the (super)soft-collinear case. 
When the gluon is attached to the light spectator quark of the B meson and one of 



the two collinear quarks in the loop (Fig. 3/4 d,e), however, this contribution is power 



enhanced by m b / Kq C d compared to the hard diagram without any gluon (Fig. |673| ). But 
for factorization to hold, we need that the logarithmic divergences have coefficients that 
are power suppressed. We therefore consider the next order in KQco/^b by keeping the 
transverse momentum of t and analyze these terms in more detail. It turns out that the 
contributions from the loop where the gluon hits one of the collinear quarks are both 
proportional to k 2 ■ The factor k 2 is contracted with the open index of the gluon emission 
from one of the B meson constituent quarks. Substituting s = ak 2 and using pt, and / for 
the momenta of the b and the spectator quark in the B meson, respectively, we obtain 
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for this term 



Oih + m b)la . 7a(- y+aft 2 )T 



(p b - ak 2 ) 2 - ml 



+ 



= s 
= 



- af 2 + m h 



a(p b - ak 2 ) 2 - ml 



-l + ak 2 ) 2 

: b _ mb ) ^/ +c ^ 2 )r 



(6.17) 



Here we employed the equations of motion for the constituent quarks of the B meson. 
It follows that the collinear-collinear divergences are also of sub-leading power. This 
completes our proof that at 0(a s ) for the 1PI contribution there are no collinear or soft 
infrared divergences at leading power in l/m^ in the hard-scattering kernel. It supports 
our claim that these diagrams are calculable using QCD factorization. As no assumption 
on the Dirac structure V of the inserted operator had to be made, this demonstration 
applies to all four-quark operators Q1...6 i n an Y possible Fierz ordering. 

For the evaluation of the 1PI diagram in Fig. |6^ we can use the building blocks 
B { Jf A)x{VTA) and £g+p)x(s-P) of flOoD and flQ3|) . We simply have to substitute g s T a -> 
Q q e, q — > ki, r — > k 2 , and can use the on-shell relations k 2 = 0, q • ki = for both 
photons. Having regard to the possible insertions of the Fierz-transformed QCD penguin 
operators, we find herewith the matrix elements 



(Qi) 


fnoc n 


[z p )e(k 1 ,k 2 ,e 1 ,e 2 ) 






(6.18) 


(Ql) 


= N(Ql) 










(6.19) 


(Qs) 


= —£{ki 

71 


k 2 ,e 1 ,e 2 ) 


q 


z q ) + Q 2 d (g(0) + g(l)) 


(6.20) 




= —£{ki 

71 


k 2 ,e 1 ,e 2 ) 


q 


) + NQ 2 d (g(Q) + 9(X)) 


(6.21) 


(Qb) 


fB q a{ 

7T V ' 


&2,ei,e 2 ) 




z q ) + Q 2 b (g(l) + 2) 




(6.22) 



-iQl 



5,(1) 



2fci • e 2 k 2 ■ e x - m\ e x ■ e 2 



(Qe) 



fB„a 



71 



e(k 1 ,k 2 ,e 1 ,e 2 ) 



Y,Q 2 Azq)+NQl{g{l) + 2) 



(6.23) 



-iNQl 



1-- 2 9(D 



2ki ■ e 2 k 2 ■ e x - m 2 B e\ ■ e 2 



with 



9W 



-2 + Az 



1 - VI - 4z + ie 



+ L 2 



1 + Vl-4z + 



(6.24) 



Here the dilogarithmic function was defined already in ( 4.33|) , z q is again m 2 /ml, and 
p = u,c. We note that the function g(z) is related to the kernel h(u, z) defined in ( |4.32| ) 
in the context of B — > by g(z) = h(l,z). The function g(z) is regular for z — > 0, 
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c/(0) = —2, and has an imaginary part growing almost linearly with z to sizeable values 
for z = z c and vanishing again for z > 1/4. The imaginary part is due to the rescattering 
process B — > qq — > 77 with an on-shell light-quark pair qq in the intermediate state. Since 
the process B — > qq is forbidden by helicity and angular momentum conservation for a 
massless quark q, the function g(z) is real at z = (that is for q = u,d, s). The helicity 
suppression of the phase will be absent at higher order in a s . 

Including the effect of 1PI diagrams in the decay amplitudes, the quantities A± from 
( |6.13| ) become 



A\ 
A v _ 



-C 7 ^- + (C 5 + NC 6 ) 



-C/-^- - \ (d + NC 2 ) g(z„) - (C 3 - C 5 ) 



- (C 4 - C 6 ) 



2^ c ) + -^(1; 



2() ^ 16^ 

H C 3 + 4C 4 C 5 

3 3 



(6.25) 



(6.26) 



where we have decomposed A± into 



Aj 



E 



p=«,c 



(6.27) 



neglected the light quark masses and used g(0) = —2. We see that the dominant effect 
is a modification of A_ via the current- current operators Q\ and Qi- The imaginary 
part of the 1PI loop diagram with a charm quark running in the loop leads to a relative 
CP-conserving phase between the X u and the A c contribution. This gives a nonvanishing 
r^p while r^ p is still zero because A" = A c + . The effect of QCD penguin operators is 
small because of their tiny Wilson coefficients. 



6.5 Long-distance contributions to B and D 



11 



Long-distance contributions were calculated for both B — > 77 ||141| , |142| , |149| , |150|| and D 



77 |[L51]J152 |. Different mechanisms were considered: 



B and D — > Vj followed by a V —>■ 7 conversion with the conversion factor supplied 
by the vector meson dominance (VMD) model. This gives long-distance contri- 
butions to the operator Q 7 . Counting powers of in the explicit result for the 
B s — > 07 — ► 77 amplitude in [|141|| shows that it is suppressed with Agc^/mb 
compared to the leading power result from Fig. |6.1| . Yet, it is just one among the 
subleading power terms for the matrix element of Q 7 . Another comes for example 
from the diagram in Fig. |5T2". 



We can apply a similar power counting to the D — > 77 VMD amplitude of ||152 |. 
In the limit of a heavy charm quark it is likewise power suppressed with Aqcd/^c 
compared to the short- distance amplitude. The limit m c 3> Aqcd is certainly ques- 
tionable, but nevertheless some properties of charmed hadrons have been analyzed 
successfully in this limit ||153||. Again, the VMD contributions are just one out of 



several possible power corrections to the Q 7 matrix element. 
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Single-particle unitarity contributions where the D° mixes with a spinless interme- 
diate particle, which decays into a photon pair. Such a process corresponds to our 
1PI contributions in Fig. |6.3| . Relative to this contribution the amplitudes in [152 
are at most of the same order in powers of the inverse "heavy" charm quark mass. 

Two-particle unitarity contributions are other hadron level counterparts of our 1PI 
diagrams. Triangle loops with intermediate Dg ||149| , 150 l and [152] mesons 
were calculated for B s — > 77 and D — > 77, respectively. As we will see in the next 
chapter, the effect of the 1PI contributions in B — > 77 on the branching ratio is very 
small. The authors of [|150|| also get a not too large effect, reestimating the result 
of [|149|| . Yet, one has to keep in mind that the effects from intermediate D and D* 
mesons are only two of many possible hadronic intermediate states. 

With obvious replacements we can use our 1PI results in (|6.18| - |6~23l) to estimate 
the quark-level analogue of the two-particle unitarity contributions to D —>■ 77 
calculated in [152]. For the current-current operators we get for the 1PI amplitude 
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(6.28) 



where we set fr> = 200 MeV and evaluated the Wilson coefficients at m c . The 
contributions of a strange quark and a down quark in the loop are accompanied by 
CKM elements which are practically equal in absolute value but opposite in sign. 
This leads to a strong GIM cancellation in the square bracket. If we neglected the 
strange quark mass, the 1PI contribution would vanish altogether. This is a crucial 
difference of the c —>■ w~y transition at one loop compared to the b — > 57 transition 
and disturbs the hierarchy of higher order QCD terms. Leading logarithmic QCD 
corrections are known to enhance the amplitude by more than an order of magnitude. 
Including the two-loop QCD contributions increases the amplitude further by two 
orders of magnitude ||154|| . Therefore we would like to compare equation ( |6.28| ) with 
the c — > wy amplitude from the dominant two-loop diagrams \A\ = 0.0047 of ||154 . 



Evaluating the factorization formula then leads to the following "short distance" 
amplitude: 



^SD 



I ^£)->77| 
I^SD - I 



Gf ol 



_ -f D Q u A^- [2ie(ei, e 2 , h, k 2 ) - 2ki-e 2 k 2 -e x + m 2 D e v e 2 ] 
V 2 7T Ad 
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-1, 



(6.29) 



We see that the 1PI amplitude is smaller than the CP-odd part of the short distance 
amplitude. This should be understood as an estimate only because we relied on the 
heavy quark limit for the charm quark. 

Burdman, Golowich, Hewett, and Pakvasa on the other hand find a large effect on 
the D — > 77 branching ratio from K + K~ intermediate states ||152|| . Their amplitude, 
however, is formally suppressed with three powers relative to our 1PI amplitude. In 
view of our above estimates the large effect could be overrated. Only one interme- 
diate state was taken into account so that important contributions in the duality 
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sum might be missing. We also expect cancellations due to the GIM mechanism as 
explicitely seen in our estimate above. 

The chain process B s — > <fnf> — > 07 — > 77 was estimated to give very small effects 
on the B — >• 77 decay amplitudes [|142|| . It is multiply power suppressed compared 



to the leading power 1PR contribution. 



Chapter 7 

The Phenomenology of B —> 77 



Here we present numerical results for the expressions derived in the previous chapter. 



Additionally to the input parameters given in table |5.1| we need the parameters related 
to the B s meson. These we summarize in table |7.1| . We will work to leading logarithmic 
order in QCD and use the running mass m 4 (m 4 ) = 166 GeV for the top quark as input for 
the Wilson coefficient C-j. For defmiteness we employ the two-loop form of the running 
coupling a s (fi) as quoted in fl2.19p , which corresponds to a s (Mz) = 0.118 for = 
225 MeV. With central values of all input parameters, at /i = m&, and using a nominal 
value of the CKM angle 7 = 58° we find for the branching ratios to leading logarithmic 
accuracy, evaluating (|6.5|), ( |6.25| - |6l2^ ): 



77) = 1.22 -10- 6 (7.1) 
B(B d ^j-y) = 3.10 -10- 8 (7.2) 

Using the leading power amplitudes (|6.13|) without the IPI contributions we get 



Adg. P ow.(^^77) = 1.17 -10" 6 (7.3) 
5i dg . pow .(5 d ^77) = 3.12 -10- 8 (7.4) 

To display the size of power corrections from IPI diagrams on amplitude level we compare 
A p + and A v _ in ( |6.25[ ) and ( |6.26| ) with the leading power result ^4± ldgpow = —C 7 m B J \ Bs 
for B s -» 77 

A u/c = 4g72 +0.045 = 4.917 

A u _ = 4.872 +0.419 -0.061 - 0.033 i = 5.230 - 0.033 i 

At = 4.872 +0.203 -0.534i -0.061 - 0.033i = 5.014 - 0.567i ^ 7 - 5 ) 

ldg. pow. Qi t 2 contribution Q3...6 contribution 



Parameters related to the B s meson 










5.37 GeV 


(230 ± 30) MeV 


1.49 ps 


(350 ± 150) MeV 



Table 7.1: The input parameters for the B s meson. 
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where we showed separately the leading power Q-j contribution and the effects from 
current- current and QCD-penguin operators to the 1PI contribution. For B d — > 77 nomi- 
nally only A p ± ldgpow = —C 7 m B J '\ Ba = 4.79 is changed. The calculable power corrections 
thus are of the expected canonical size of 0(10%). The precise value of these power- 
suppressed effects depends strongly on the renormalization scale \i. Typically, the relative 
size of the 1PI terms diminishes when the scale is lowered, both due to an increase in C7 
as well as the simultaneous decrease of the combination a^- The effect of QCD penguin 
operators to the 1PI diagrams is small. 

The main errors in ( |7. l| - [7 5|) come from the variation of the nonperturbative input 
parameters Xb and the decay constants fs d and fs s which are all poorly known and 
enter the branching ratios quadratically. The residual scale dependence is sizeable as 
well because we calculated at leading logarithmic accuracy only. The variation of the 
branching ratio with the scale would be less severe if next-to-leading QCD corrections 
were known. We expect the NLL corrections to increase the branching ratio as was the 
case for both the inclusive |nj and exclusive | |110| , |118| | b — > 57 decays. 



The effect of next-to-leading logarithmic corrections for the branching ratio of b — > S7 
can be reproduced by choosing a low renormalization scale (fi ~ rrib/2) in the leading 
logarithmic expressions. If such a low scale were also relevant for B — > 77 the branching 
ratios would read 

B{B S 77W/2 = 1-52 -HT 6 (7.6) 
5(^^77)^/2 = 4.06 -HT 8 (7.7) 

However, we prefer to use the nominal \x = and to quote the standard scale ambiguity 
with rrib/2 < u. < 2m b theoretical uncertainty. 

For comparison we also show the results for the case where short-distance QCD effects 
are neglected altogether. This amounts to taking fi = M w in the leading-order Wilson 
coefficients and leads to 

5(S s ^7 T ) M=Mw = 0.60 -10- 6 (7.8) 
B(B d -+ 77)m=Mw = 1-22 -10- 8 (7.9) 

Clearly, the leading logarithmic QCD corrections lead to a substantial enhancement of 
the branching ratios ([71|), ( |7.2| ), similar to the case of b — > S7. 

The numerical evaluation of fl6.10| ) gives for central values of the input parameters and 
at /i — rrih the following predictions for the CP- violating ratios r^ p : 

rcp(B s ) = 0.35% (7.10) 
rcp(B d ) = -9.66% (7.11) 

The effect is negligible for B s — ► 77 as expected. For B d — ► 77 CP violation occurs at 
the level of about 10%. It is interesting to note that in B — ► 77 a non-vanishing CP 
asymmetry appears already at 0(oP s ). This is possible via a modified BSS mechanism: 
Instead of the usual QCD penguin, already the electroweak loop in Fig. [U^ itself leads to a 
CP-conserving phase. This phase difference defined in ( p.lOQ is nonvanishing for the CP- 
odd amplitude only: a~ — a~ = —6.1°. Including the effects of QCD penguin operators 
does not lead to a CP-conserving phase difference a+ — and therefore r^ p = 0. 
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Figure 7.1: The CP- violating ratios r^ p for decays of neutral and B s mesons to two 
photons as a function of the CKM angle 7, each for three values of the renormalization 
scale /1 = rrib/2, and 2rrib- 



The scale dependence of the CP asymmetry is rather strong, as can be seen in Fig. |7Tl . 
The ratio rep also depends sensitively on fundamental CKM parameters such as the CKM 
angle 7 (Fig. |7.1|) . The extremal value of r^ P (Bd) is obtained for 7 = 50°. The sensitivity 
of the branching ratios and CP asymmetries to variations in the relevant input parameters 



is summarized in table 7.2 



Our predictions for the central values of the branching ratios are roughly two orders of 
magnitude below the current experimental bounds. The recent upper limit for B{B<i — > 
77) from BaBar in (ETT) improved the previous limit from the L3 collaboration ||155|| 
by a factor of 20. With more integrated luminosity accumulated the upper bound will 
further be reduced. Yet, even with the planned total integrated luminosity of 500 fb -1 
for each BaBar and BELLE, the number of — > 77 decays one expects to observe is 
only a single one. The high-luminosity option SuperBABAR proposes to collect a total 
integrated luminosity of 10ab _1 . Therefore the number of observed events for Bd — > 77 
would be two dozens [|156|| . As long as the e + e~ B factories run at the T(4s) resonance 
they produce no B s mesons and therefore cannot profit from the larger branching fraction 
to observe B s — > 77 decays. At hadron machines as the Fermilab Tevatron or the LHC 
at CERN all sorts of 6-flavoured hadrons are produced. CDF, for instance, expects to 
observe at the present run Ha at Tevatron in 2fb _1 of data 28 million b — > J/ipX decays 
The number of B s mesons will therefore be roughly 7 million. It will consequently 
be very hard to distinguish B s — > 77 events from the immense combinatorial background 
events. 

As already mentioned, the dominant uncertainty for our prediction of the branching 
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B(B d -> 77) 


5(5, - 77) 


r cp( B d) 


r cp( B s) 




[10- 8 ] 


[10- 6 ] 


[%] 


[%] 


central 


3.09 


1.22 


-9.66 


0.35 


Ab 


+6.40/-1.58 


2.44/-0.61 


+4.13/-4.02 


+0.15/-0.15 


/b 


+1.00/-0.86 


+0.34/-0.30 






/i 


+0.97/-0.66 


+0.31/-0.21 


+9.03/-8.13 


+0.30/-0.33 


m c 


+0.10/-0.12 


+0.03/-0.04 


+1.79/-1.58 


+0.06/-0.06 


7 = (58 ± 24)° 


+1.28/-0.95 


+0.01/-0.02 


+1.66/-0.25 


+0.06/-0.12 



Table 7.2: Predictions for branching ratios and CP asymmetries with the errors from the 
individual input uncertainties. 

ratios comes from the variation of the hadronic parameter A# (factor 2). It governs the 



uncertainty of the leading power amplitude. So far, our value for Ab quoted in table [54] is 
only an educated guess with the uncertainty taken appropriately large. As enters the 
leading power branching ratio quadratically a large uncertainty in this input parameter is 
translated into an even larger uncertainty of the branching ratio output. It would be highly 
desirable to get a better understanding of A#. We hope that this problem can be attacked 
soon with the help of QCD sum rule or lattice QCD techniques. As As parametrizes the 
first negative moment of the B meson wave function it is a universal quantity and therefore 
independent of the precise decay under consideration. An extraction of it therefore might 
be possible from radiative semileptonic B decays, for instance. 

Leaving the uncertainty of A^ aside, the variation of the remaining parameters changes 
the B d — > 77 and B s — > 77 branching ratios by order ±50% and ±35%, respectively. If 
both the Bd — > 77 and B s — > 77 branching ratios were measured one could consider their 
ratio where the dominant uncertainties from the Xb and Jb cancel to a large extent. 

The main theoretical uncertainty for the CP asymmetries at present comes from the 
variation of the renormalization scale \i. This could be reduced when QCD corrections 
are included. But as the CP asymmetries in double radiative exclusive decays are ex- 
tremely hard to measure experimentally we see at present no urgent need to perform this 
calculation. 



Part IV 
Conclusions 



Chapter 8 
Conclusions 



While waiting for the completion of LHC and new results from Tevatron Run II, B physics 
is among the most active fields in recent particle physics. The main goal of B physics is 
a precision study of the flavour sector to extract Standard Model parameters with high 
precision in order to reveal potential New Physics effects. Theoretically, the complication 
is that three fundamental scales are involved in B decays: the weak interaction scale Mw, 
the b quark mass mj, and the QCD scale Aqcd- Via the machinery of operator product 
expansion and renormalization group equations in the framework of an effective theory 
one can factorize perturbatively calculable short distance Wilson coefficients from the long 
distance operator matrix elements. The latter pose the main difficulty as all low-energy 
contributions below the factorization scale fi = 0{m^) are collected into them. Yet, one 
can take advantage of the fact that ^> Aqcd and achieve a further separation of short- 
and long-distance part of these matrix elements. Employed for exclusive final states this 
application of QCD factorization amounts to a systematic expansion in Aqcd/to;, and the 
strong coupling parameter a s . 

Within this thesis we have analyzed in detail the exclusive radiative decays B — > V7 
(V = K* or p) and B — > 77 in QCD factorization. This enabled us to separate per- 
turbatively calculable contributions from the nonperturbative form factors and universal 
meson light-cone distribution amplitudes. A power counting in AQ CD /m b implies a hi- 
erarchy among the possible transition mechanisms and allows to identify leading and 
subleading contributions. 

Our framework has major advantages when compared with previous approaches to 
exclusive radiative decays. They all had to rely on hadronic models to estimate bound 
state effects. These models, however, did not allow a clear separation of short- and 
long-distance dynamics or a clean distinction of model-dependent and model-independent 
features. QCD factorization is a systematic and model-independent approach based on 
first principles and the heavy-quark limit ^> Aqcd- It holds for ^00. In real- 
ity, however, the bottom quark is not infinitely heavy so that we have to worry about 
Aqcd/^6 power corrections. For the radiative B decays we considered, some of these 
power corrections - possibly the numerically most important ones - are calculable in QCD 
factorization. The remaining ones are not expected to be very large. In particular, no 
"chirally enhanced" contributions, as in non-leptonic modes with pseudoscalar mesons in 
the final state, can arise in our case. 

Sizeable uncertainties come from the variation of the nonperturbative input parame- 
ters, in particular form factors, decay constants, and information on the B meson wave 



92 



8. Conclusions 



function (Xb), which are all poorly known. But we await progress on these quantities 
from lattice QCD and analytical methods as well as from an extraction from experimen- 
tal results. 

An important conceptual aspect of our analyses is the interpretation of quark loop 
contributions. We showed that these effects are calculable in QCD factorization. Before, 
such quark loops were considered as generic, uncalculable long-distance contributions. 
Non-factorizable long-distance corrections may still exist, but they are power-suppressed. 
Let us in the following address our findings for B — > V'-f and B — > 77 separately. 

The radiative decay with a vector meson K* in the final state is of particular promi- 
nence as its branching ratio was already measured. We presented in this work the first 
really complete calculation with next-to-leading logarithmic accuracy. Compared to our 
publication [|110|| we also included the effect of QCD penguin operators and added an 
explicit proof, at one-loop order, for the factorizability of annihilation contributions. At 
next-to-leading order in QCD interesting effects occur: the operators Qi..s and Q 8 start 
contributing, quark loops allow for CP-violating effects, and the hard-spectator correc- 
tions appear. We have seen that weak- annihilation amplitudes are power suppressed, but 
calculable in QCD factorization and numerically important because they enter with large 
coefficients. Therefore, we included their effects in our phenomenological discussion. We 
estimated other power corrections to be numerically small. 

Our NLO predictions for the central values of the branching ratios are 

B(B° -> i?*° 7 ) = (7.4^-6) . io-5 B(B~ -> K*~j) = (7.3±1;|) • KT 5 

B{B o = (o.8l8 : |). ID" 6 B{Br-+p~i) = (1.6^) ■ 1(T 6 

They are substantially larger than the leading logarithmic values when the same form 
factors are used. The dominant correction comes from the hard-vertex contributions, 
which similarly increased the result for the inclusive B — > X s j decay. Adding the effect 
of varying the input parameters in quadrature leads to an error estimate of (9(35%) and 
(9(40%) for the B — > K*^ and B — > p7 branching ratios, respectively. For B(B — > K* r f) 
by far the dominant uncertainty comes from the B — > K* form factor, which alone already 
amounts to 30%. This situation, however, can be systematically improved. In particular, 
our approach allows for a consistent perturbative matching of the form factor to the 
short- distance part of the amplitude. The largest uncertainty in B(B — > fry) is also the 
form factor. Yet, variation of the hadronic parameter A 5 and the CKM angle 7 affect the 
prediction for the branching ratio also considerably. 

Compared to the experimental measurements of the B — > K*j branching ratios our 
predicted central values are somewhat high. Although within the uncertainties still com- 
patible, the discrepancy of the central values could be a hint that the form factors at 
q 2 = are smaller than the value obtained from QCD sum rules. At the present stage we 
also cannot exclude that power corrections or even a failure of the QCD factorization ap- 
proach are responsible for the difference. Because the experimentally measured inclusive 
b — > 57 rate is in good agreement with its theoretical prediction we do not expect sizeable 
New Physics contributions to the exclusive B — > K*^ decays either. For the b — *> d'y 
transitions on the other hand there is still room for physics beyond the Standard Model 
as experimentally at present only upper bounds for these decays exist. They are only a 
factor of two above our predictions so that we expect the B — > p7 branching ratios to 
be measured soon. From such a measurement we can obtain important information on 
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our input parameters, in particular the form factor, the hadronic parameter Xb, and the 
CKM angle 7, and also on the validity of the QCD factorization approach to exclusive 
radiative B decays. 

For B —>■ fry both CKM sectors of the effective Hamiltonian have the same order 
of magnitude. This allows a sizeable CP- violating asymmetry of A{p ±r y) = 10.1%. In 
B — > -fT*7, on the other hand, |Al s ^| <C [ Ac | and therefore the CP asymmetry is tiny, 
typically —0.5%. Because the CP asymmetry arises at 0(a s ) for the first time, we have 
a large dependence on the renormalization scale p. This could be reduced by going one 
order further in perturbation theory, which would amount to a three-loop calculation. 

We included the numerically most important power corrections - weak annihilation - 
into our analysis. Within our approximations weak annihilation is the only effect that is 
sensitive to the charge of the B meson and therefore to isospin breaking. In B — > if* 7 
annihilation contributions from are CKM suppressed such that their influence on 
the branching ratio is small but a curiosity is a rather large effect of the QCD penguin 
operator Qq ||121|| . So for both B — > K*^ and B — > p7 the isospin asymmetry is predicted 



to be sizeable: A(K*-f) = (-7.5j££)% and A(p 7 ) = (2.0t?£?)%. The isospin breaking in 
B — > p7 depends very sensitively on the CKM angle 7 and can thus in principle serve to 
constrain this quantity once measurements become available. 

An explicit estimate of U-spin breaking effects in B — > V'y decays quantified the 
limitations of the use of U-spin symmetry in exclusive decays as a Standard Model test. 

Our formalism can also straightforwardly be applied to other radiative rare B decays, 
such as B s — > V'-f where V = <fr, K*, B& — > uj'-f, or radiative decays to higher resonances. 
Some of these modes are discussed in | ]120| , |157| |. Yet, for QCD factorization it is again the 



form factors, which pose the problem. To our knowledge they were never calculated. This 
would be highly desireable because for example the B — > 7^2(1430)7 decay was already 



experimentally observed [p8||158| . 



The double radiative decays B s — > 77 and B& — > 77 were not seen in experiment 
yet. Within QCD factorization we predict the branching ratios to leading logarithmic 
accuracy: 

B(S S ^77) = (1.2j$) • 10- 6 
2?(fl d -77) = (3.1±1;D • 10- 8 

In the heavy quark limit the dominant contribution at leading power comes from only a 
single diagram, but an important class of subleading contributions can also be calculated. 
We used these corrections to estimate CP asymmetries in B — > 77 and one-particle 
irreducible two-photon emission from quark loops in D — ► 77. The CP-violating ratio 
r^ p (B s ) = 0.4% is negligible and r^ p (B d ) = —9.7%, albeit not very small, will not be 
measured in the near future. 

The absolutely dominant uncertainty in the branching ratios comes from the variation 
of the hadronic parameter A#- This universal quantity is theoretically very poorly known 
and enters the branching ratio quadratically. We propose to further investigate Xb which 
is of relevance also for other applications of QCD factorization to exclusive B meson 
decays. 



In table |S.l| we summarize the decay reach of the modes we discussed. The information 
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Hadron Collider 


e + e i? Factories 






Experiments 






Decay 


Branching 


CDF 


BTeV 


BaBar 


Super- 


Mode 


Fractions 


DO 


LHC-b 


Belle 


BaBar 






(21b- 1 ) 


(10 7 s) 


(0.5 ab _1 ) 


(10 ab" 1 ) 


B -> K*-f 


5 ■ 1(T 5 


170 


25K 


6K 


120K 


B -> p7 


2 • 10~ 6 






300 


6K 


- 77 


3- 10~ 8 






1.3 


25 



Table 8.1: Decay reach of B experiments for exclusive radiative decays [|156|| . 



is taken from the Super-BABAR proposal ||156 |. Although in our estimate the branching 
fraction for B® — ► 77 is by a factor of three larger compared to the value in ||156|| , these 
decays are indeed very rare. The B® — > 77 modes are roughly 40 times more frequent, but 
they cannot be produced at the T(4s) resonance. At hadron colliders, on the other hand, 
also B s mesons are produced. But the problem with the pj and 77 final state in a hadron 
collider environment is that the combinatorial background is extremely large. Photons 
and p mesons are produced copiously. Therefore, in table |S.1| no estimates exist how many 
of these decays could be really observed. Nevertheless we trust in the excellence of our 
experimental colleagues that they manage to measure the B — > p7 or even B — > 77 modes 
at hadron collider experiments. Running a prospective high-energy e + e~ linear collider 
on the Z resonance could also improve the bounds on the B — > 77 branching ratios or 
even turn them into measurements. 



According to Niels Bohr, the "goal of science is to gradually reduce our prejudice" 
159| . We think we made one step in reducing the prejudice against a theoretical calcu- 
lation of exclusive decays. Many exclusive B decays can be treated successfully in QCD 
factorization and the radiative ones are part of this set. QCD factorization has put the 
fear of God into hadronic matrix elements of exclusive B decays. Still a lot of work has 
to be done, especially concerning the power corrections. But we have for the first time 
also for exclusive decays a tool which is based on a systematic framework, namely the 
heavy-quark limit. We will see if the phenomenological impact of QCD factorization is 
equal. 



Part V 
Appendices 



Appendix A 

Wilson Coefficients for b —> 57 



For completeness we want to give the explicit formulas for the Wilson coefficients of the 
operator bases ( |2.5Q - p37| ) and ( |2.58| - |2T65| ) in the naive dimensional regularization and 
MS renormalization scheme. 

Let us consider the magnetic penguin diagrams in Fig. |A.1| The loop integral of 



these diagrams contains two 7 matrices. Therefore, insertions with the Dirac structure 
7mi • • • 7/i2n+i (1 — 75) ® 7 Ml • • • 7 At2n+1 (1 — 75) vanish. This is in particular the case for Q1...4 
and is eventually a consequence of chirality conservation in the quark loop and QED 
gauge invariance. Non-vanishing results come from Q5 and Qq only in the NDR scheme. 
The magnetic photon penguin diagrams gives 



(Q 5 ) (a) 
(Q 6 ) (a) 



1 e 
~3 4^ 

N e 

T 4^2 



s i im h {\ + 75)6 
s i 4m h {\ + 75)6 



y5(si\Q 7 \b) 
VeisilQrlb) 



(A.l) 



For the insertion in the magnetic gluon penguin diagram we replace the photon vertex 
with the gluon vertex and get 

■^T a s ^m b (l + l5 )b =: z 5 (s-f\Q 8 \b) 
47r (A.2) 

(Qe) (b) = =: ^ 6 <^|Q 8 |&> 



(Q 5 ) 



(6) 



The above matrix elements vanish for an on-shell photon in any 4-dimensional regulariza- 
tion scheme and also the HV scheme. Such a peculiar regularization scheme dependence of 




(a) 




Figure A.l: Magnetic photon and gluon penguin diagrams. 
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leading order matrix elements has to be cancelled by a corresponding scheme dependence 
of the LO anomalous dimension matrix [47|. This is disturbing because 7^°- ) is usually 



scheme independent. It is convenient to introduce the so-called "effective coefficients" for 
the operators Q7 and Q$ ||160|| 



cf «"(/,) = Cf +5>Cf(A«) (A.3) 

i=l 

5>^°V) ( A - 4 ) 

i=l 

(0, 0, 0, 0, 1, 0) in the NDR scheme. Correspond- 



a 



(0) 



with y = (0, 0, 0, 0, -|, -f ) and z 

ing values for the CMM operator basis ( p.5^ - p^65| ) are y = (0, 0, — ~, — |, — 22, — §2) and 



3 ' 



z — (0, 0, 1, — |, 20, — y). The leading order effective coefficients C\ with the corre- 
sponding anomalous dimension matrix ^^ e ff then are regularization and renormalization 
scheme independent. In addition, the leading order b — > S7 and b — > sg matrix elements 
are proportional to C^ e ^ and C^ e ^ , respectively. In order to simplify the notation we 
will omit the label "eff" throughout this work. Whenever C-j or Cs appear, the effective 
coefficients are understood. 

Let us expand the Cj in powers of a s \ 



1(1), 



(A.5) 



The evolution of the (effective) Wilson coefficients is driven by the (effective) anomalous 
dimension matrix. Its leading order in the expansion ( 2.40[) is given for our operator basis 
( ^50H237|) by 
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3 



\ 



(A.6) 



Here we have already set the number of quark colours to 3. The number of active up- and 
down-type flavours is denoted by u and d, respectively and f = u + d. The corresponding 
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ADM for the Chetyrkin, Misiak, Miinz operator basis is [51 



(0) 
7CMM 
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-4 


8 
3 





2 
9 








208 
243 


173 
162 


\ 




12 








4 
3 








416 
81 


70 
27 















52 
3 





2 


176 
81 


11 

27 












40 
9 


100 

9 


1 

9 


5 
6 


152 
243 


587 
162 















256 
3 





20 


6272 
81 


6596 
27 












256 
9 


56 
9 


10 

9 


2 
3 


4624 
243 


4772 
81 
























32 
3 
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32 
9 


28 
3 


/ 
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These two matrices are related via the leading order expansion of equation ( |B.7| ) 

tSm = W^W- 1 (A.8) 

with W given in ( |B.8|) . At next-to- leading order the anomalous dimension matrices not 
only look complicated but are very difficult to calculate. The complete 8x8 ADM requires 
performing three-loop renormalization of the effective theory. This calculation was done 
in the CMM basis with the result 



7cMM 



/ 


355 


502 


1412 


1369 


134 


35 


818 


3779 


9 


27 


243 


243 


243 


162 


243 


324 




35 


28 


416 


1280 


56 


35 


508 


1841 




3 


3 


81 


81 


81 


27 


81 


108 










4468 


31469 


400 


3373 


22348 


10178 




81 


81 


81 
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243 


81 










8158 
243 


59399 
243 


269 
486 


12899 
648 


17584 
243 


172471 
648 










251680 
81 


128648 
81 


23836 
81 


6106 
27 


1183696 
729 


2901296 
243 










58640 


26348 


14324 


2551 


2480344 


3296257 




243 


243 


243 


162 


2187 


729 






















4688 
27 
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2192 
81 


4063 
27 



(A.9) 



/ 

For the standard operator basis the 6x6 submatrix governing the mixing of current- 
current and QCD-penguin operators among each other is known to be 



7 



(i) 
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21 2/ 
2 9 


1+ 


2/ 
3 


202 
243 


1354 
81 


1192 
243 


904 \ 
81 ' 




7,2/ 
2 ~ 3 


21 
2 


2/ 
9 


79 
9 


7 
3 


65 
9 


7 

3 












5911 , 71/ 
486 ' 9 


5983 , / 
162 ' 3 


2384 71/ 
243 9 


1808 £ 
81 3 












379 , 56/ 
18 ' 243 


91 _|_ 808/ 


130 502/ 
9 243 


14 , 646/ 
3 + 81 












61/ 
9 


11/ 

3 


71 , 61/ 
3 ~ 9 




V 










682/ 
243 


106/ 
81 


225 , 1676/ 
2 243 


1343 , 1348/ , 
6 1 81 ' 



(A. 10) 

More on the transformation relations between the two operator bases can be found in the 
subsequent Appendix O. 
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As explained in section |2.4| , the initial conditions for the Wilson coefficients are found 
by matching the effective theory amplitudes with the full Standard Model ones. At 
Ho = Mw the leading-order results are 



C?\M W ) 



I 




0^3 _ o, r .2 

t \nx t + 



5x? + 7x t 



4(x t - f) 4 
-3x f 2 

4(x t - f 



lnxj + 



24(x t - i) 3 
3 + hx\ + 2x t 



for z 
for z 

for z 



i 

2. 



.6 



(A.ll) 



for z = 8 



where 



.7 ; 



Ml 



(A.12) 



Up and charm quark mass are much smaller than the top quark mass and were therefore 
neglected in ( |A.11| ). The only difference for the initial conditions in the CMM basis is 
that the indices 1 and 2 are interchanged. Equating out (|2.41|) we get for our standard 
operator basis 



and 



kjiif 1 for j = 1 ... 6 

i=l 

^Cf\M w ) + \(rfi- rfi) Cf\M w ) + £ k 7iV »< 

8 



i=l 



i=l 



(A.13) 
(A.14) 
(A.15) 



y~] hji-if 1 for j = 1 ... 6 



8=1 



z<°V) = cf ( M ) 

4V) = Cf\») 



(A.16) 

(A.17) 
(A.18) 



for the CMM basis. 

At next-to- leading order we only need the Wilson coefficient CV. As explained in 



Appendix [FJ we have C\ (a) = Z\ (a) and can directly take the formula from JO 



/ 297664 16 7164416 m 
+ ( 1 ^000 v 23 — nrn ^ r V 23 



256868 37 6698884 39 \ (0) 



14283 



357075 



-?/- 



37208 / 39 16 
H — ( ?y 23 — r) M 



476 i 



14283 ' 357075 

8 

cf\M w ) + ^(e^x) + ft + g iV )v ai 



(A.19) 



i=l 
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The initial conditions for the NLO Wilson coefficients CV and Cg are 

^1}/,, x -16x 4 -122a; 3 + 80a; 2 -8a; / 1\ 6a; 4 + 46a; 3 - 28x 2 , 2 

C^»'> = L < 1 -;j + 3(,-l)» 

-102x 5 - 588x 4 - 2262x 3 + 3244x 2 - 1364a; + 208 



\l(x- 1 



1 5 



hix 



1646x 4 + 12205a; 3 - 10740x 2 + 2509a; - 436 . A rt , 

+ 486(131)^ < A ' 20 > 

, , , , r -4x 4 + 40a; 3 + 41a; 2 + x r A 1 \ -17a; 3 - 31x 2 , 2 



c ^>w = — — L i 1 -*) + 2(,-i)» ln * 

-210a; 5 + 1086a; 4 + 4893a; 3 + 2857x 2 - 1994a; + 280 , 

H ; In a; 

216(ar- l) 5 

737x 4 - 14102a; 3 - 28209a; 2 + 610a; - 508 



1296(x- l) 4 



(A.21) 



All the "magic numbers" for the Wilson coefficients necessary in a NLL analysis of b — > 57 
processes are summarized in table [A.l. 
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i 


1 


2 


3 


4 


5 


6 


7 


8 




14 
23 


16 
23 


6 
23 


12 

23 


0.4086 


-0.4230 


-0.8994 


0.1456 


hi 








1 

2 


1 

2 














hi 








1 

2 


_ 1 

2 














hi 








1 

14 


1 
6 


0.0510 


-0.1403 


-0.0113 


0.0054 


k^i 








1 

14 


1 
6 


0.0984 


0.1214 


0.0156 


0.0026 


hi 














-0.0397 


0.0117 


-0.0025 


0.0304 


hi 














0.0335 


0.0239 


-0.0462 


-0.0112 


hi 


626126 
272277 


56281 
51730 


3 
7 


1 

14 


-0.6494 


-0.0380 


-0.0186 


-0.0057 


hi 


313063 
363036 











-0.9135 


0.0873 


-0.0571 


0.0209 


hu 








1 


-1 














hi 








2 
3 


1 

3 














hi 








2 

63 


1 

27 


-0.0659 


0.0595 


-0.0218 


0.0335 


h 4i 








1 

21 


1 
9 


0.0237 


-0.0173 


-0.1336 


-0.0316 


hi 








1 

126 


1 

108 


0.0094 


-0.100 


0.0010 


-0.0017 


hi 








1 

84 


1 

36 


0.0108 


0.0163 


0.0103 


0.0023 


e* 


4661194 
816831 


8516 
2217 








-1.9043 


-0.1008 


0.1216, 


0.0183 


fi 


-17.3023 


8.5027 


4.5508 


0.7519 


2.0040 


0.7476 


-0.5358 


0.0914 


9i 


14.8088 


-10.809 


-0.8740 


0.4218 


-2.9347 


0.3971 


0.1600 


0.0225 



Table A.l: The "magic numbers" for the Wilson coefficients in b — > S7. The a, are the 
eigenvalues of 7^ divided by 2/3 - 



Appendix B 

Operator Basis Transformations 



Here we address the transformation relations of the operator bases ( p. 50 - 2. 57 ) and ( [2.58 - 



|2.65| ). In general, if we change the operator basis from Q to P = WQ, to translate the 
scheme-dependent NLO quantities like the NLO Wilson coefficients or the NLO anomalous 
dimension matrix, we need the finite parts of the one-gluon exchange corrections to the 
operator matrix elements where the scheme dependence is manifest. If the one- loop matrix 
elements in the starting basis are 



(Q) = (i + ^ r )<Q>W (B.l) 



then the ones in the transformed basis can be written as 

a .±A /P\W 



(P) = \l + ^)(P)V» (B.2) 

= W (l + ^Ar) <Q> + 0{a 2 s ) (B.3) 

with 

Ar = W~VW - r (B.4) 

Denoting the Wilson coefficients in the transformed basis with Z, we can derive the 
transformation rules starting from the fact that the full amplitude has to be the same in 
both bases 

C T (Q) = Z T (P) (B.5) 

which leads to 

Z=(W-Y(l-^r T )c + 0(al) (B.6) 

For the sake of completeness we also give the transformation rule for the anomalous 
dimension matrix 7' in the new basis which can be derived for example via differentiation 
of the left and right hand side of ( |B.6| ) with respect to In /i using the definition of the 
anomalous dimension matrix in ( |2.38| ). 

l' = W {7 + ^[Ar, 7] + 2/3 (^) 2 Ar} W' 1 + 0(a 3 s ) (B.7) 

Here /3q was defined in (|2.16|) and the square brackets denote the commutator Ar 7 — 7 A 



r. 



The operator basis (|2 . 5 0KP~57|) we use mostly is apart from an overall factor of 4 and 



the fact that we reversed the numbering of the current- current operators Qi^ the same 
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B. Operator Basis Transformations 



as the one of Greub, Hurth, and Wyler [S0j. Chetyrkin, Misiak, and Miinz [4T| and also 
Buras, Czarnecki, Misiak, and Urban however, use another operator basis ( p. 58 



2.65|) , which is better suited for multi-loop calculations. The transformation matrix W in 



four space-time dimensions is given by 

/ 



W 
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1 
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1 
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1 





1 

















1 

6 


1 

2 


1 
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1 

2 














16 





4 

















8 
3 


8 


2 
3 


2 


























1 


























1 



(B.8) 



and transforms our operator basis ( |2.5U| - |2"37| ) into the CMM basis ( |2.58| - |2TBo] 



D=A 



(B.9) 



k=i 



Note the block diagonal structure of W. The most important consequence hereof is that 
W~ x as well has block diagonal form with the inverse of the submatrices as entries. 

To determine Ar we have to insert all the operators into current-current and penguin 
diagrams with one additional gluon. For the 6x6 submatrix of Ar concerning the current- 
current and QCD penguin operators we calculated the diagrams (b), (c), (d), and (e) in 
Fig. Ol with insertions of Q\...% and P L .. 6 . We get for N = 3 



Ar 



6x6 



/ 


128 
27 


52 
9 


7 
81 


7 

27 


7 
81 


7 

27 


\ 




2311 
81 


173 
27 


35 
243 


35 
81 


35 
243 


35 
81 












178 
27 


34 
9 


164 
27 


20 
9 












1- L 

1 9 


-25+/ 
3 


-18-/ 
9 


6 + 1 












160 
27 


16 

9 


146 
27 


2 
9 




V 








Z + 9 


u 3 


3 + 1 


-11-/ 
3 


/ 



(B.10) 



which indeed transforms the NLO anomalous dimension matrix for our operator basis 
( |A.10| ) into the 6x6 submatrix of ( |A.9| ). Without calculation we already can say that the 
first six entries of the seventh and eighth row of Ar vanish because of the block diagonal 
structure of W and W~ l and the fact that the magnetic penguin operators do not mix 
into the four-quark operators. Yet, determining the entries in the seventh and eighth 
columns of r and r' would amount to calculating the finite parts of the 0(a s ) virtual 
corrections to the matrix elements for b — ► 57. This challenging two-loop calculation is 
what Greub, Hurth, and Wyler pDl did for Qi 2 and (^7,8 of our operator basis. Recently 



Buras, Czarnecki, Misiak, and Urban performed this calculation for the complete CMM 
basis ( |2.58| - |276l)| ) pT,|B2|. But the NLO matrix elements of the penguin operators in our 
operator basis have not yet been calculated so that we cannot determine the complete Ar. 



Therefore we lack the information to translate the results of [pl] , |62| into our operator basis. 
The use of only one operator basis would perhaps be interesting from an esthetical point 
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of view, but it is not necessary. We use the two different operator bases for physically 
distinct processes in B — > Vj, namely hard vertex and hard spectator corrections which 
is perfectly alright if the corresponding Wilson coefficients are used. For the Wilson 
coefficient of the electromagnetic penguin operator we need the next-to-leading value, 



which was calculated in |4T| in the CMM basis. 

To get C-j and C§ no transformation is necessary because the electro- and chromo- 
magnetic penguin operators are the same in both bases. Furthermore Q-j and Q 8 cannot 
be expressed as linear combinations of Q1...6 or Pi. ..6 and they also do not mix into these 
operators. Using eq. ( |B.5| ) we therefore can identify the coefficients of (Q7,s) and {Pr t s) 
to be the same. So we can directly use the NLO result for Z7 from |^l| for our operator 
basis. 

Chetyrkin, Misiak, and Miinz use a different strategy to derive the transformation 
properties between the 6x6 submatrices of the NLO anomalous dimension matrices |42| . 
They pass from one operator basis to the other by a series of subsequent redefinitions of 
physical and evanescent operators. 
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Appendix C 

Weak Annihilation in B — > V7 



In this appendix we explicitely prove to 0(a s ) that weak annihilation contributions to 
B — ► V7 decays are calculable within QCD factorization. To this end we have to show 
that all collinear and soft divergences cancel to leading power when one-gluon exchange 
to Fig. [4.5| is considered. By leading power we here mean the lowest nonvanishing order in 
the power expansion, where the entire annihilation contributions start only at subleading 
power (cf. section [O ). Our proof here is analogous to the one given in subsection |3.2.4 . 
We again expect the soft divergences to cancel because of the colour transparency argu- 
ment and the collinear ones via collinear Ward identities. 

Let us first show these cancellations for the weak annihilation involving operators 
Q1...4 and the case where the photon is emitted from the light quark in the B meson 
(see Fig. |C1| ). For the soft region all components of the gluon momentum s become 
small simultaneously: s ~ A. Power counting shows that every single of the six possible 
diagrams from Fig. (J.l is logarithmically IR divergent. The sum of the two contributions 
in the second row of Fig. |C.1| , however, vanishes when the equations of motion for the 
vector meson constituent quarks q and q' are used: 



j a (v ft- fir* 1 r»(-vjt+,/)Y 

s 2 — 2v s ■ k 



Q 



2vk a F» 



s 2 — 2v s ■ k 
2vk a Y^ ~ 



—2v s ■ k —2v s ■ k 
0(A) 



q + 0(X) 



(C.l) 



where we abbreviated T M = 7 M (1 — 75). Hence, the soft divergences cancel to leading 
power in the total amplitude. 

For s collinear with the light-cone momentum k of the vector meson, i.e. s = ak + . . ., 
we have the scaling 

s + ~ A , s± ~ A, 8~ ~ A 2 (C.2) 



so that 



d A s ~ A 



s-k~ A 2 , 



s 2 ~ 



A 



(C.3) 



Using these scaling laws we find that the divergence of the single diagrams in the collinear 
region is also logarithmic. We examine again the sum of the two contributions in the 
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C. Weak Annihilation in B — > V'j 








Figure C.l: One-gluon exchange corrections to the annihilation contribution to B — > V^y 
decay from Q1...4 operator insertions. The six possible diagrams are obtained through 
combining the three diagrams in the first row with each of the two in the second row. 



second row of Fig. C.l: 



q 

= 4 

oc k a 



s 2 — 2v s ■ k 
2(v - a)k a Y^ 
s 2 — 2v s ■ k 



s 2 — 2v s ■ k 
2{y - a)k a Y^' 
s 2 — 2v s ■ k 



(C.4) 



The k a is contracted with the 7 a from the gluon vertices in the first row of Fig. |C.l 



k a 



_(s 2 - 21 ■ s)(s 2 - 21 ■ s + 2s • q - 21 ■ q) 



21 ■ q(s 2 - 21 ■ s + 2s ■ q - 21 ■ q) 



+ 



21 ■ q(s 2 - 2p b ■ s) 



(C.5) 



The above expression vanishes when the equations of motion for the B meson quarks are 
used. So there are no collinear divergences at leading power either. 

The other annihilation contribution which appears at the same level of power counting 
comes from the insertion of operators Q^^ and photon emission from one of the final state 
vector meson quarks. The one-gluon-exchange diagrams for the case where the photon is 
emitted from the antiquark in the vector meson are depicted in Fig. |C.2| . Both constituent 
quarks in the vector meson are light, collinear, and carry a momentum of the same size. 
Photon emission from the other vector meson quark therefore can be treated completely 
analogously. 

A power counting as above shows that the two diagrams where the gluon hits the 
off-shell light antiquark in the vector meson (third diagram in second row of Fig. |C.2j 
combined with the two diagrams in the first row) are infrared finite both in the soft and 
collinear region. All other combinations are logarithmically divergent. 
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Figure C.2: One-gluon exchange corrections to the annihilation contribution to B — > V'y 
decay from Q 5 e operator insertions. The six possible diagrams are obtained through 
combining both diagrams in the first row with each of the three in the second row. Com- 
binations with the third diagram in the second row are infrared finite. Photon emission 
from the other vector meson quark is not shown. 



For the soft gluon the sum of the first two contributions in the second row of Fig. |C.2| 
again cancels to leading-power when the equations of motion for the vector meson quarks 
q and q' are used: 



T(-vft+ 4- () 4(-vft+ 4) la 
s 2 — 2vs ■ k + 2vk ■ q — 2s ■ q)(s 2 — 2vs ■ k) 
2vk a T{-v 2vk a T{-v fi- q) 4 



+ 



4v 2 k ■ qs ■ k 



4vvk ■ qs ■ k 



la{v ft- 4)^{-v ft- £ 
(s 2 — 2vs ■ k)2vk ■ q 

■ + 0{\) 



0{\) 



(C.6) 



Here Y represents the Dirac structure of the quark currents of Q^fi. 

Similarly, in the collinear region the sum of the same two diagrams is again propor- 
tional to k a , which is contracted with the 7 a from the gluon vertices in the first row of 



there 



Fig. |C.2| . For on-shell B meson quarks this once more vanishes. 

This completes our proof. For all annihilation contributions of Q\,„&, to 0(a s 
are no collinear or soft infrared divergences at leading non- vanishing power in l/m&. 
Annihilation contributions therefore are, although power suppressed, calculable in QCD 
factorization. 
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